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ABSTRACT 


This report presents the results of the second year of a three year 
design project on the automation of the Environmental Control and 
Life Support System (ECLSS) of the Space Station Freedom (SSF) . The 
results are applicable to other space missions that require long 
duration space habitats. The report begins with a description of 
conceptual controls which are developed for the Water Recovery and 
Management (WRM) subassembly. Mathematical modeling of the Air 
Revitalization (AR) subassembly is presented next. The report 
concludes the work done by the Kansas State University NASA/USRA 
interdisciplinary student design team with a discussion of the 
expert system which was developed for the AR subassembly. 
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1.0 INTRODUCTION 


1 



1.1 PROJECT DESCRIPTION 


The development of Environmental Control and Life Support Systems 
(ECLSS) for the Space Station Freedom, future colonization of the 
moon and Mars missions presents new challenges for present 
technologies. ECLS Systems which operate during long duration 
missions must be semi-autonomous to allow crew members 
environmental control without constant supervision. A control 
system for the ECLSS must address these issues as well as being 
reliable. The Kansas State University Advanced Design Team is in 
the process of researching and designing controls for the 
automation of the ECLS system for Space Station Freedom and beyond. 

The ECLS system for Freedom is composed of six subsystems. The 
Temperature and Humidity Control (THC) subsystem maintains the 
cabin temperature and humidity at a comfortable level. The 
Atmosphere Control and Supply (ACS) subsystem insures proper cabin 
pressure and partial pressures of oxygen and nitrogen. To protect 
the space station from fire damage, the Fire Detection and 
Suppression (FDS) subsystem provides fire sensing alarms and 
extinguisher. The Waste Management (WM) subsystem compacts solid 
wastes for return to earth, and collects urine for water recovery. 
The Atmosphere Revitalization (AR) subsystem removes C0 2 and other 
dangerous contaminants from the air. The Water Recovery and 
Management (WRM) subsystem collects and filters condensate from the 
cabin to replenish potable water supplies, and processes urine and 
other waste waters to replenish hygiene water supplies. 

These subsystems are not fully automated at this time. 
Furthermore, the control of these subsystems is not presently 
integrated; they are largely independent of one another. A fully 
integrated and automated ECLS system would increase astronaut ' s 
productivity and contribute to their safety and comfort. 
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1.2 THREE PHASE DESIGN PLAN 


A three phase approach was implemented by the Kansas State 
University Advanced Design Team to design controls for the ECLSS. 
The first phase, completed within one year, researched the ECLSS as 
a whole system and then focused on the automation of the Atmosphere 
Revilatization (AR) subsystem. 

During the second phase, the system control process was applied to 
the Atmosphere Revitalization subsystem. To aid in the development 
of automatic controls for each subsystem and the overall ECLSS, 
mathematical models have been developed for system simulation on a 
computer. Expert system control as well as conventional control 
methods are being tested on the models. Using the Atmosphere 
Revitalization subsystem control system as a "proof of concept" , 
the other ECLSS subsystems will be automated. 

Finally, during phase three, the control system of the six 
subsystems will be combined to form a control system for ECLSS. 
The expert system developed for the AR will be expanded to control 
the ECLSS as well as provide fault diagnosis and isolation to the 
astronauts . 

The Kansas State University Design Team has completed phase one and 
two. Mathematical models for the C0 2 Removal Assembly, C0 2 
Reduction Assembly and Oxygen Generation Assembly, as well as an 
expert system, have been developed for the AR. The mathematical 
models are written in the C-Language. At this time, the models 
function independently at the assembly level. The expert system, 
using a written shell program called CLIPS, monitors and controls 
the AR System assemblies in a hierarchical manner. 
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1.3 DESIGN TEAM DESCRIPTION 


The Kansas State University Advanced Design Team is composed of 
engineering students from several disciplines, a student from 
general science and education, two graduate student assistants, and 
engineering faculty members. Architectural, chemical, computer, 
electrical, and mechanical engineering disciplines are represented 
by both students and faculty. 
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1.4 ORGANIZATION OF THE DOCUMENT 


Chapter two of the report presents the work which was done by the 
conceptual controls group during the first semester of the design 
project. Conceptual controls were developed for the Water Recovery 
and Management System. These control schemes will be tested next 
year when the appropriate mathematical models are developed, and 
will aid in the development of an integrated ECLSS control scheme. 

The mathematical models which were formulated for the Air 
Revitalization (AR) are presented in chapter three. These models 
are being used to test the conceptual control schemes which were 
developed in the first year of the project and also the expert 
system which was developed this year. 

An expert system was created which monitors the AR system. Chapter 
four details the development and use of the expert system within 
the AR. 
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2.0 CONCEPTUAL CONTROLS 
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2.1 INTRODUCTION 


This chapter outlines the progress made by the conceptual 
controls group toward the design of a control system for the 
Water Recovery and Management (WRM) System of the ECLSS. The WRM 
System was divided into four subsystems which are 1) Water 
Quality Monitor and Process Control Water Quality Monitor, 2) 
Urine Processor, 3) Hygiene Water Processor and 4) Potable Water 
Processor. 

In designing the conceptual control scheme of the WRM System, 
several steps were taken before the final conceptual control 
scheme was developed. The first step was to research existing 
control schemes. Using this information, a tentative control 
scheme for each subsystem was developed. The final step involved 
critiqueing and revising each control scheme by group members and 
the faculty advisor until the best conceptual control scheme 
emerged . 

Each section of this Chapter contains a description of the 
current control scheme, the proposed control scheme, the system 
monitoring scheme and a detailed sensor layout. A comparison 
between the current and proposed control schemes is also 
included. 


7 



2.2 WATER QUALITY MONITOR 
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2.2.1 WATER QUALITY MONITORING SYSTEMS DESCRIPTION 


Two separate systems monitor the water quality from the various 
processors. The systems are called the Process Control Water 
Quality Monitor and the Water Quality Monitor which are discussed 
below in that order. 


PROCESS CONTROL WATER QUALITY MONITOR 

The Process Control Water Quality Monitor (PCWQM) is a system of 
monitors that continuously monitor the quality of product water 
from a particular processor — either a potable processor, a 
hygiene processor, or a urine processor. Each processor has a 
dedicated PCWQM through which all of its product water flows. 

The PCWQM has two distinct sections — the Continuous Section and 
the Intermittent Section (See Figure 2.2-1). 

In the Continuous Section, product water from the particular 
processor flows into the PCWQM and on out through a series of 
four monitors. The monitors measure pH, Temperature, Iodine, and 
Conductivity levels every 15 seconds [5:1], [8:6]. Conductivity 
and pH are measured using standard conductivity and pH cells. 
Iodine is measured using absorption of UV/VIS light at 476 
nanometers [5:1]. It is not known how the temperature is 
measured. The temperature reading is used to adjust the 
conductivity and pH readings. 

In the Intermittent Section, Total Organic Carbon (TOC) is 
analyzed and pH and TOC are calibrated. Every 15 minutes, a 10 
ml sample is diverted from the product water flow in the 
Continuous Section to be analyzed in the Intermittent Section 
[5:1-2], [8:6]. 

The TOC is measured using an ultraviolet reactor. The water 
sample is mixed with dilute phosphoric acid from the Acid 
Injection Module (AIM) . The acid reacts with any inorganic 
carbon present. Any gases present in the water sample or 
resulting from the reaction of acid and inorganic carbon are 
vented. Oxygen is mixed with the water sample and ultraviolet 
light oxidizes the organic carbon to carbon dioxide. The 
infrared detector measures the amount of infrared light (4.28 
micrometers) absorbed by the carbon dioxide and correlates this 
measurement to TOC content [5:1]. The TOC analyzer is calibrated 
every 24 hours by an unknown process using fluid from the 
Calibration Injection Module (CIM) [8:6], 

The pH Calibration Monitor calibrates the pH monitor in the 
Continuous Section. The pH Calibration Monitor measures the pH 
of the 10 ml sample and compares its reading to that of the pH 
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monitor. If the readings are different, the controlling 
processor calibrates the pH monitor accordingly. The pH 
Calibration Monitor itself is calibrated every 24 hours using 
fluid of known pH from the Base Injection Module (BIM) [8:6]. 
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Figure 2.2-1 Current PCWQM 

(adapted from [2:1]) 


The acid-tainted 10 ml sample and the used calibration fluid are 
sent to the hygiene waste water storage. The gasses from the 
gas/ liquid separator are vented to the cabin. This minuscule 
amount of gas is eventually removed by the Trace Contaminant 
Control System. 

Every 15 seconds, the PCWQM sends the pH, conductivity, raw 
iodine, and adjusted iodine readings to the Data Management 
System (DMS) [5:2], [8:6]. The raw iodine is the iodine 
measurement at the pH of the water. The adjusted iodine is the 
iodine measurement that would be present at a pH of 7 [5:2]. 

TOC measurements are sent to the DMS as they become available. 
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WATER QUALITY MONITOR 


A Water Quality Monitor (WQM) analyzes each full tank of product 
water from a Potable Water Processor (PWP) and Hygiene Water 
Processor (HWP) pair. Only one product water tank from either 
the PWP or the HWP can be analyzed at any one time. The WQM also 
accepts manual samples when it is not already conducting an 
analysis. See Figure 2.2-2 for system layout. 

The Data Management System (DMS) sends a signal to start the 
analysis process. The process begins by circulating water to the 
WQM from a selected product water tank. Circulation continues 
until a representative sample of the processed water can be 
extracted for analysis. This small sample is distributed to the 
various analyzers using an unknown internal sample distribution 
system. Each analyzer performs its own tests. As results from 
each analyzer become available, they are sent to the DMS. Some 
analyses are available almost instantaneously, while other 
analyses take up to 12 hours [7:2]. 

The various analyzers and the corresponding measured parameters 
for the WQM are listed in Table 2.2-1. In addition to these, the 
Gas Chromatograph/Mass Spectrometer can identify thousands of 
other compounds by sending the Mass Spectrometer scan to earth 
for further analyses. 
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Figure 2.2-2 Current WQM 

( adapted from [3:1]) 


Table 2.2-1 WQM Measured Parameters 

(based on information from [7:1-2] and [3:1]) 


ANALYZER 

PARAMETER 

UV-Visible 

Spectroscopy 

pH Electrode 

Conductivity Cell 

Spectrophotometer 
(UV-Vis . Spectroscopy) 


BLOCK NAME MEASURED 

(In Fig. 1.1-2) 


Turbidimeter 

unknown 

unknown 

Spectrophotometer 


Turbidity 

pH 

Conductivity 

Color 

Iodine 


Gas Analyzer 
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continued 

Table 2.2-1 WQM Measured Parameters 


MEASURED 


ANALYZER BLOCK NAME 

PARAMETER 


Gas Chromatograph/ 
Alcohol 

Mass Spectrometer 


Gold Film Analyzer 

UV Oxidation 

High Performance 
Liquid Chromatograph 


GC/MS Aliphatic 

Carboxylic Acid 
Halogenated 
Hydrocarbons 

Mercury Analyzer Mercury 

Total Organic Carbon TOC 

HPLC Ammonium 

Barium 

Cadmium 

Calcium 

Copper 

Iron 

Lead 

Manganese 

Magnesium 

Nickel 

Potassium 

Silver 

Arsenic 

Chloride 

Chromium 

Cyanide 

Fluoride 

Iodide 

Nitrate 

Selenium 

Sulfate 

Sulfide 

Phenol Alcohol 


Bact. Sensor Assembly Sensor Assembly Bacteria 

(Bacteria Concentrator, 

Incubation Chambers, 

Conductivity Sensor, 

Chromogenic Limulus, 

Reagents & Reaction 
Chamber, & Colorimeter) 
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2.2.2 WATER QUALITY MONITORING CONTROLS 


Two separate systems — the Process Control Water Quality Monitor 
(PCQWM) and the Water Quality Monitor (WQM) — analyze the water 
quality of the product water from the potable, hygiene, and urine 
water processors. A PCWQM constantly monitors the quality of the 
continuous flow of product water from its associated potable, 
hygiene, or urine processor. A WQM thoroughly analyzes one 
sample at a time from each processor storage tank as the final 
test before the water is used. A WQM is responsible for 
monitoring a potable water processor and hygiene water processor 
pair. The PCWQM and the WQM will be discussed separately below. 


PCWQM CONCEPTUAL CONTROLS 
CURRENT PCWQM CONTROL SCHEME 

The current PCWQM control scheme is a three operational mode 
system that has a sensor system consisting only of valve position 
indicators monitoring each valve. 

The three operational modes are ON, STANDBY, and OFF [4:1]. In 
the ON mode, all systems are powered. In the STANDBY mode, only 
the electronic control system is powered. In the OFF mode, all 
electrical power to the PCWQM is off. The PCWQM operational mode 
is determined by the operational mode of its associated 
processor. For example, if the associated processor changes from 
the STANDBY mode to the ON mode, the processor sends a signal 
telling the PCWQM to do the same. 

When the PCWQM is in the ON mode and product water is available, 
the PCWQM monitors the product water. If the monitored values 
are within the normal range, the results are sent to the Data 
Management System (DMS) . If the analyzed value in a particular 
monitor is out-of -range, the calibration of the monitor is 
checked. If the calibration is acceptable, the results are sent 
to the DMS. If the calibration is unacceptable the astronauts 
are notified. 

Valve position sensors provide feedback to determine if the 
valves are functioning properly. Readings from the sensors are 
monitored by the DMS. 


PROPOSED PCWQM CONTROL SCHEME 

Several features where added to the initial control scheme to 
allow for greater fault detection and isolation (Fig. 2.2-3). 
The features added are: 
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■ calibrating analyzers that are out of calibration before 
notifying the astronauts. 

■ adding seven sensors. 

■ adding a manual shut-off option to the solenoid valves on 
the Calibration Injection Module (CIM) , the Base Injection 
Module (BIM) , and the Acid Injection Module (AIM) . 


aow 

IN 


aow 

our 



Figure 2.2-3 Proposed PCWQM 

(adapted from [2:1]) 


PCWQM MONITORING SCHEME 

System components are listed below along with a description of 
the monitoring scheme for each component. 

Monitors 


The pH and Total Organic Carbon (TOC) monitors are calibrated in 
the TOC section of the PCWQM. The iodine, conductivity, and 
temperature monitors do not need calibration. 


15 




Valves 


Position indicators are used to determine if a valve is in the 
correct position. Whenever power is not available to the valve, 
the valve moves to a fail-safe position. 

Check Valves 


Flow sensors are placed after these valves to detect any leaks. 
CIM. BIM and AIM 

Fluid level sensors are used to indicate if a fluid tank is 
almost empty. Liquid volume sensors are used to ensure that the 
proper amount of fluid is drawn from the tank. 


COMPARISON OF PCWQM CONTROL SCHEMES 

The proposed control scheme has several advantages. It has more 
feedback because of the additional sensors. This increases 
knowledge that the system is working properly. The automatic 
calibration of out-of-range analyzers saves the astronauts from 
constantly monitoring and adjusting them. Finally, the valves 
with a manual override for the CIM, BIM, and AIM can be used to 
prevent damage to the PCWQM. 

The disadvantages of the proposed system are the increased cost 
and complexity, the small increases in power consumption, and the 
increase in associated microprocessor size. 


WQM CONCEPTUAL CONTROLS 


CURRENT WQM CONTROL SCHEME 

The current WQM control scheme is a three operational mode system 
with five sensors, in addition to valve position sensors, to 
monitor the process. 

The three operational modes are ON, STANDBY, and OFF [6:1]. In 
the ON mode, all subsystems and analyzers are powered. In the 
STANDBY mode, only the ion pump, embedded controller, and related 
electronics are powered. In the OFF mode, all electrical power 
to the WQM is off. The normal WQM operational mode is the 
STANDBY mode, but the WQM changes to the ON mode when a sample 
becomes available. The OFF mode is used for maintenance of the 
WQM. 

When the WQM is testing a sample (ON mode) , each particular test 
result is sent to the DMS as it becomes available, as long as the 
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result is in the normal range. If it is not in the normal range, 
the calibration status of the particular analyzer is determined. 
If calibration is acceptable, the test results are sent to the 
DMS. If the calibration is unacceptable, the astronauts are 
notified. 

Sensors provide feedback to determine if various system 
components are functioning properly. Readings from the sensors 
are monitored by the DMS. 


PROPOSED WQM CONTROL SCHEME 

This scheme modifies the current control scheme by: 

■ calibrating analyzers that are out of calibration before 
notifying the astronauts. 

■ adding ten sensors. 

■ stirring the particular processor storage tank before a 
sample is sent to the WQM. 

■ discontinuing any remaining WQM tests if one test is out of 
range . 

See Figure 2.2-4 for the proposed WQM system layout. 


WQM MONITORING SCHEME 

System components are listed below along with a description of 
the monitoring scheme for each component. 

Analyzers 

Each analyzer block has its own calibration mechanism within the 
block. 

Valves 


Position indicators are used to determine if a valve is in the 
correct position. Whenever power is not available to the valve, 
the valve moves to a fail-safe position. 

Check Valves 


Flow sensors are placed after these valves to detect any leaks. 




Figure 2.2-4 Proposed WQM 


Pumps 

Pressure sensors are placed directly before and after each pump 
to make sure the pump has the proper pressure differential across 
it. 

Also, current sensors are placed on each pump to determine if the 
pump motor is drawing the proper amount of current. This will 
help determine if a pipe is plugged or if a pump fin is broken. 

Pressure Regulator 

A pressure sensor is placed after the regulator to determine if 
the pressure is being regulated properly. 
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COMPARISON OF WQM CONTROL SCHEMES 

The proposed control scheme has several advantages. It has more 
feedback because of the additional sensors. This increases 
certainty that the system is working properly. The automatic 
calibration of an analyzer before notifying the astronauts saves 
the astronauts' valuable time by not calling them unless 
absolutely necessary. Stirring the processor storage tanks 
before sending a sample to the WQM ensures that a representative 
sample is analyzed. Stopping the remaining WQM tests if one test 
is out-of-range saves power. 

The disadvantages of the proposed system are the increased cost 
and complexity, the increased power consumption for some of the 
changes, and the increase in associated microprocessor size. 
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2.3 URINE PROCESSOR 
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2.3.1 URINE PROCESSOR ASSEMBLY DESCRIPTION 


The main objective of the Urine Processor Assembly (UPA) is to 
recover water from the urine and flush water from the 
urinal. [2:1] The UPA works in conjunction with the hygiene water 
processor by sending product water to be further processed and 
with the potable water processor by accepting brine to be 
processed. 

Two urine processors will be located on the Space Station 
Freedom. A schematic of the UPA is shown in Figure 2.3-1. [4:1 of 
1] One will be in the laboratory module and another in the 
habitat module. Each UPA can accept urine and waste water 
partially pretreated with Oxone (manufactured by Dupont) in the 
urinal on a 24 hour basis. [2:1] To further inhibit microbial 
growth and ammonia formation, a small amount of concentrated 
sulfuric acid will be injected by the acid injector pump into 
each flush from the urinal as it proceeds to the pretreated urine 
collection tank . [2:1] 

The UPA units must be in a preheated state before a batch process 
is initiated. Electric heaters will operate for a sufficient 
period of time before each batch process, to achieve the 
necessary temperature for processing. 

Processing will begin when the pretreated storage tank is 60-85% 
full depending on whether the other processor is running. [3:Fig 
I] Only one processor can run at a time because of power 
constraints. A UPA will also begin processing if it has been 
more than 48 hours since the last batch process. Pretreated 
urine should not be stored longer than 48 hours. [3: Fig I] 

The pretreated urine is drawn from the storage tanks through a 
regenerative heat exchanger which preheats the urine and flush 
water. This mixture is then sent through a filter to trap any 
large particles in the waste water. The heated waste water is 
pumped into the Thermoelectric Integrated Membrane Evaporative 
System (TIMES) . It first enters the TIMES heat exchangers which 
receive heat transferred from the Thermoelectric Devices 
(TED). [2: 2] Heated waste water is pumped through bundles of 
hollow fiber membrane tubules whose external surfaces are exposed 
to a saturated gaseous environment maintained at a pressure of 2 
psig.[2:2] Relatively clean product water permeates through the 
membranes and evaporates from their outer surface due to the high 
temperature of the waste water and the reduced pressure outside 
the tubes. [2:2] Accompanied by the evaporation is a decrease in 
temperature and an increase of the brine solute in the recycle 
flow. [2: 2] The water vapor enters a condenser where it is cooled 
and partially condensed, releasing latent heat to the TED on 
either side of the condenser. The TED'S act like heat pumps by 
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Waste Gas Vent 


transferring the latent heat from the cooler condenser to the 
warmer waste water heat exchanger . [2:2] 



Figure 2.3-1 Urine Processor Assembly Schematic 













Due to the waste heat generated by the electrical input to the 
TED's, only part of the steam condenses and a mixture of vapor and 
liquid leaves the condenser [2:2]. The mixture exits the TIMES and 
more vapor is condensed as the mixture passes through the 
regenerative heat exchanger and an air-cooled heat exchanger .[ 2 : 2 ] 
The mixture of water vapor and liquid water enters a rotary water 
separator where water is removed from the stream and the 
noncondensible gases are vented to the Trace Contaminant Control 
System (TCCS).[2:2] A compressor removes enough gas to maintain a 
pressure of approximately 2 psig at the water separator gas 
outlet. [2:2] The water leaving the separator is checked by the 
Process Control Water Quality Monitor (PCWQM) for conductivity, pH, 
total organic carbon (TOC) and microbes. A listing of water 
quality requirements is shown in Table 2 . 3-1. [2 :Table 1] 


Table 2.3-1. 

Product Water Quality Requirements 

conductivity 

< 

150 micro-mho/ cm 

pH 


3-8 

total organic 

carbon (TOC) < 

50 ppm 

microbes 

< 

1 cfu/100 ml 


If the product water does not meet the PCWQM water quality 
requirements it remains in the recycle loop and is reprocessed. At 
about 25% dissolved solids the flow rate through the hollow fiber 
membrane tubules decreases below a minimum set point and the 
decreased flow rate initiates a brine dump to the brine storage 
tank. [2:1] Once the brine tank is full it is replaced with an 
empty one and the full tank is sent back to earth. Pretreated 
urine will displace the old recycle loop brine when processing 
begins again. [2:1] 
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2.3.2 URINE PROCESSOR ASSEMBLY CONTROLS 


CURRENT URINE PROCESSOR ASSEMBLY CONTROL SCHEME 

The current control scheme used for the Urine Processor Assembly 
(UPA) subsystem is a batch process or timed process depending on 
the situation. When the pretreated urine tank is 60-85% full or 
the time elapsed since the last processing is 48 hours then 
processing will begin. Pretreated urine should not be stored 
longer than 48 hours. See Figure 2.3-2 for a schematic of the 
current system. [4:1 of 1] 


PROPOSED URINE PROCESSOR ASSEMBLY CONTROL SCHEME 

A timed control scheme was considered for controlling the UPA, but 
the efficiency of the equipment would be decreased if small batches 
were run on timed intervals. The batch process was found to be the 
best method for controlling the UPA. A few additions were made to 
improve the system and these are explained below. 

An Emergency Process mode was added to the existing modes of On, 
Off, Standby, and Emergency Shutdown .[ 5 : 1 ] The Emergency Process 
mode would be initiated when there is a need for waste water to be 
processed immediately because of another system failure. This mode 
would be included in the Standby mode, so if directed the UPA will 
begin processing even though the pretreat tank is not 60-85% full. 

Another system improvement would be to add a pH sensor to the inlet 
of the UPA subsystem. Currently there is a small quantity of 
sulfuric acid being injected into the urine and waste water stream 
as it enters the waste water bus and flows into the pretreated 
urine tank. Since the pH of human urine varies from 4.8 to 8.4 
depending on changes in the body metabolism or diet, the urine may 
need to be treated with more or less sulfuric acid than the 
measured amount. A pH sensor located in the waste bus between the 
2 -way valve and the pretreat urine storage tank could indicate how 
much acid pretreat is required for that particular flush. This 
would decrease the bacterial growth and thus decrease the amount 
of recycling required through the UPA. 

When the recycle water has reached 25% solute, the efficiency of 
the water permeating through the hollow fiber membrane tubules 
decreases and a brine dump is initiated. There are no sensors in 
the current system to indicate when the recycle water has reached 
25% solute. A flow sensor was considered since the flow through 
the hollow fiber membrane tubules decreases as the solute 
increases. Measuring the conductivity of the recycling water would 
indicate when there is 25% solute in the recycling water. Once 
this set point is reached a brine dump would be initiated to the 
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Figure 2.3-2 Current Urine Processor Assembly Control Schematic 
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brine storage tank. Other recommendations consist of placing 
additional sensors for system monitoring and failure detection. 

The present UPA subsystem consists of a wide range of sensors. 
However, a few additional sensors could help track down specific 
problems that may occur. The recommended additions along with the 
current sensors are discussed in the following sections. 


UPA MONITORING SCHEME 

The following information is a list of all the major system 
components. Sensor methods currently being employed and suggested 
are included in these descriptions. For each sensor there will be 
a comparison made with a range of the desired output and an error 
signal will be generated. The controlling processor unit will 
determine if the sensor is giving correct output or if a failure 
has occurred. In the event of failure, the astronauts will be 
alerted. See Figure 2.3-3 for a schematic of the proposed system 
with the additional sensor placements. [4 : 1 of 1] 

Acid Pretreat Tank 

A quantity sensor is located in the tank to indicate when the 
sulfuric acid is getting low. 

Pretreat Urine Tank 


A quantity sensor is located in the tank to indicate when to begin 
processing. A liquid sensor detects any pretreated urine leaks. 

Heat Exchangers 

There were no sensors placed near the heat exchangers located 
external to the Thermal Integrated Membrane Evaporative System 
(TIMES) . It is suggested that temperature sensors be placed at the 
outputs and inputs of the regenerative heat exchanger to monitor 
any temperature deviations that would indicate that the heat 
exchanger was not working properly. Flow sensors were placed in 
the air-cooled heat exchanger to monitor any changes in the flow of 
air. 

Pumps 

The two pumps used in the UPA are the acid injector pump and the 
recycle loop pump. Each one has a speed sensor and a current 
sensor. A differential pressure sensor was included with the 
recycle loop pump and it is suggested that a differential pressure 
sensor be added to the acid injector pump. 
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TIMES 


The TIMES includes evaporators, heat exchangers, thermoelectric 
devices, and a condenser. There are currently three temperature 
sensors located on each evaporator to monitor the efficiency. 

A temperature sensor was added to the output of each heat exchanger 
to check the efficiency of the thermal electric devices to transfer 
latent heat to the recycling waste water. Flow sensors were added 
to the inlets of the TIMES to monitor any obstructions that may 
occur in the pipes. 

Air /Water Separator 

A speed sensor is located on the air/water separator to monitor the 
compressor speed. In the product water line there is a pressure 
sensor which monitors the separator liquid pressure. A relief 
valve is included to relieve any excess pressure buildup. In the 
waste gas line there is a liquid sensor to sense any leaks in the 
separator. A pressure sensor to monitor the separator gas pressure 
is also in the waste gas line. 

Compressor 

A speed sensor is located on the compressor to determine if it is 
running . 

Valves 


All of the valves, excluding the 2-way valves located between the 
filters of the acid pretreat tank and brine tank, have manual back- 
ups. Motor driven valves are indicated by a semi-circle setting on 
top of the valve in Figure 2.3-3. If the valve is not motor driven 
then it is a solenoid valve which is electrically actuated. All 
the 2-way and 3-way valves have position indicators. 

Check Valves 


The check valves prevent backflow from contaminating other lines. 
PCWOM 

The Process Control Water Quality Monitor is explained in section 
2.2 in this chapter . 

Flow Sensors 


It is suggested that flow sensors be placed in select locations 
along the recycle loop to indicate any obstructions in the piping. 
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COMPARISON OF CONTROL SCHEMES 


The basic control scheme of the current control system and the 
proposed control system are the same. The batch process will 
continue to be the controlling system of the UPA subsystem. The 
proposed pH sensor in the waste water bus line would allow a more 
accurate amount of acid pretreat to be added. This would decrease 
the amount of bacterial growth and recycling necessary through the 
processor. The additional Emergency Process mode would be 
advantageous in emergency situations. Measuring the conductivity 
of the recycling water to initiate a brine dump would also be 
highly effective. Also, the recommended sensor additions can allow 
for more fault protection, detection, and isolation. Although the 
additional sensors of the proposed system could increase the cost 
of the system, it would be more productive and efficient. 
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2.4 HYGIENE WATER PROCESSOR 
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2.4.1 HYGIENE WATER PROCESSOR DESCRIPTION 


The Hygiene Water Processor (HWP) consists of four major sections: 
gas pressurization, waste water storage, water processing and 
product water storage. 

The Gas Pressurization section pressurizes the air on the air side 
of the bellows tanks (water busses) . This pressure is necessary in 
order for the water to flow through the processor. 

Water enters the HWP at a pressure of 34.5-82.74 kN/m 2 (5-12 psig) . 
This pressure range allows the waste water storage section to 
maintain a waste water bus pressure between 34.5-82.74 kN/m 2 (5-15 
psig) during both static and flow conditions. The gas 
pressurization section maintains this pressure by monitoring the 
water inlet pressure and by opening a valve on the gas side of the 
tank to lower pressure if the pressure exceeds 82.74 kN/m 2 (12 
psig) and opening a valve from the higher pressure section if the 
pressure falls below 34.5 kN/m 2 (5 psig). If the waste water 
storage tank level exceeds 80% capacity. The inlet valves will be 
closed to prevent potential damage to the water storage tanks. A 
check valve is placed before the inlet valve on the waste water 
storage tank. This prevents the delivery of water from the waste 
water storage tank back to the waste water bus in the event of a 
waste water bus depressurization due to external leakage .[ 1 : 1 ] 

When water is available in the waste water storage section, and the 
product water tanks are not full, the processing section will 
process and deliver water from the waste water storage section to 
the product water storage section. A water pump draws the waste 
water from the waste water storage tank to the water processing 
section. In this water processing section the water goes through 
a sterilization loop that sterilizes and filters the waste water. 
Sterilization is accomplished by heating and maintaining the water 
at a temperature greater than 121° C (250° F) for 20 minutes. For 
energy efficiency, a regenerative heat exchanger is used in 
conjunction with the particulate filter electrical heater to 
achieve the sterilization temperature. After heating, the water 
goes into the combination filter-residence tank, which uses a 40 
micron filter and provides for a residence time greater than 20 
minutes. Water leaves the filter and exits from the sterilization 
section by flowing through the hot side of the regenerative heat 
exchanger. The pressure of the cold inlet and cold outlet water 
sides of the regenerative heat exchanger are monitored to detect an 
expended filter and insure proper pressure maintenance to maintain 
water in the liquid state. [5: 4] 

The water then enters the recycle loop which circulates the water 
through two membrane-based filtration modules. The first membrane 
filtration is done in the Ultrafiltration (UF) module. Suspended 
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solids and dissolved macromodules (i.e. skin cells and soaps) are 
rejected while lower molecular-weight species pass through the 
membrane and become feed for the second stage of membrane 
filtration. The Reverse Osmosis (RO) module, is capable of 
reversing this feed flow within the UF module periodically to 
backflush the module and prolong membrane life. A valve also 
directs the outflow from the UF module to the brine processor to 
flush out the concentrated brine solution and replace it with fresh 
waste water. The reservoir tank placed before the recycle loop 
reduces the brine concentration in the recycle loop. This will 
allow an increase in the time between backf lushes. [1: 2] 

The permeate water from the UF module passes through an air/water 
heat exchanger to cool the water to a temperature below 45° C (113° 
F) and prevent the delivery of water with a temperature in excess 
of 45° C (113° F) to the RO module. The permeate water then enters 
the second membrane-based filtration module, the RO module. The RO 
module removes low-molecular weight species (i.e. salts and low 
molecular weight organic material) which passed through the UF 
module, while allowing water to pass through and return to the 
recycle flow upstream of the reservoir tank. A check valve 
prevents backflow of fresh waste water through the RO module. In 
this way the RO module separates a relatively cleaner (compared to 
the UF permeate) water, which is then delivered to the unibeds. 
The RO permeate water is monitored by both a pressure sensor and a 
conductivity sensor. The pressure sensor will determine if flow 
blockage downstream of the RO has occurred. The conductivity 
sensor is used to determine when to backflush the recycle loop. A 
check valve prevents the backflow of water from the unibeds into 
the RO module. This backflow could damage the RO module. 

The RO permeate then flows to the post-treatment loop. The post- 
treat loop consists of two unibeds that use the principle of 
adsorption to remove contaminants. Adsorption is a process in 
which contaminant particles are deposited on the surface of solid 
particles. The solid particles known as adsorbents, are chosen to 
preferentially adsorb the contaminants expected in the waste 
stream. The amount and type of material adsorbed varies from 
adsorbent to adsorbent. The unibeds in the HWP are arranged in a 
series configuration, such that when one bed is expended, the next 
bed is used. A new bed is then placed in the spot vacated by the 
expended bed. These beds will be replaced in staggered intervals 
of 30 days. [5:4] 

The Process Control Water Quality Monitor evaluates the quality of 
the processed water and determines if it meets the required 
specifications. If it does not, the water must be redirected to 
the water pump inlet to be reprocessed. The PCWQM is discussed in 
detail elsewhere in this report. 
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The product water storage section stores the processed water. If 
the product water tanks are full, the inlet valves are closed and 
the HWP enters a standby state. The only other time the inlet 
valves are closed is during sterilization. If the product water 
tanks are empty, the outlet valves are closed. Otherwise, the 
outlet valves are kept open, and the water bus pressure is 
controlled between 124.1-206.9 kN/m 2 (18-30 psig) by the 
pressurization section. This pressure range allows the product 
water storage section to maintain a product water bus pressure of 
between 102.4-206.9 kN/m 2 (15-30 psig) during both static and flow 
conditions. The pressurization section maintains this pressure by 
monitoring the water pressure and by operating a valve to ambient 
pressure if this pressure exceeds 206.9 kN/m 2 (30 psig) and by 
turning on the compressor if this pressure falls below 124.1 kN/m 2 
(18 psig) . The product water tanks can accept processed water and 
at the same time deliver stored water to the hygiene water 
bus .[1:3] 
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2.4.2 HYGIENE WATER PROCESSOR CONTROLS 


CURRENT HYGIENE WATER PROCESSOR CONTROL SCHEME 

The present system uses a continuous cycling of waste water. The 
water processing section will process water when the waste water 
storage tank is full and there is room in the product water tanks 
for the processed water. The mode of operation is determined from 
a computer used by the astronauts. The water is sterilized and 
filtered in the processing section. All valves and major system 
components are controlled by electrical signals that arise from 
sensors placed throughout the system. The hierarchy of this 
control scheme is yet to be done. Detection of failure in the 
current scheme is limited to pressure, temperature and conductivity 
sensors and the PCWQM system. Again, the hierarchy of this control 
scheme is for future work [3:1-2]. 


PROPOSED HYGIENE WATER PROCESSOR CONTROL SCHEME 

This control scheme will be based on the amount of water to be 
processed and the amount of water in the product water storage 
tanks. The waste water will be processed when the tank is 60% full 
and the product water tanks are only 60% full. The mode of 
operation will be determined by the ratio of the amount of water to 
be processed to the amount of water in the product water storage 
tank, thus allowing the system to operate independently of an 
external command. The modes of operation are standby, on, off and 
emergency shutdown . [1:4] 

Hygiene Water Processor Standby Mode 

When the waste water storage tank is less than 60% full or the 
product water tank is full, the processing section will be in 
standby mode. This mode is used for equipment preheat, tank 
pressurization and setting the valves to the ready position. 

Hygiene Water Processor On Mode 

When the waste water storage tank is over 60% full and the product 
water tank is less than 60% full the processing section will begin 
processing waste water. The processing occurs as discussed in the 
operational description of the HWP. A submode of the on mode is 
backflush. During backflush normal processing is stopped and the 
filters are backf lushed to clean and extend their useful life. 
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Hygiene Water Processor Off Mode 

9 

When the waste water storage tank is less than 3 0% full and the 
product water tanks are greater than 20% full the processing 
section will go into the off mode. During off mode processes are 
gradually stopped and the power is removed from the system. While 
the system is in this mode the system is isolated and can be 
serviced. 

Hygiene Water Processor Emergency Shutdown Mode 

When either a part of the hygiene water processing system breaks or 
the power consumed by this system is needed elsewhere the processor 
will go into emergency shutdown mode. This mode is similar to off 
mode except the system is stopped as quickly as possible. 


SYSTEM MONITORING SCHEME 

The major components and the sensors associated with each component 
are listed in the following paragraphs. The monitoring scheme will 
allow the system to be controlled, and allow component failures to 
be detected. The output of the sensor will be compared with a 
desired output and an error signal will be generated if the output 
is not within the determined limits. A microprocessor will monitor 
and determine, through the use of an artificial intelligence 
program, if a component failure has occurred. The astronauts will 
be alerted via a message displayed on a computer terminal in the 
event of a failure. The necessary steps to correct any failure can 
then be taken. [4:1-12] Figure 2.4-1 shows a schematic of the 
proposed system . [1:4] 

Gas Pressurization Section 


Pressure sensors are placed on the waste water inlet line to 
determine if the water pressure is between 34.5-82.74 kN/m 2 (5-12 
psig) . If the pressure is either high or low then air is either 
released or added to the gas side of the storage tank. Quantity 
sensors are located on the water busses (bellows tanks) to 
determine if they are pressurized. Quantity sensors are used to 
measure the amount of gas or liquid contained in a tank. 

Solenoid Valves 


Position indicators are located on each valve. The position of the 
valve is compared to the position required during the process. If 
a valve becomes stuck open or shut the system will alert the 
astronauts and begin emergency shutdown mode. 
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Figure 2.4-1 Proposed Hygiene Processor Schematic 
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Figure 2.4-2 Legend of the Proposed Hygiene Water Processor 


Check Valves 

A flow sensor will be placed on either side of a check valve to 
determine if the valve is leaking. Check valves assure that water 
will not flow backwards in the system. Backwards flow could 
possibly damage the system. 



A quantity sensor is located on the storage tank to determine if 
the processing cycle should begin, the system should go to standby 
mode or if the system should proceed to off mode. The quantity 
sensor also determines if a solenoid valves placed on the inlet and 
outlet sides of the tank should be open or closed. These valves 
either allow the water to be processed or stored. 
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Water Pump 


Pressure sensors are placed on either side of the pump to determine 
if it is working properly. Values for the pressure on either side 
of the pump are checked with predetermined values in the 
microprocessor program. If these values are not within the 
predetermined range the system will go into emergency shutdown 
mode. 

Heat Exchanger 

Pressure sensors are placed on the cold inlet and outlet sides of 
the heat exchanger to detect an expended filter and ensure the 
water is at the proper pressure to remain in the liquid state. 

Filter /Residence Sterilization Tank 

Temperature sensors are located in the tank to ensure that the 
water has reached its sterilization temperature of 121.1° C (250° 
F) . A current sensor is also located on the power supply to the 
heater to determine if it is working properly. 

Ultrafiltration Module 


Pressure sensors are located on either side of the UF module to 
determine if the filter is expended. Values for pressure still 
need to be determined. A flow sensor is also located before the UF 
module to determine if it is being backf lushed or not. 

Heat Exchanger 

Flow sensors are placed on both the inlet and outlet sides of the 
heat exchanger to determine if it has become clogged. 


Reverse Osmosis Module 


Pressure sensors are located on either side of the RO module to 
determine if the module is functioning properly and determine if 
blockage has occurred. A conductivity sensor will be used to 
determine when to flush the recycle loop. 

Unibeds 

A conductivity sensor is placed between the two unibeds to 
determine if they are expended. Each unibed will last 
approximately 60 days according to current data. 
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Process Control Water Quality Monitor 


This is a separate subsystem of the WRM and is discussed at length 
in section 2.2 of this report. 

Product Water Storage Tank 

A quantity sensor is placed on the product water storage tank to 
determine if the system should begin processing, go to standby 
mode, or switch to off mode. The quantity sensor determines the 
position of the solenoid valve on outlet of the tank. If the tank 
becomes empty these valves would be closed. 

Gas Pressurization Section 


Pressure sensors are located on the exit side of the product water 
storage tank to determine if the pressure is between 34.5-82.74 
kN/m 2 (15-30 psig) . Quantity sensors are located on the water 
busses to determine if they are pressurized. 

COMPARISON OF CONTROL SCHEMES 

The new control system has advantages over the current control 
scheme. The primary advantage is that the new control system has 
feedback to determine more accurately the condition of the system. 
It operates on an as-needed basis, processing water only when there 
is a sufficient amount to be processed and when water is needed. 
This enables the new system to operate more efficiently by 
consuming less power. The new system also frees up the astronauts 
to perform their duties and not be bothered with controlling the 
system. This control scheme also contains many fail safe features. 
This allows the system to shut itself down if a problem arises that 
would endanger the astronaut's safety. 

The proposed control scheme has three main disadvantages; it will 
be more complex, consume more space and be more costly to 
implement. The new method also requires further research into the 
proper operating pressures of each point and a complex hierarchy of 
control. [3:1-2] 
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2.5 POTABLE WATER PROCESSOR 
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2.5.1 POTABLE WATER PROCESSOR DESCRIPTION 


Figure 2.5-1 [2:8] is a diagram of the Potable Water Processor 
(PWP) . The purpose of this subsystem is to produce potable water 
for consumption by space station astronauts. Condensate from the 
Temperature and Humidity Control (THC) and Atmosphere 
Revitalization (AR) subsystems provide feed water for the PWP. The 
Potable Water Distribution Bus (PWDB) distributes processed water 
about the space station. Two processors, one in the Habitation 
Module and one in the Laboratory Module, comprise the PWP. Each 
has a capacity based on 2.84-5.16 kg H 2 0/man-day with each 
processor operating for forty hours [3:1]. Two twenty hour 
processing periods exist in the forty hour cycle. The break allows 
the Condensate Storage Tank to refill should the Condensate Storage 
Tank empty before the end of the forty hour cycle. 

During a normal processing cycle, condensate from the Condensate 
Storage Tank enters the Filter/Heater at a nominal rate of 1.81 
kg/hr [3:1]. The Condensate Storage Tank and Main Pump bring system 
pressure to 103.4-206.8 kPa[3:l]. The Filter/Heater raises the 
water from ambient temperature to 121°C. The water remains at this 
temperature for at least 20 minutes[3 : 1] . The Heat Exchanger pre- 
heats water entering the Filter/Heater and allows the effluent to 
move to the Unibeds at ambient temperature. The Filter/Heater 
destroys microorganisms and filters particulate. 

The Unibeds are multifiltration beds composed of ion-exchange 
resins and activated carbon. They remove inorganic and high 
molecular weight organic compounds. They also maintain a residual 
iodine level of 2-4 ppm[2:3]. The Unibeds have a 30 day life- 
time^:!]. Replacing the filters involves rotating the second 
filter into the lead position and placing a new filter in the 
second position. The pressure and conductivity sensors monitor the 
Unibeds for clogging. The controlling CPU can automatically 
backflush the beds, based on bed age and sensor information. The 
Urine/Flush Water Processor receives backflush brine. 

The water moves on to a Volatile Removal Apparatus (VRA Reactor) . 
Hydrogen peroxide is the oxidant for a catalytic oxidation reaction 
in the VRA Reactor. Alcohols and low molecular weight organic 
compounds are oxidized for removal by a set of Polishing Ion- 
Exchange beds. The oxidation reaction takes place at 121°C or 
higher [3: 2]. A regenerative Heat Exchanger pre-heats water 
entering the VRA Reactor and brings VRA Reactor effluent to ambient 
temperature. A Phase Separator removes oxidation process gas by- 
products and vents them to the Gas Rack. A gas check valve 
prevents backflow from the Gas Rack to the PWP. The Phase 
Separator is assumed to be passive and does not need monitoring. 
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Figure 2.5-1 Multifiltration Potable Subsystem Schematic 




Polishing Ion-Exchange beds follow the Phase Separator and remove 
the alcohols and low molecular weight organic compounds oxidized by 
the VRA Reactor. Pressure and conductivity sensors monitor for 
filter clogging. The controlling CPU can backflush these filters 
automatically with the Urine/ Flush Water Processor receiving the 
brine . 

The Process Control Water Quality Monitor monitors water quality 
during processing. Extensive batch testing by the Water Quality 
Monitor provides calibration information to the controlling CPU 
about the PCWQM. The PCWQM samples every 15 seconds to find pH, 
conductivity, temperature and total organic carbon content [Table 
2.2-1]. These factors determine if the processed water will enter 
the Potable Water Holding Tanks or if the water should be recycled 
through the PWP. To recycle the water, a three-way solenoid valve 
redirects flow to the Main Pump inlet. 

On its way to the Main Pump, the water passes through a Microbial 
Check Valve (MCV) that at present adds iodine or another biocide to 
the water. The MCV is assumed to be passive and lacks monitoring. 
A check valve on the MCV outlet prevents water flowing into the PWP 
from the Condensate storage Tank from entering the MCV. 

Each PWP has four Potable Water Holding Tanks. For each half cycle 
(20 hours) they have different functions. One tank is filling with 
processed water for later astronaut use while a second tank is 
filling with processed water for batch testing by the WQM. A third 
tank fills with potable water for standby and emergencies and the 
fourth tank supplies potable water to the astronauts. The tank 
supplying the astronauts is the tank that was filling for standby 
in the previous half cycle. These tanks have bellows to maintain 
an outlet pressure of 103.4-206.8 kPa on the PWDB. Control of 
bellows pressure is similar to that used on the Condensate Storage 
Tank. The two-way solenoid valves attached to the inlet and 
outlets of the Potable Water Holding Tanks have fail-safe positions 
reachable by spring force. They will close in the event of power 
failure or power removal. Manual movement of the valves is 
possible. 

Each PWP processor has a controlling CPU. The controlling CPU 
makes decisions on processor mode, start and stop times and 
emergency control based on sensor information. The controlling CPU 
reports to a master CPU, which then communicates with the space 
station Systems Software. Providing each PWP with its own 
controller relieves the Systems Software of processor control. 
This allows a more sophisticated control scheme for the processor, 
capable of near autonomous operation. 
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2.5.2 POTABLE WATER PROCESSOR CONTROLS 


CURRENT POTABLE WATER PROCESSOR CONTROL SCHEME 

The current Potable Water Processor (PWP) has five modes of 
operation with sensors to report on the system status. The five 
modes of operation are: Initialization, Standby, Operate, Shutdown 
and Stop[3:3]. Space Station Freedom will have two PWPs, one in 
the habitation module and one in the laboratory module. They 
operate one at a time. Each processor operates for forty hours. 
The operating processor will continue to process condensate stored 
in a tank until the tank is empty or the Systems Software signals 
a stop. 


PROPOSED POTABLE WATER PROCESSOR CONTROL SCHEME 

The proposed PWP control scheme uses additional sensors for 
monitoring of storage tank pressure and three-way solenoid valves 
to allow backf lushing of individual filters. Figure 2.5-2 shows 
placement of sensors and control valves. In the proposed system, 
each PWP processor will be controlled by a dedicated CPU with a 
third CPU monitoring both processors and reporting to the Systems 
Software. Changes and additions to the processor modes provide 
greater control of the system during normal operation and emergency 
situations. The proposed modes of operation are: Initialization, 
Update, Standby, Process, Normal Stop, Emergency Stop, Backf lush, 
Manual Valve Control, Receive, and Off. 

The Initialization mode is the primary start-up mode for the PWP. 
The controlling electronics of the PWP and the heating element of 
the filter/heater receive power. If pumps that can be engaged and 
disengaged while maintaining a constant motor speed are present, 
then they will be powered. Bellows Storage Tank Compressors are 
powered. The valve positions will not be changed during this mode. 

Update mode will usually follow Initialization and the Controlling 
CPU will find all valve positions to prepare for moving into 
Process, Backf lush or Manual Valve Operation modes. The CPU also 
checks fluid levels in the Condensate Storage Tank, Potable Water 
Storage Tanks, and the H 2 0 2 Bladder Tank. If there is no condensate 
to operate on, or if the Potable Water Storage Tanks are full, then 
the system will move to the Off and Receive modes to conserve 
power . 

In the Standby mode, all electronics are powered, as are heaters, 
pumps and compressors. The Controlling CPU will know all valve 
positions and fluid levels in condensate storage, potable water 
storage and H 2 0 2 bladder tanks. Standby is a staging point for the 
processor to move to Process or Backf lush modes. The processor 


44 




45 


Figure 2.5-2 Proposed Potable Water Processor Control Schematic 



will not remain in this mode for longer than it takes the 
Controlling CPU to change valve positions for the next mode of 
operation. 

The PWP will spend most of its forty hour duty cycle in the Process 
mode. During this mode, the PWP is actively processing Condensate 
Storage Tank water. The Condensate Storage Tank will not receive 
condensate from the Atmosphere Revitalization (AR) or Temperature 
and Humidity Control (THC) sub-systems while Process mode is 
active. The pressure of the Condensate Storage Tank is monitored 
and maintained such that the outlet pressure is sufficient for 
passing water through the system. When the bellows pressure is 
within the proper range, the outlet valve will open and the main 
pump will engage to begin processing. The temperature of the 
filter/heater is monitored and adjusted as necessary. Line 
pressures, water conductivity, quality testing PCWQM and holding 
tank levels are monitored. 

The Normal Stop mode of the PWP will be selected at the end of the 
processor's forty hour cycle, if all water from the Condensate 
Storage Tank has been processed, or if the Potable Water Holding 
Tanks are approaching capacity. The Controlling CPU will close the 
Condensate Storage Tank outlet valve but continue to power the Main 
Pump until the Potable Water Storage Tank quantity sensors report 
no addition of water to the receiving tank(s) . The three-way valve 
following the PCWQM is then set to recycle water back to the main 
pump inlet and power will be removed from pumps, heaters, and 
compressors. If the processor must be restarted after the Normal 
Stop mode, the Initialization and Update modes must be activated. 
The Off mode will generally follow the Normal Stop mode. 

Emergency Stop provides a way to stop flow from the processor to 
the Potable Water Storage Tanks quickly. Upon receiving the signal 
to move to the Emergency Stop mode, the Controlling CPU will change 
the three-way valve after the PCWQM to recycle the processed water 
back to the main pump inlet. The outlet valve of the Condensate 
Storage Tank also will be closed. Emergency Stop removes power 
from the pumps, heaters and compressors. The Emergency Stop is 
nearly the reverse of the Normal Stop, except that water will be 
trapped in the system. Most of the trapped water can be removed by 
the Backflush mode or by operating the PWP in the Manual Valve 
Operation mode. To begin processing water after an Emergency Stop, 
the Initialization and Update modes must be used. 

The Manual Valve Operation mode allows the astronauts to start and 
stop pumps and compressors manually, to have manual control of 
heater temperature and manual control of valve positions. The 
valve positions can be changed at the valve or via the Systems 
Software, Master CPU and dedicated CPU. This mode allows the 
astronauts complete control over the system during repairs or 
routine maintenance (i.e., changing of filters). This mode can 
only be reached after Initialization. If the system fails and 
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cannot be initialized, then the Systems Software must inform the 
Master CPU to use only the undamaged PWP processor. 

The Condensate Storage Tank cannot receive water from the AR and 
THC sub-systems during Update, Process, Normal Stop, or Emergency 
Stop modes. The Receive mode is a background task for the 
Initialization, Standby, Backflush and Off modes. When the 
processor is in the Receive mode, the bellows pressure is reduced 
below the inlet line pressure (but no lower than ambient pressure) 
and the inlet valve to the Condensate Storage Tank will be opened 
allowing condensate to flow into the storage tank. Process mode 
requires repressurization of the bellows. 

The Off mode removes power from the electronics and generally will 
be accompanied by the Receive mode. The Off mode follows the 
Normal Stop and Emergency Stop modes. In an emergency or loss of 
system control, the Off mode can be executed by the astronauts via 
Systems Software. Executing the Off mode out of sequence (i.e., 
because of astronaut intervention) , removes all power from the 
system. All solenoid valves in the PWP have fail-safe positions 
(inhibit flow into and out of system). In the event of abnormal 
power removal, these fail-safe positions are reached using springs. 
The Off mode cannot be executed by the Controlling CPU except 
following a Normal Stop mode or Emergency Stop mode. 

The Controlling CPU will monitor sensors by multiplexing them onto 
fewer lines. If a sensor begins to send information requiring the 
attention of the CPU (i.e., pump speed dropping below a set level) 
the CPU will poll that sensor more frequently until the problem has 
been resolved. If a problem is not resolvable by the Controlling 
CPU, it will report to the master PWP CPU, which will decide 
whether to switch processing to the dormant PWP, continue 
processing with the ailing PWP, halt all processing, signal the 
astronauts via the Systems Software or some combination of the 
above . 


MODES FOR PROPOSED CONTROL SCHEME 

Potable Water Processor Power-up Mode 

The normal Power-up Sequence begins with the master PWP CPU 
selecting the Habitation or Laboratory PWP. Processor condition 
(damaged or healthy) and last processor to operate determine PWP 
activation. If both processors are healthy, the dormant processor 
will begin a power-up sequence as the operating processor is 
following a normal stop sequence. The mode order for a normal 
power-up sequence is: RECEIVE (background mode, operating since 
last normal stop sequence) -> INITIALIZATION -> UPDATE -> STANDBY - 
> PROCESS. If the filters in the processor are nearing the end of 
their lifetime (approximately 30 days) , the normal power-up 
sequence could be modified to contain a Backflush mode before 
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Initialization. The Backflush could be performed on each section, 
or on separate sections. The modified mode sequence is: RECEIVE -> 
BACKFLUSH -> INITIALIZATION -> UPDATE -> STANDBY -> STANDBY -> 
PROCESS . 

Potable Water Processor Power-down Mode 


The normal power-down sequence execution occurs at the end of a 
processor's forty hour processing cycle, if the Condensate Storage 
Tank is empty or if the Potable Water Storage Tanks are reaching 
capacity. The mode order for a normal power-down sequence is: 
PROCESS (current mode) -> NORMAL STOP -> OFF -> RECEIVE. If the 
filters are near the end of their lifetime and a backflush was not 
performed during the normal power-up sequence, a Backflush mode may 
be inserted following the Normal Stop mode. The modified normal 
power-down sequence is: PROCESS -> NORMAL STOP -> BACKFLUSH -> OFF 
-> RECEIVE. 

Potable Water Processor Emergency Power-down Mode 

In an emergency (i.e., power must be diverted to another space 
station function or there is an uncorrectable failure in the PWP) 
the Controlling CPU can initiate an emergency power-down sequence. 
The sequence of modes is: PROCESS (current mode) -> EMERGENCY STOP 
-> OFF. No Backflush or Receive mode is part of the emergency 
power-down sequence. The master PWP CPU also will alert Systems 
Software to the emergency power-down and provide a reason for the 
power-down . 


SYSTEM MONITORING SCHEME 

Water Storage Tanks (Condensate and Product) 

A pressure sensor monitors the bellows pressure of each tank. If 
the bellows pressure is not in range for the current mode of 
operation, two solenoid valves provide a means for altering the 
pressure. One solenoid valve connects the bellows to a compressor 
to increase pressure. The fail-safe position for this valve is 
closed. The second valve can open to vent air to release pressure. 
The Gas Rack receives vented air. The fail-safe position for this 
valve is open. Level sensors monitor the quantity of water in the 
tanks. The Controlling CPU places this value in memory so that as 
the water level in the tank begins to decrease rapidly (i.e., leak) 
the CPU can change bellows pressure accordingly and, if necessary, 
initiate the Emergency Power-down Sequence. 

Condensate Bus Inlet 


A pressure sensor on the Condensate Bus Inlet monitors the water 
pressure so that the Condensate Storage Tank bellows pressure is 
lower than the bus pressure. If the bellows pressure has been 
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reduced to ambient pressure but the bus pressure is still lower, 
the Controlling CPU will alert the Master CPU. The Master CPU will 
then attempt to increase the bus pressure. A check valve placed 
between the two-way solenoid valve and the Condensate Storage Tank 
prevents backflow from the Condensate Storage Tank to the bus. 
Integral flow sensors in the check valve monitor for valve failure. 

Condensate Storage Tank Outlet 

The Condensate Storage Tank Outlet has a pressure sensor to monitor 
outlet pressure. Information from this sensor is used to adjust 
the bellows pressure. A two-way solenoid valve controls outlet of 
water from the Condensate Storage Tank. A check valve follows the 
two-way solenoid to prevent backflow into the storage tank. A 
second check valve appears on the Recycle loop of the PWP to 
prevent flow of water from the Condensate Storage Tank from 
entering the outlet of the recycle loop. The check valves have 
integral flow sensors to monitor for valve failure. 

Pumps 

Both the Main Pump and the H 2 0 2 Pump have speed sensors. These 
sensors monitor for proper operation of the pumps and provide an 
estimate of the pump outlet pressure. The estimated outlet 
pressure provides a check on both the pump performance and the 
pressure sensors following the pump. Three-way solenoid valves 
surround the main pump. Forcing potable water backward through the 
pump to the Urine Processor Assembly backf lushes the pump. 

Filter /Heater and 1* Heat Exchanger 

A set of pressure sensors surrounds the Filter /Heater and Heat 
Exchanger. These monitor input and output pressures to alert the 
Controlling CPU to leaks or the need for a backf lush of the 
Filter/Heater . Three-way solenoid valves surround the 
Filter/Heater for backf lushing. A temperature sensor in the 
Filter/Heater provides information for maintaining the 250° 
temperature required. A mechanical thermostat backs up the 
temperature sensor. The heating element has a current sensor to 
monitor power consumed. 

Multifiltration and Polishing Unibeds 

Conductivity sensors surround the two Multifiltration and two 
Polishing Unibeds followed by a pressure sensor. These sensors 
inform the Controlling CPU about the condition of the Unibeds. A 
marked increase in water conductivity from previous sampling 
suggests Unibed clogging. Each Unibed can be backf lushed 
individually. 
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Volatile Removal Apparatus (VRA) Reactor and 2 nd Heat Exchanger 

A set of pressure sensors surrounds the VRA Reactor and Heat 
Exchanger. The Bladder Tank has a controllable pump to release 
known amounts of H 2 0 2 into the VRA Reactor. A level sensor on the 
HjOj Bladder Tank and speed sensor on the pump inform the 
Controlling CPU about H^ release. The CPU releases to match 

the flow of water into the VRA Reactor, as monitored by a flow 
sensor. Pressure sensors surround the VRA Reactor and Heat 
Exchanger to monitor for proper line pressure. Insufficient line 
pressure suggests a clogged VRA Reactor. Three-way solenoid valves 
surround the VRA Reactor for backf lushing. 

Phase Separator 

The Phase Separator vents gasses to the Gas Rack. The Phase 
Separator is assumed to be a passive system and currently lacks a 
control scheme. 

Process Control Water Quality Monitor 

This is a separate subsystem of the PWP and is discussed at length 
in section 2.2 of this chapter. 

Potable Water Storage Tanks 

During potable water processing, each of the four tanks has a 
different function. One tank is filling with processed water for 
later distribution, another is filling with processed water for 
further testing by the Water Quality Monitor subsystem, another 
tank is on standby (empty) and the fourth tank is delivering water 
to the Potable Water Distribution Bus. Each tank has a level 
sensor to monitor fluid level inside the bellows. Bellows pressure 
is monitored and used to control the inlet two-way solenoid valve 
and vent two-way solenoid valve. The same compressor that supplies 
air to the Condensate Storage Tank is used to supply the Potable 
Water storage Tank bellows with air. The tank designated as the 
Potable Water Supply Tank maintains a higher pressure to deliver 
water to the Potable Water Distribution Bus at 103.43 to 206.86 kPa 
gage. A pressure sensor on the distribution bus is also used to 
control the Potable Water Supply Tank bellows pressure. The two 
filling tanks maintain a pressure below line pressure (not to fall 
below ambient pressure) to allow filling. Check valves on the 
outlet lines of the Potable Water Storage Tanks prevent backflow 
from the distribution bus into the tanks. The check valves contain 
integral flow sensors to monitor for valve failure. A two-way 
solenoid valve on the distribution bus can be closed to isolate the 
storage tanks from the bus. The valve's fail-safe position is 
closed (isolate tanks) . There are two of these valves present in 
the event of device failure. 
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Microbial Check Valve 


The Microbial Check Valve (MCV) is used to add an additional 
biocide to the processed water as it is recycled to the main pump 
inlet. A check valve after the MCV prevents flow of water from the 
Condensate Storage Tank into the MCV. The MCV is assumed to be a 
passive system and lacks a control scheme. 

Check Valves 


All check valves in the PWP have integral flow sensors (not shown 
on the system diagram) . These flow sensors notify the Controlling 
CPU of a valve failure. In the event of a pressure sensor failure, 
the check valve flow sensors indicate high line pressure. 
Depending on the sensor location, number of sensors reporting out 
of range, water waiting to be processed, Potable Water Storage Tank 
levels and time remaining in the processing cycle, the Controlling 
CPU can opt to continue processing or initiate the Emergency Stop 
Sequence. The Master CPU will notify Systems Software in the event 
of Emergency Stop Sequence. 

Two and Three-Wav Solenoid Valves 

All solenoid valves in the PWP contain position sensors. These 
sensors report to the Controlling CPU the current valve position. 
If a valve is stuck, the Controlling CPU will initiate the 
Emergency Stop Sequence. The Master CPU will notify Systems 
Software of the Emergency Stop Sequence. 


COMPARISON OF CONTROL SCHEMES 

The proposed control scheme places the PWP under the control of a 
dedicated CPU which reports to a Master CPU. The Master CPU then 
reports to the Systems Software and ultimately the astronauts. By 
placing these "buffers” between the PWP and the astronauts, the PWP 
can 1) be self-contained, 2) self -repairing (in a limited way) and 
3) portable. Because this control scheme is portable, the same PWP 
can be implemented on different craft, each with different Systems 
Software. The software of the Master CPU can be rewritten to 
interface with the host Systems Software, rather than modifying the 
host Systems Software and running the risk of losing control of 
another vital system. The current control scheme implies that the 
same software that is monitoring other life support systems also 
will be required to monitor the PWP. 

Because the PWP will have limited self-repair capabilities, it can 
continue to run and ensure an emergency supply of potable water for 
the astronauts. The inclusion of a Manual Valve Operation mode 
makes the testing and repair of the system easier for the 
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astronauts. With this mode, they can perform dynamic testing on 
the system. The current control scheme only allows for static 
testing. 

A self-contained system eases repair and upgrading of the system. 
Addition of sensors requires only the Controlling CPU software to 
be modified rather than the Master CPU software and/or the Systems 
Software. 

Drawbacks to this system result from increased system complexity. 
The proposed control scheme has ten modes of operation, twice the 
current scheme. The addition of sensors to storage tank bellows 
and check valves will increase system cost. The writing of control 
software for Controlling CPUs and Master CPUs also will increase 
the cost. The number of sensors will produce a system that is 
impossible to test exhaustively. Computer simulation of the system 
must supplement physical tests. Reporting insufficient information 
to the Controlling CPU<, the Master CPU or the Systems Software 
could result in a system failure and notification of the astronauts 
too late to correct the failure. 


52 



3.0 MATHEMATICAL MODELING 
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3.1 INTRODUCTION 


The modeling of the Air Revitalization (AR) system is divided 
into three distinct areas of concentration. First is the Carbon 
Dioxide Removal Assembly. Next, the Carbon Dioxide Reduction 
Assembly, and finally the Oxygen Generation Assembly (OGA) . 

Carbon Dioxide Removal Assembly uses a series of molecular sieve 
beds to extract C0 2 from the cabin air of the Space Station. The 
C0 2 is removed from the air and is pumped to a C0 2 storage tank. 

Stored C0 2 is drawn from the storage tank into the Carbon Dioxide 
Reduction Assembly for conversion to methane (CH 4 ) and water 
vapor. The methane/ water vapor mixture is cooled and water is 
removed by a centrifugal separator. The methane is placed in a 
storage tank for later disposal. 

The Oxygen Generation Assembly converts water to oxygen and 
hydrogen gas in electrolysis cells. An electrical current passes 
through electrodes in each cell, causing the "cracking" of water 
into hydrogen and oxygen gas at the electrode surface. 

Math modeling serves two purposes. The first purpose is to 
provide a means of learning how each system operates. Second it 
provides the Expert System with data for its knowledge base. 

The following subsections give an analysis of the three Air 
Revitalization subassemblies. 
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3.2 C0 2 REMOVAL ASSEMBLY MODEL 


DESCRIPTION 

The C0 2 Removal Assembly is part of the Air Revitalization (AR) 
Subsystem. Its purpose is to remove C0 2 from the cabin 
atmosphere, deliver C0 2 to the C0 2 Reduction Assembly and return 
humidified air to the cabin. This is done using a four bed 
molecular sieve consisting of: two desiccant beds to remove water 
vapor from incoming air, two C0 2 adsorption beds, a blower to 
force air through the system, a C0 2 pump, a C0 2 accumulator, a 
pre-cooler and five multiple-flow selector valves [1:8]. Figure 
3.2-1 illustrates these major components. 



Figure 3.2-1 co 2 Removal Assembly 


HATH MODEL 

Assumptions 

1. All beds modeled as lumped systems. 

2. For simplicity, we have neglected the fact that the 
adsorbing/ desorbing processes vary as a function of 
distance through the bed. 

3. Assumed thermal equilibrium which negates the dependence 
on bed length. 

4. Thermal equilibrium assumed for the C0 2 desorbent bed. 

5. Pump operation is 100% efficient, isentropic, and 
adiabatic. 
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6 . 


The model can be improved by sub-dividing each bed into 
many smaller beds. This method results in plug flow 
operation in which the simulation takes place several 
times for each bed. 

7. To implement the improved modeling scheme, the overall 
bed volume and sorbent mass will be divided by the 
number of plugs. 

All four beds have been modeled as lumped systems. For 
simplicity, we have neglected the fact that the adsorbing/ 
desorbing processes vary as a function of distance through the 
bed. In our model, we have assumed instantaneous thermal 
equilibrium, which in effect, negates the dependence on bed 
length. We intend to improve the model by sub-dividing each bed 
into many smaller beds. Plug flow operation results, in which 
the simulation takes place several times for each bed. To 
implement this scheme, the overall bed volume and sorbent mass 
will be divided by the number of "plugs". The simulation will 
most likely be modified using 2-D arrays. 

The relationship equations and state variable equations are as 
follows: 

Equations 

1) For the C0 2 desorbent bed assuming thermal equilibrium, 
one obtains 
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( 6 ) 
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The required definitions are given as: 

P b = C0 2 equilibrium pressure of bed (kPa) , 
ki = constant (picked to be 0 < K, < 1) , 
m d = mass of C0 2 in sorbent material (kg) , 

mb = mass of sorbent in bed (kg) , 

T b = temperature of the bed (K) , 

T ref = reference temperature (K) , 

T g = temperature of C0 2 gas (K) (substitute T b ) , 

P g = pressure of C0 2 gas in bed void space (Kpa) , 

mg = mass of C0 2 gas in void space (kg) , 

R = C0 2 gas constant (kPa.m 3 / kg.K) , 

V g = void space of bed, also volume of C0 2 (m 3 ) , 

^ = transfer coefficient (0 < K 2 < 1) , 

m Q = C0 2 gas mass flow rate, determined by pump (kg/sec) , 

S c = heat of sorption of C0 2 (J/kg C0 2 ) , 

Power = power applied to bed (J/sec) , 
c vb = heat capacity of sorbent material (J/kg.K) , 
dm d /dt = rate of C0 2 desorbed (kg C0 2 /sec) , 
dnig/dt = change in mass of C0 2 in void space (kg C0 2 /sec) , 
dT b /dt = change in temperature of bed (K/sec) . 

Pump operation has been assumed to be 100% efficient, 
isentropic, and adiabatic. In actuality, no pump exhibits 
such behavior due to friction, windage, heat loss and 
pressure losses. For our first attempt, however, we will 
neglect these inefficiencies to simplify the model. At this 
time, this is an acceptable model as it provides the correct 
overall response. 

2) For the rotary vane pump assuming 100% efficiency, 
isentropic and adiabatic operation, we have 


mrm 0 , 

(7) 


(8) 



( 9 ) 
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( 10 ) 






The symbols are defined as: 

m 0 = mass flow rate of C0 2 gas stream (kg / sec) , 
w = angular velocity of pump (rad / sec) , 

k 3 = transfer coefficient of pump (0 < K 3 < l) , 

p g = density of C0 2 gas (kg / m 3 ) , 

R = C0 2 gas constant (kPa.m 3 / kg.K) , 

P g = entering C0 2 gas pressure to pump (kPa) , 

T b = entering C0 2 gas temperature to pump (K) , 

P p = exiting C0 2 gas pressure from pump (kPa) , 

T p = exiting C0 2 gas temperature from pump (K) , 

k = specific heat ratio of C0 2 . 

3) The accumulator, during continuous operation, is governed by 
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The needed quantities are defined by: 

T„ = temperature of accumulator (K) , 

U„ = internal energy of accumulator (J) , 
c v = specific heat of C0 2 gas (J / kg.K), 

T 0 = thermodynamic reference temperature for internal 
energy (K) , 

P p = C0 2 pressure at pump exit (kPa) , 

R = C0 2 gas constant (kPa.m 3 /kg.K) , 

V a = volume of accumulator (m 3 ) , 
m a = mass of C0 2 in accumulator (kg) , 

T p = exiting temperature of pump (K) , 
dnij/dt = change in total mass of C0 2 in accumulator (kg/ sec) , 
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dU a /dt = change in internal energy of accumulator (J/sec) , 
m 0 = C0 2 mass flow rate through pump (kg/ sec) , 

m^ = C0 2 mass flow rate to C0 2 Reduction Assembly (kg/ sec) , 
Qioss = heat loss of accumulator (J/sec). 

TRANSFER COEFFICIENTS AND INITIAL CONDITIONS 

The final values of the adsorbing cycle determine the initial 
conditions T b [0] and m d [0] for the desorbing cycle. From these 
values, the initial conditions are found. From Equation (1) one 
obtains 

fJO]"*,— (rjoi-r^ (15) 

m b 


Where the index of zero refers to the time relative to start of 
cycle. At t^ = 0 with P b [0] = P g [0], using equations (1), (3) 
and (15) yields the initial condition m^O] as 



(16) 


Since the initial amount of C0 2 in the desorbing bed is known and 
will be desorbed in about 60 minutes, the coefficient k 2 can be 
estimated by (5) to obtain 

^ s . m ^ . °l (1) = 0.0006 (17) 

60- min 


Next, to get k 3 , calculate p g [0] from (9) as 
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Taking the constant mass flow rate m 0 to be dm d /dt for w=200 
rad/sec. , k 3 can be found from (8) as 
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mM 

( 1 ' 

60- min 

Ipi" 


(19) 


*3*0.0023 


The k values are physical properties of the components and do not 
change. The values calculated above are not exact, but they are 
of the proper order of magnitude to allow model operation over a 
wide range of parameters. 
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MODELING TECHNIQUES 


Solution of the governing equations relies on solving 
differential equations. This simple integration procedure gave 
results accurate only for small time steps. The original 
simulation program used a six second time step. In regions where 
the solution has a steep slope, the large time step forced the 
program to over or under estimate the answer. The over and 
under-shoot appeared as transients on graphs of state variable 
equations. These transients also caused negative absolute 
pressures and negative absolute temperatures. The erroneous 
negative numbers were fatal errors in program statements with 
exponential functions. Decreasing the time step to one second 
for the AT&T 3B1 model and 0 . 1 second for the IBM models 
alleviated the problems caused by steep solution slopes. 

Decreasing the time step introduces new difficulties. Because 
data was saved in arrays, decreasing the time step increased the 
number of array elements. The limited memory of the AT&T 3Bl's 
prevented modeling of more than about ten minutes of the cycle. 

To overcome this problem, temporary registers were used to hold 
state equation variables and other variables until about thirty 
seconds of the model has been simulated. After the thirty 
seconds, the temporary register values are placed in plotting 
arrays. In the case of CLPMODEL.C, the temporary registers are 
used to collect data for five (length is arbitrary) minutes after 
which time it is written to CLIPS sensor files. 

Temperatures of the desorbing bed, C0 2 accumulator and C0 2 pump 
reached very high (>1000 K) levels in initial modeling attempts. 
Simple on/off control of the desorbing bed heater and cooling of 
the accumulator were used to keep the temperature within 
reasonable values. The controllers produce oscillations in bed 
temperature, accumulator temperature and associated components. 
Future versions of the program should incorporate control schemes 
reflecting actual controllers. 

SAMPLE OUTPUT 


Several simulations were made using the AT&T 3B1 version of the 
model, each time changing a single parameter. The model seems 
most sensitive to pump speed and reference temperature. The 
following set of graphs shows the model running with the listed 
parameters : 


Mass flow rate of air into model 
H 2 0 concentration into model 
C0 2 concentration into model 
Temperature of air into model 
Angular velocity of C0 2 pump 
Reference temperature of model 


0.2 kg/s 

0.01 kg H 2 0/kg air 
0.001 kg C0 2 /kg air 
300 K 
200 rad/s 
250 K 
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Figure 3.2-2 gives indication that a half-cycle time of 55 
minutes is sufficient for the desiccant bed to stabilize at or 
near the ambient temperature. The temperature of the adsorbing 
C0 2 bed follows the same curve as the desiccant bed. Because of 
simplifications made to the model, the maximum adsorbing bed 
temperature was nearly the same as the maximum desiccant bed 
temperature or approximately 322 K. 



Figure 3.2-3 Concentration of C0 2 Leaving 
Sorbent Bed 
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Figure 3.2-3 shows that most of the C0 2 is adsorbed in the first 
third of the cycle. The graph actually shows the return air 
concentration, but the amount of C0 2 adsorbed would be the 
initial concentration minus the graph value for a given 
simulation time. 



The C0 2 Removal model assumes that the C0 2 Reduction assembly is 
drawing C0 2 from the storage tank at a constant rate. Figure 
3.2-4 shows the accumulation of C0 2 in the tank. It was 
important to note that the Removal and Reduction assemblies could 
both place demands on the storage tank. As long as the Reduction 
assembly does not draw C0 2 from storage at a greater rate than 
the Removal assembly could produce C0 2 , the system will function 
properly. The control system must be able to resolve conflicts 
such as an empty C0 2 storage tank or a Reduction assembly that is 
drawing too much C0 2 from the tank. 
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The accumulator tank was modeled with a cooling jacket in place 
to remove a constant amount of heat when the tank reaches a set 
temperature. The jacket is controlled by a thermostat as the 
ECLSS control system should not bother with this low level 
control requirement. An interesting result of this on-off 
control of the tank temperature is the reflection of the graph 
shape in any other variable which has accumulation tank 
temperature as part of its controlling equation. The pump outlet 
pressure shown in figure 3.2-6 below illustrates this point. 



60 65 70 75 80 85 90 95 100 105 110 

Cycle Time (minutes) 


Figure 3.2-6 Outlet Pressure of C0 2 Pump 
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3.3 C0 2 REDUCTION ASSEMBLY MODEL 


DESCRIPTION 

The C0 2 Reduction Assembly of the ECLSS reduces the excess C0 2/ 
provides oxygen for use in the cabin and hydrogen for use in the 
Oxygen Generation Assembly. This is accomplished with a system 
comprised of a Sabatier Methanian Reactor (SMR) which converts 
C0 2 and H 2 to methane gas (CH 4 ) and water vapor (H 2 0) , a 
Condenser/ Separator (C/S) to remove product water, and a Carbon 
Formation Reactor (CFR) to break down methane into Carbon and 
Hydrogen. A block diagram of the C0 2 Reduction Assembly is shown 
in Figure 3.3-1. 



Figure 3.3-1 C0 2 Reduction Assembly 

HATH MODEL 

The following assumptions are made: 


1. Steady-state behavior. 

2. Plug flow behavior. 

3. No convection. 

4. No work is done by the system. 

5. Heat capacity is constant over the temperature range. 

6. No CH 4 or H 2 0 enters the gas stream. 

7. The reactions occur in stoichiometric amounts. 

8. Catalyst does not deactivate. 
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The C0 2 Reduction Assembly takes carbon dioxide and reacts it 
with hydrogen over a catalyst to produce methane and water. The 
reaction is given by the Sabatier reaction: 

C0 2 + 4H2 <==> CH 4 + 2H 2 0. (1) 

To find the volume of the catalyst (i.e. the reactor) necessary 
to carry out this reaction, a mass balance needs to be 
formulated. 


Mass Balance 



In words, the mass balance is as follows: 

Rate of accumulation rate at rate at rate of production 

of A in the = which A - which A + of A by chemical 
system enters leaves reaction 

In this derivation, A will refer to C0 2 . 

At this point, two assumptions are made. The system exhibits 
1) steady-state behavior and 2) plug flow behavior. Assuming 
plug flow behavior allows radial variations to be neglected. 

These requirements give 


0 = F a\ V ‘ ^l(v ♦ Av) + r A* V ’ 
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where 
F a ! v 


F a 


i 

I v+Dv 


r A 

dV 


= mqlar flow rate of C0 2 entering reactor 
(moles/sec) , 

= molar flow rate of C0 2 leaving the reactor 
(moles/ sec) , 

= reaction rate (moles/cm 3 ) , 

= change in volume of the catalyst (liters) . 


The molar flow rate of C0 2/ F A , can also be represented in terms 
of the conversion of C0 2 as 

F a = F^d - X A ), ( 4 ) 

where 

X A = conversion of C0 2 and F Ao = inlet flow rate of C0 2 . 

This gives 


dF A = d(F Ao - F Ao (X A ) ) = -F Aa (dX A ) . 


(5) 


Substituting this expression for dF A back into equation (4) 
yields 


dV = - 




( 6 ) 


Integrating both sides of equation (6) gives 



(7) 


Equation (7) is the mass balance integral that gives the volume 
of the catalyst necessary to carry out the reaction to be 
determined. In order to find the volume of the reactor, a 
reaction rate r A for this reaction must be determined. A 
reaction rate was found in a SAE technical paper [l:all]. The 
reaction rate used is given by 
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( 8 ) 


. 4 Jfc] 


pO-25 p 
r CO i r H 2 


pO-25 

r CH. 


p o.s 

r B,0 I 




where 

A = 1.112393 x 10 6 


— = 8 553°K 
R 


(9) 


K eq = exdO.50 


with 

T = Temperature 

P = Partial Pressure (atm) . 

In order to evaluate the integral in the mass balance (7) , r A 
must be exclusively in terms of X A . Thus, temperature and 
partial pressures must be in terms of X A . Temperature can be 
represented in terms of conversion through the evaluation of an 
energy balance. 

Energy Balance 

In words the energy balance is as follows: The rate of 
accumulation of energy equals the rate of flow of heat from 
surroundings, minus the rate of work done by system on 
surroundings, plus the rate of energy added by mass flow in, 
minus the rate of energy out by mass flow out. The standard 
energy balance for this system can be represented by the 
expression 


33 |5^0 + 34^3 _ 164ln(7 ) + 0.005577 


+ 33.165 


atm 


-2 


( 10 ) 


(f L - Q - W *- F ^ /J 9 ‘ C ‘ dT - 

(a»/ Pa * /J AC, dT j F m ■ X A , 


(ID 
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where 

Q = Heat flow from surroundings to system (cal/sec) , 
W, = Work done by system on surroundings (cal) , 

9; = Ratio of the number of moles of species i 

entering initially to the number of moles of 
A entering initially ( = FJF Ao ) , 

C pi = Heat capacity of species i (cal/gmoleK) . 

The following expression represents the heat of reation: 

T 

A H r ° [rj + f A C p dT = Heat of reaction ( cal/gmole ). 

t r 


This derivation requires the following approximations: there is 
no convection (Q = 0) , there is no work done by the system on the 
surroundings (W, = 0) , the heat capacity, Cpj, stays constant over 
the temperature range in question, and the system is in 
steady-state . 


The heat of reaction for this system was found in [1:1] to be 
given by 


A/7J7] = 16.4T - 0.00557-7* + - 1 -- — - - + 34633- cal 


gmole 


( 13 ) 


For this system the term. 


n r 

E / ^ C p,dT 

i=1 4 


( 14 ) 


expands to: 




I ^C pA dT 


F 

. T Ao 

\ L Ao 





( 15 ) 


where 
A = C0 2 

B = H 2 

C = CH 4 

D = H 2 0. 
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It is assumed that there is no CH 4 and H 2 0 in the entering stream, 
thus and are zero and the last two terms in equation (8) 
are also zero. The heat capacities for C0 2 and H 2 were repre- 
sented by linear functions with respect to temperature. The heat 
capacity for C0 2 is given by [3:1] 

C pA = a + bT, (16) 

where 

a = 5.316 

b = 0.014285 

for a temperature range of 273 K to 1800 K. 

The heat capacity for H 2 is represented by 

C pB = a + bT, (17) 

where 

a = 6.952 

b = -0.0004576 

for a temperature range of 273 K to 1800 K. 

It is also assumed that the reaction occurred in stoichiometric 
amounts, that is. 



(18) 


Substituting Equations (16) and (17) into (15) and integrating 
from T Ao = Tj* = 293 to T yields (19) : 

Fj(5.316T + 0.0071425-7* - 2170.8)+(27.808-r - 0.0009152-7* - 8069.2)]. H9) 

Substituting Equations (13) and (19) into the energy balance 
(14) , along with the assumptions made, yields: 

-FJ33.124-T + 0.0062273-7* - 10240] + 

16.4T - 0.00557-7* + 112 °°° + 34633 Jf^ =0. (20 > 

Solving for X A then gives 

x = 33.124-7’ + 0.0062273-7* - 10240 

16.4-T - 0.00557-7* + 112Q °° + 34633 < 21 > 

T 


After conducting a linear regression analysis, equation (21) 
reduces to: 

T = 1088X a + 297.59 (22) 
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Equation (22) now allows the temperature to be expressed in terms 
of X A in the mass balance. 

The next step is to relate the partial pressures of each species 
in terms of X A . This can easily be done by assuming ideal gas 
behavior and using the ideal gas law. The ideal gas law for 
component A (C0 2 ) is given by 


P A = C A RT, 

where 

P A = partial pressure of C0 2 (atm) 

C A = concentration of C0 2 (gmole/cm 3 ) 

R = gas constant 

T = temperature (°K) . 


The C0 2 concentration C A can be represented by 



(23) 


(24) 


where 

N a = number of moles of C0 2 

V = volume. 

The number of moles N A can then be represented by 

- M' - X A> 


(25) 


Since there is a change of moles in the Sabatier reaction, the 
volume can best be represented by 


V = 




(26) 


Epsilon, e, is defined as the fraction change in volume per mole 
of A reacted resulting from the change in the total number of 
moles. Epsilon can be expressed by 

n 

-E v, 

i=1 (27) 

e = , 

V* *To 
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where 

V; = stoichiometric coefficients which are negative 
for the reactant and positive for the products, 

N Ao = initial moles of CO 2 , 

N Xo = total moles initially present. 

Again the assumption that the reaction occurs in stoichiometric 
amounts is made so N Ao /N To is equal to 1/5 and e becomes -0.4. 
Substituting Equations (25) and (26) into (24) and using e = -0.4 
yields: 


"Jt - *a)T. 

V1 1 - 0.4X,)7- 


( 28 ) 


Substituting Equation (28) into (23) and noting that N Ao /V G is 
equal to C Ao yields 


Pa 



1 ) 

" 0.4X a> 



( 29 ) 


Noting that C Ao RT 0 is equal to P Ao and also noting that the 
temperature cancels out. Equation (29) becomes 




0AX a ‘ 


( 30 ) 


With the assumption of reactions occurring in stoichiometric 
amounts, P Ao is equal to 0.2. Thus, the partial pressure of A 
(C0 2 ) is now a function of conversion (X A ) . Following a similar 
derivation, the partial pressures for the remaining species can 
be represented as follows: 


Pb 




-4X a 

- oax A/ 


t 


Pc 




Pd 



(31) 


(32) 


(33) 
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Since all partial pressures and temperatures for all the species 
are in terms of conversion, the mass balance integral can now be 
integrated where r A is exclusively in terms of conversion. 


In equilibrium reactions, a maximum conversion can be achieved. 
This maximum conversion becomes the upper limit when evaluating 
the mass balance integral. The maximum conversion can be 
determined by setting the reaction rate, r A , equal to zero and 
solving for X A . Doing this yields a maximum conversion of 
approximately 0.52. This maximum conversion corresponds to a 
reactor operation temperature of approximately 860 K. Thus, the 
volume of the catalyst can be found by: 


V = 



(34) 


where 

r A = equation (8) 

temp. = given in terms of conversion by Equation (22) 
pp. = given in terms of conversion by Equations (30) 
through (33) . 


As can be seen by the complexity of the equations, this integral 
is difficult to evaluate analytically. Therefore a different 
method was employed. To evaluate this integral, a plot of l/-r A 
vs. X A was made and the area under the curve measured from X A = 0 
to 0.52. Using a software package, this area was found to be 
1.652 x 10 7 . Therefore the volume is 


V = 


1.652 x 10 7 

\ 


sec -cm 3 
gmoles 



(35) 


Equation (35) now allows us to determine the catalyst volume 
necessary to carry out the reaction for a specified inlet C0 2 
rate. Below is a table representing the reactor volume necessary 
for various man-loadings. 

Table 3.3-1 Reactor Volumes for Various Man-loadings 


Man Loading 

Inlet C0 2 Flow Rate 

Catalyst Volume 

3 - man 

7 x 10 4 mol/sec 

3.1 gallons 

5 - man 

1.3 x 10' 3 mol/sec 

5.7 gallons 

8 - man 

2.1 x 10' 3 mol/sec 

9.2 gallons 

10 - man 

2.8 x 10' 3 mol/sec 

11.3 gallons 


* Numbers obtained from SAE technical paper, # 840936. 
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Throughout this derivation there was one significant assumption 
that was made but not incorporated in the mass balance. This 
assumption was that the catalyst did not deactivate (i.e. the 
catalyst remained active indefinitely) . 

Based on this analysis, a reactor size of 11.3 gallons is 
recommended operating at a temperature of 860 K. A size of 11.3 
gallons was chosen because any man-loading below 10-man could 
easily be handled by this reactor. 

At this point in the modeling, a computer program was developed 
that would simulate the Sabatier reactor and it is listed in 
Appendix 3b. Output from this program correlated with the 
previously recommended reactor size. For a set reactor volume, 
any flow rate of C0 2 equal to or below the flow rate that 
corresponds to maximum conversion is fine. Yet, for flow rates 
that are just slightly larger, conversion drops exponentially. 

One item to note is that small changes in flow rates in moles/ sec 
correspond to fairly big changes in grams/ sec. Never-the-less, a 
reactor size corresponding to a 10-man loading will be sufficient 
for most worst-case scenario's. 


MODELING REVIEW 

- There are no controls in our model. 

- We did not consider transient responses in our model. 

This is because the Sabatier reactor has heaters to 
heat the catalyst at the entrance. We assumed that the 
reactor would be heated up to the temperature relating 
to maximum conversion before any species entered the 
reactor . 

- Obviously, our model is only a part of the overall WRM. 

We decided to assume that the C0 2 and H 2 coming into the 
Sabatier was actually from an external tank and not 
from another sub-assembly. Therefore, our model does 
not "communicate" with anything else in the WRM. 

- As far as the expert system is concerned, they can control 

bur whole model since the only output that is vital is 
conversion, which relates to the production of H 2 0 
forthe OGA. Depending on the needs for water for the 
OGA, or conversely the amount of C0 2 that the C0 2 
Removal sub-assembly is producing, the expert system 
can vary conversion to account for different flow rates 
of C0 2 in or of H 2 0 out. 


73 



3.4 OXYGEN GENERATION ASSEMBLY MODEL 


The Oxygen Generation Assembly (OGA) of the AR provides the daily 
requirements of oxygen to the cabin compartment of the space 
station. It accomplishes this by using electrolysis to convert 
water into its two components, oxygen and hydrogen. This water 
is supplied to the OGA from the C0 2 reduction sub-assembly and 
the hygiene water supply. A by-product of this reaction is 
hydrogen gas which is used in the C0 2 reduction sub-assembly of 
the ECLSS. 



The OGA assembly consists of two main components, a water storage 
and supply system and an electrolysis module. The storage and 
supply system of the OGA consists of two water storage tanks 
pressurized by nitrogen. These tanks supply water to the 
electrolysis module in alternating cycles. As one tank drains to 
the electrolysis module, the other tank fills with water. This 
is accomplished by means of a 4-way, 3-position directional 
control valve (DCV) . See Figure 3.4-1. The Differential 
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pressure across the bellows of the water tanks force the water to 
the electrolysis module membrane. When the draining tank reaches 
a low level set-point, detected by a photo-cell indicator, the 
DCV switches to drain the other tank and allows this tank to 
refill. This cycle continues throughout the normal operational 
mode of the oxygen generation process. 

When the water reaches the electrolysis module, it encounters a 
permeable membrane that separates the water supply from the water 
vapor region of the 18 electrolysis cells. Upon reaching this 
permeable membrane of the electrolysis module, the water diffuses 
across it and forms water vapor on the other side. The 
electrolysis module then produces H 2 from the water vapor at the 
cathode and passes OH" ion through a KOH matrix. 0 2 is then 
formed on the other side of this matrix at the anode of the 
module. Both gases then leave the electrolysis unit and pass 
through the pressure control assembly. See Figure 3.4-2. This 
pressure control assembly, consisting of pressure control valves 
on both the 0 2 and Hj outgoing lines, regulates the differential 
pressure of the two departing gases. 



Figure 3.4-2 Electrolysis Module 
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The reation of interest at the cathode is 

4ff 2 0+4e'=2H 2 +40ir f 

which for the anode one obtains 

40ir=0 2 +2tf 2 0+4e'. 


HATH MODELING 

The modeling of the OGA has been separated into two components: 
the water storage tanks and the electrolysis module. Each of 
these have their separate governing equations which are covered 
in their respective areas. 

Water Tank 


The mathematical modeling for the water tanks can be done by 
modeling the water portion, and the nitrogen portion of each 
tank. 

Assumptions 

1. Water has constant physical properties. 

2. Nitrogen has a constant heat capacity. 

3. The water flow into the storage tank can be considered 
turbulent since the tank fills in a four to five 
minute period through a small diameter pipe. 

4. The water flow out of the storage tank can be considered 
laminar since it flows out over a three hour period 
through a similar sized pipe. 

5. The nitrogen flow into and out of the storage tank can 
be considered as turbulent during the fill mode; and is 
laminar during the drain mode. 

6. Nitrogen follows the ideal gas law. 

7. All heat transfer by conduction and radiation will be 
considered negligible because of minimal temperature 
differences in the system (i.e. no heat loss) . 

Math Modeling of Water Portion of Tank 

When modeling the water portion of the water tank, the dependent 
variables of interest are the change in the volume and 
temperature, and the inlet and outlet water flows. 
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A material balance can be written to relate the change in volume 
of the tank with respect to time. This material balance can be 
written as: 

Rate of Mass Accumulation = [Rate of Mass - [Rate of Mass] 0Dt . 

Symbolically this becomes 


P«j 

at 


where 

V X1 = volume of water contained in the storage tank at 
any given time (m 3 ) , 

w H2o,i = mass flow rate of water in (Kg/ sec) , 
w h 2 o, 2 = mass flow rate of water out (Kg/sec) , 

Pmo = density of water (Kg/m 3 ) . 


During the tank fill mode, the Vl mo2 term in (1) is equal to zero. 
In the drain mode, the W ffi0 i term is equal to zero. 


Since the tank takes four to five minutes to fill, it is a good 
assumption that the water flowing into the tank is turbulent 
because the Reynolds number (Re) is approximately 40,000. In 
deriving an equation for the turbulent flow of water through a 
tube, a good starting point is the Blassius Formula for friction 
factors. The Blassius Formula is as follows: 


f _ 0.0791 
J?e 1/4 


( 2 ) 


for, Reynolds numbers in the range 2.1(10) 3 < Re < 10 5 , 
where 

Re = the Reynolds number given by Re= D <v> Ph 2 o/Mh 2 o» 
D = the tube diameter (m) , 

<v> = the average velocity of the water (m/sec) , 

Mh 2 o = the viscosity of water (Kg/m sec) . 


When the definition for the Reynolds number is substituted into 
the Blassius Formula, Equation (2) becomes 


f= 


0.0791 


D<V> 


Ph 2 o 

Ph 3 o 


(3) 


The friction factor can also be written as 


£_ D(P gouzce ~Pvi) 


2Lp„o<v>‘ 


( 4 ) 
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where 

P source = the pressure of the external water source (N/m2) , 
P V1 = the pressure in the water side of the tank (N/m2) , 
L = the tube length (m) . 


Setting Equations (3) and (4) equal to each other gives 

0.0791 


f=- 


D ( Pgour ce~Pvi^ 


D<V> 


2Lp* j0 <v> 5 

Ph 2 o 


( 5 ) 


Multiplying each side by the tube cross sectional area(A=7rD 2 /4) 
and then rearranging gives 


P h 2 o< v >A= 


(0 .1582Lp„ j0 1/4 ) 4/7 


( 6 ) 


Substituting Equation (6) into Equation (1) gives 

P h 2 c @ V T1 _ (13 5 ^ 4 p ffj0 (P S ource - Pvi^ )*^ 7 A 
dt (0.1582 L[1 Hj0 1/4 ) 4/7 


(7) 


since 


^ h 3 o~ P A . 


With the equation derived for the change in volume of the water 
tank during the fill mode with respect to time, an equation that 
models volume change when the water is flowing out of the tank is 
needed. The water in the tank flows out over a period of three 
hours, so laminar flow is assumed. The laminar flow of a fluid 
through a tube is given by the Hagen-Poiseuille law as 


q . _ 71 (Pyi~£ > elec') r * 


( 8 ) 


where 

Peiec = the pressure in the electrolysis module. 

Altering Equation (8) can be changed into mass flow rate by 
multiplying by p mo to obtain 

rr 71 (P VI “Peiec) ■ Z ’ 4 P»j<? 

2 flu.. JL 


Substituting Eq. (9) into Eq. (7) , gives the change in volume of 
the water portion of the water tank during the drain mode as 
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( 10 ) 


Ph 2 C@Vti _ 71 elect Z i Pn 2 0 

at 8 1 x Hi0 l 

With these equations that model the change in volume of the tank 
during the fill and drain modes, an equation needs to be derived 
to relate the change in temperature of the water with respect to 
time. This can be done by writing an energy balance on the 
system. The energy balance for the water portion of the water 
tank is given by: 

[Rate of Energy] = [Rate of Energy ];„ - [Rate of Energy] out 

-[Net Rate of Heat Loss to Surroundings] . 

Symbolically, this becomes: 

d(V T1 (T H20 -T zef )) _ 

dt ~ (ID 

W H 2 0, 1 C P„ 2 o ( ^H 2 0 , 1 - Tzef ) - 2 C P rf2 o ^ ^O, 2 ~ T ref> ' 


where 

c P H 2 o = heat capacity of Water (KJ/Kg K) , 

t h 2 o,i = temperature of water into the compartment (K) , 

Th 2 o ]2 = temperature of water out of the compartment (K) . 


Note that the net rate of heat loss to the surroundings equals 
zero. Assumptions were made that the system was well insulated 
and any heat loss to the surroundings is negligible. 


Equations for the mass flow rate of water in and out of the tank 
have already been derived. By manipulating the left hand side of 
Eq. (11) a relationship can be drawn between the rate of change 
in temperature and volume given by 


, , dV T1 

C Ph 2 <Ph 2 o ( t h 2 o T zeft -~rr- + C } 


Ph 2 o 


?h 2 o v ti~ 


dT 

dt 


( 12 ) 


When solved for <ST/£t, Eq. (12) relates the change in temperature 
of the water in the tank with time. This term should be 
negligible if the temperature of the water and the temperature of 
the surroundings are approximately the same. Assuming this, the 
only variables of interest in the system become mass flow rate of 
the water and the rate of change of volume in the water tank. 

Nitrogen Portion of the Water Tank 

When modeling the Nitrogen portion of the water tank, the 
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variables of interest are as follows: 


w N2,i = ®ass flow rate of Nitrogen in (Kg/sec) , 

Wn 2^2 = mass flow rate of Nitrogen out (Kg/ sec) , 

Tj^i = temperature of Nitrogen in (K) , 
t n^ 2 = temperature of Nitrogen out (K) , 

p n2!i = pressure of the Nitrogen in the Water tank (N/m 2 ) , 
P N2^ = pressure of the Nitrogen in the supply tank (N/m 2 ) , 
V^’ = volume of the Nitrogen side of the tank (m 3 ) . 


These variables can be related by differentiating the Ideal Gas 
Law with respect to time to give 


( 13 ) 


where 

dV-re/dt equals -dV T1 /dt since 

= V T1 + Vtc 

or: 

^ tank _ dV T1 &V T 2 

dt dt dt ' 


( 14 ) 


( 15 ) 


The Sn m / S t term can be determined by making a mole balance of the 
Nitrogen compartment given by 

[Accumulation of Moles] = [Rate of Moles in] -[Rate of Moles out] 

Symbolically this can be written as 

= <1/M, U«> 


where 

n^j = moles of Nitrogen, 

M = molecular weight of Nitrogen. 


When the water side of the tank is filling, the Nitrogen side of 
the tank is draining. This process occurs over a four to five 
minute period. Assuming turbulent flow for the nitrogen in this 
case is a safe assumption. By analogy, Eq. ( 6 ) applies to the 
flow rate of nitrogen leaving the tank given by 




V-P P/fr ^N2,X~^H2.2> > 
(0 . 1582L[i w 1/4 ) 4/7 


'A 


( 17 ) 


where 
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p m = density of the nitrogen at the given pressure (Kg/m 3 ) , 
= viscosity of the nitrogen (Kg/m sec) . 


The flow of nitrogen into the tank can be considered laminar due 
to the extended period of time of the nitrogen fill. Similar to 
the case above, a simple manipulation of the equivalent water 
side produces the equation for nitrogen flow into the tank of 


W, 


N t . 1 


_ 1 1 ( P N 2 2 ~ P N 3 l) r *PN 2 
8\l„L ~ 


( 18 ) 


Knowing these two flow rates, and assuming that the heat loss to 
the surroundings is negligible, an equation similar to Eq. (11) 
can be derived from the energy balance of the nitrogen portion of 
the tank to give 


PV 


S<V„{T„ r T r . t )) 

at 




where 

CpN 2 = Heat capacity of Nitrogen (KJ/Kg K) 

Expanding the left-hand side of Eq. (19) gives 

dV_, dT„ 

T *t ~ T zef ) -q£- + c Plf 2 v raPif 2 -Jjr 


( 20 ) 


However , knowing that 


. dVr2 _ dn * 
P ^~dT~~dF' 


( 21 ) 


the left side of Eq. (19) reduces to 


dT, 


n 2 


W N 2 .l~ W N 2 ,2 + C Pla V T2Ptr 2 dt 

Solving for dT m /dt give 

^ t n 2 _ W N2il ( c p ia (T Ni ' 1 -T zgf ) - 1 ) +W N3 ' 2 (l-C pia (T Ni 2 -T zaf ) ) 

C Pla V T 2 Plf 2 


dt 


( 22 ) 


(23) 


With equations for dV^/dt, dn m /dt , and dT m /dt, an equation for 
dP m / dt can be derived by substitution into Eq. (13) . This 
equation will yield the temperature, pressure, moles of nitrogen 
and the volume of the tank at any time. The temperature term is 
constant because of the negligible difference in the nitrogen and 
ambient temperatures. 
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Electrolysis Module 


The assumptions for this model are as follows: 


1. Water has constant physical properties. 

2. Water is an incompressible fluid. 

3. The heat capacities of the gases are constant. 

4. There is no accumulation of gases in the electrolysis 
module . 

5. "The non-electrolytic 'parasitic' energy used by the 
electrolyses unit per kilogram of water converted to 
gas" is negligible. 

6. The amount of water vapor leaving in the product gases 
is negligible. 

Initially, the electrolysis module will be modeled as a single 
cell. Modeling could be accomplished by dividing the module into 
a feedwater compartment and the actual electrolysis compartment. 
Variables of interest in modeling the water compartment are: 


Wfco.i = mass flow rate of water into module (Kg/ sec) , 
w h 2 o !2 - mass flow rate of water out (Kg/ sec) , 

Pmo - density of water (Kg/m3), 

V E1 = volume of water contained in the water compartment 
of the electrolysis module (m 3 ) . 


A material balance on the feedwater module produces the following 
equation: 


PhjC^^ei 

dt~ 




(24) 


The volume of the water compartment in the electrolysis module 
remains constant at steady state, therefore the above equation 
can be reduced to 


^H20,l — ^H20,2 * (25) 

The flow out of the water tank is considered laminar, so it 
follows that the flow into the water compartment of the 
electrolysis module is laminar as well. This allows Eq. (9) to 
be used for the flow rate of water into the feedwater compartment 
which takes the form 
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*( P tank~ P elec 


>* 4 P*, 


H a o, 1 




At steady state, the vapor rate (W m0t2 ) is equal to the rate of 
water out of the feedwater module. The flow rate out could also 
be estimated by deriving an expression for the diffusion of a 
vapor through a porous membrane. However, very little 
information exists about the membrane and the above equations 
provide a better representation. 

Writing an energy balance for the water compartment produces and 
equation similar to those used for the water tank. However, the 
equation for the water compartment contains a term for the 
vaporization of water. The energy balance is written as follows: 

dT no 

Ph 2 0 C p HjO V E1 = W H 2 0, 1 C p !h0 ( T H 2 O t 1 " T zef> ~ "up, 2 C p H2 o ( T H 2 0. 2 ~ T zef> 

~@cond. El ~ Qrad. El + ( ~^) ^ H 2 0 . 2' ( 26 ) 


where 

t h 2 o,i = temperature of water into the compartment (K) , 
Ticcu = temperature of water out of the compartment (K) , 
Q CO nd,Ei = heat loss by conduct ion (KJ/ sec. ) , 

Qrad,m - heat loss by radiation (KJ/sec. ) , 

-h = heat of vaporization of water (KJ/Kg) . 


The equations for flow rate in and out of the water compartment 
can be replaced by Eq. (27) . The heat loss by conduction through 
the tank walls can be expressed as 


Qcond. El A El 


So^H 2 C 


lnr B ifr Eo 

Krr 


■t-Bl^gurr 




( 27 ) 


where 

r^ = inside radius of the tank (m) , 
r E1 = outside radius of the tank (m) , 

A m = surface area of the water compartment, excluding side 
with membrane (m 2 ) . 

Heat loss by radiation from the tank can be written as follows: 

^rad.El^^Bl^El 0 ^ ‘ P H 2 0~' P suzr^ • ( 28 ) 


83 


where 

e E1 = emissivity of the outer compartment walls. 

The temperature in the feed water compartment is known at any 
particular time. Also, the flow rate into and out of this 
compartment can be calculated. 

When looking at the electrolysis compartment of the module, there 
are three variables that need to be considered, i.e., the 
production of hydrogen, oxygen, and the change in temperature in 
the module. It becomes necessary at this point to work on a mole 
basis. Writing a mole balance for hydrogen and oxygen yields the 
form 


[moles accumulated ]= [moles in] - [moles out] + [moles produced]. 


By dividing Eq. (26) by the molecular weight of water, the number 
of moles entering the electrolysis compartment are obtained. 

This number of moles equals the number of moles of H 2 and half 
the number of moles of 0 2 leaving the system. The equations for 
these two variables can be stated as 


^h 2 o. i 


^ ^ tank ^ elec ^ ^ Ph 2 0 

8(18 )\i H j0 L 


( 29 ) 


where 

qico.i ~ The molar flow-rate of water into the electrolysis 
module (moles/sec) . 

Knowing that the moles of hydrogen and oxygen produced are 
dependent upon the electrolysis rate given by 


q=I/F*z , 


( 30 ) 


where 

q = moles of gas produced (mols/sec) , 

I = current supplied to the electrodes (col/sec) , 
F = Faraday's constant (9.64 8 7 x 10 4 ) , 
z = equivalents /mole (H2=l, 02=2 ) , 

the moles of hydrogen produced can be written as 

I(col/sec) 

Hl ' 2 (9 . 6487 xl0 4 col/ eqv. ) (2eqv.H 2 /mole) 


( 31 ) 


Here hydrogen has 2 equivalents per mole. The stoichiometric 
relationship of the dissociation of water into hydrogen and 
oxygen is 


h 2 o => h 2 + %0 2 . 


( 32 ) 
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Therefore the number of moles of oxygen produced equals one-half 
that of the amount of hydrogen produced. One can say that the 
rate of production of moles of oxygen equals one-half that of 
hydrogen which is written as 

qo 2 ,out = %qm,ouf (33) 

With this relationship, the flow rate of water into the 
electrolysis module for a given current can be calculated, 
assuming the accumulation of moles in the module equals zero. By 
knowing the flow rate of the water, a pressure drop from the 
water tank to the electrolysis can be also found. 


An energy balance will produce the temperature changes in the 
electrolysis module. This energy balance is similar to the 
energy balances calculated earlier. However, it also accounts 
for the energy supplied to the system through the electrodes 
which operate at approximately 137 °F. The energy balance can be 
expressed as 


P avcr C Pam V Bl' 


dT, 


*** = ^h 2 o, i c p„ 2 o ( t h 2 o, i ~ T ref ) 2 (2) 


at 


C PH 2 ( ' Th 2 Tref ' > 


-^2 02) c^T^-T^) -Q, 


cond, El 


-Q 


rad, Bl + E elec^O^O, 




(34) 


where : 

T H20,1 

T H2,2 

^02,2 

Pavg 

C pavg 

C p02 

C pH2 

Qcood,E2 

Qrad,E2 

®“elec 


= temperature of water into the compartment (K) , 

= temperature of hydrogen out of the module (K) , 

= temperature of oxygen out of the module (K) , 

= average density of the three gases (Kg/m 3 ) , 

= average heat capacity of the gas (KJ/Kg K) , 

= average heat capacity of the oxygen (KJ/Kg K) , 

= average heat capacity of the ^ (KJ/Kg K) , 

= heat loss by conduction (KJ/sec) , 

= heat loss by radiation (KJ/sec) , 

= electrical energy added to the system by the 
electrodes (KJ/Kg) , 

= the chemical energy added to the system due to 
the heat of the reaction. 


The heat loss through conduction can be expressed as follows: 


? Cond,B2 & B2 ' 


r E2^H, 


^E2 

k2 


r B3^-guTz 


] _1 (^-r surr ) 


(35) 


where 

r^ = inside radius of the electrolysis module (m) , 
r^ = outside radius of the electrolysis module (m) , 

Aj ?2 = surface are of electrolysis module, excluding side 
with the membrane (m 2 ) . 
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( 36 ) 


Heat loss by radiation from the tank can be written as 

Qrad, E2~^B2^B2° ( “ -^surr ) • 


where 

6 122 = emissivity of the electrolysis module. 

The energy added to the system by the electrodes can be 
represented by the following equation: 

((297 3.-?) +E p ) M w 

jg = £ 

Blac ( M w +M v ) 2 . 7 87 8 xlCT 4 ' 


(37) 


where 

Ep = "the non-electrolytic 'parasitic' energy used by 

the electrolyses unit per kilogram of water converted 
to gas, 

V = cell Voltage, 

I = electrolyses unit current efficiency, 

M„ = mass of water electrolyzed to gas (Kg) , 

^ = mass of water leaving the system (Kg) . 

The value of Ep is less than 0.05% of the total energy, and is 
therefore negligible. The My term is a function of temperature 
and pressure and is essentially constant at the temperature and 
pressure of interest. This allows a safe assumption that a 
negligible amount of water vapor leaves the system with the 
product gases. The previous equation can then be simplified to 

E^ = 10702 V/I. (38) 

The chemical energy removed from the system due to the heat of 
reaction can be written as 

pzoductB ~M£ (H°-Ho+AHf 0 ) reactants' (39) 

where 

H° - H° 0 = the enthalpy of the element at 298K and the system 
temperature, respectively (KJ/mole) , 

AH f ° 0 = the heat of formation of the compound from its 
elements . 

Oxygen Generation Assembly 

The rate of production of 0 2 in the OGA is directly related to 
the electrolysis rate, which in turn is directly dependent upon 
the amount of current supplied to the electrolysis module. This 
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assumption is the basis for the modeling of the entire OGA. 

The rate of electrolysis depends upon the following relationship: 



( 40 ) 


where 

q = rate of gas produced (moles/sec) , 

I = current supply to the electrodes (col/sec) , 

F = Faraday's Constant (9.6487 x 10 4 (col/eqv) ) , 
z = equivalents /mole (H 2 = 1, 0 2 = 2) . 


For example, the moles of hydrogen produced can be written as 


2, out - 


I(col/sec) 

(9 .6487 xio* (col/ eqv) ) ( 

mole 


( 41 ) 


According to the stoichiometric relationship between H 2 and 0 2 is 


H 2 0=>H 2 +±0 2 . 


( 42 ) 


Therefore the moles of 0 2 produced equals half the moles of H 2 
produced . 

MODELING TECHNIQUES 

The above assumptions that the rate of 0 2 and H 2 productions do 
not take into account any temperature or pressure effects on the 
electrolysis rates. After studying the complexities of the 
pressure and temperature effects and seeing how they related to 
the production results, it was determined that for simplicity, 
they would be ignored. In essence, this assumption creates a 
simple black box with known inputs and outputs for the 
electrolysis module. It also provides a simple linear 
relationship between the current and the production rate of 0 2 . 
The model eliminates any time based differential equations that 
would otherwise govern the output of the system. 

Based on the above assumptions and knowing the requirements for 
0 2 production, sizing of the water storage and supply side of the 
OGA can be done. During normal (4 -man) operational mode, the OGA 
needs to produce 9.081b of 0 2 /day. Knowing this 0 2 is generated 
from 10.22 lb of H 2 0/day, a rough tank size was determined to 
hold 1 lb of water. During the emergency (8 -man) operational 
mode, the system needs to generate 15.6 lb of 0 2 /day [3:1-11]. 
Since the water tank fills in only 3 min. and drains in over 70 
min. , this tank size meets water capacity needs and space 
requirements . 
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4.0 EXPERT SYSTEM 
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4.1 INTRODUCTION 




Monitoring and maintaining the Environmental Control and Life 
Support System (ECLSS) on the Space Station Freedom (SSF) 
requires the effort and expertise of many supervisors, engineers, 
and technicians. If the station personnel is to provide this 
expertise, then they will be less able to provide the expertise 
and labor needed to perform their mission. However, if ECLSS can 
intelligently-monitor itself and perform self-diagnosis, then the 
crew's limited resources can be better applied to the mission. 
This year, the design team investigated the application of expert 
systems to provide for intelligent automated control and 
monitoring of ECLSS. The prototype expert system the teams 
developed is discussed in the following sections. 

Besides the ultimate effect of providing automated life support 
expertise for the SSF crew, an expert system controller provides 
special benefits to the system designer. The use of an expert 
system shell allows developers to program the system's knowledge 
base using a high-level specification, object oriented, or rule 
based style language. The format of the knowledge base can be 
customized to suit the knowledge particular to a given 
application. Furthermore, expert systems can be developed with 
incomplete and inexact information. This feature permits expert 
systems to function in poorly characterized applications such as 
non-linear modeling and natural language interpretation. Because 
an expert system can operate on incomplete information, it is 
easier to produce and improve working prototypes. Improving the 
prototype is a matter of adding additional and more specific 
knowledge . 

The final goal of the expert system design teams is to provide an 
intellect controller for the whole of ECLSS. This year, the team 
developed a prototype expert system for controlling the 
Atmosphere Revitalization called ARES. As an expert system, ARES 
has two major components, a knowledge base and an inference 
engine. The three design teams developed knowledge bases for 
controlling and monitoring the C0 2 Reduction, C0 2 Removal, and 
Oxygen Generation subassemblies. Additionally, one of these 
teams developed the inference engine to drive the expert system. 

A more detailed description of ARES follows. 
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4.2 ATMOSPHERE REVITALIZATION EXPERT SYSTEM (ARES) 


The function of ARES is to intelligently monitor and control the 
Atmosphere Revitalization (AR) assembly. Because this is a 
hardware control application, ARES must receive sensor signals 
from the assembly, process these input signals, and provide 
control signals to the assembly. Because ARES is an expert 
system, the sensor signals are processed after the manner of an 
expert. That is, as an expert, ARES translates the sensor values 
into terms compatible with its knowledge base, applies its 
knowledge to develop responses and diagnostics as an expert 
would, and implements its responses by providing the proper 
control signals to the assembly. An explanation facility enables 
ARES to justify its conclusions to the crew. 

As an expert system, ARES has two parts, an inference engine 
shell and a knowledge base. The inference engine of an expert 
system is a shell that separates the general knowledge and 
support software from the application specific knowledge. It is 
the inference engine that does the '•thinking" for the expert 
system. The knowledge base of an expert system contains the 
knowledge that represents the expertise for the given 
application. 

ARES' s inference engine supports and interprets a near-english, 
rule-based knowledge base. ARES "thinks” by deductively chaining 
conditional statements and thereby connects general sensor values 
to specific conclusions. ARES's general line of reasoning is as 
follows: "What sensors should I read? What numerical values do 
these sensors currently offer? What linguistic values may I use 
to describe these values, how are these linguistic values 
defined, and how accurately can they be applied? What 
conclusions can I draw from the linguistic descriptions I just 
developed? Based on these conclusions, what instructions should 
I send to the assembly. How should I resolve possible conflicts 
in the instructions?" These questions are generally answered by 
facts and rules found in the knowledge base. However, numeric 
sensor values are read from data files. The general "meta- 
knowledge" that is needed to interpret and chain rules, evaluate 
the "fuzzy" logic embedded in the rules, and resolve conflicts 
between instructions is encode into the inference engine shell. 

ARES's knowledge base contains the knowledge needed to control 
the Atmosphere Revitalization assembly. Rather than expressing 
the control function as sets of feedback equations as is done for 
conventional controls, the control function is express as a set 
of conditional rules. Composed in a near-english format these 
rules are written for ARES as one would write instructions to a 
human performing the same job. The general format for a 
conditional rule is 
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(cclassif ier>: if <statement> then <conclusion> <weight>) , 
where <statement> and <conclusion> are generally of the form 

<noun> is/is_not <modifier>. 

Note <statement> and <conclusion> can really be of any form but 
it is preferred that they convey meaning to both developers and 
users. 

The inference engine interprets conditional rules by searching 
through the knowledge base for a fact that matches the 
<statement>. The <statement> has a weight value (0-1) associated 
with it called the confidence value. The confidence value 
represents the confidence an expert would have in the <statement> 
for the given circumstances. ARES is constructed so that the 
<statement> can always be found (eventually) although its 
confidence value may be zero. Having found the <statement>, the 
inference engine asserts the <conclusion> into the knowledge base 
and gives it a confidence value. The confidence value for the 
conclusion is calculated by multiplying the confidence value 
found for the <statement> by the <weight> value associated with 
the conditional rule. A future version will support multiple 
layers of conclusions and drop the structural restrictions 
imposed by the <classif ier>: convention. 

Once asserted by the inference engine, the <conclusion> can be a 
valid <statement> for another rule. In this manner, conditional 
rules can be chained. But how are initial statements provide to 
start the chains? By membership rules that translate numeric 
sensor values into descriptive linguistic values. The general 
form of a membership rule is 

(if <sensor> is X then <noun> is <modifier> 

<min> <low> <high> <max> <units>) , 

where <sensor> is the name of a sensor listed in the knowledge 
base. The "<noun> and <modifier>" is the conclusion that the 
inference engine asserts. The numeric values <min> <low> <high> 
and <max> define a triangular "fuzzy" membership function used to 
calculate the confidence value for membership rule's conclusion. 
If the value read from the <sensor> is between <low> and <high>, 
then the confidence in the conclusion is absolute (its confidence 
value is 1) . If the value read is less than <min> or greater 
than <max>, then the confidence in the conclusion is null (its 
confidence value is 0) . All other possible values for the 
<sensor> are linearly interpolated to produce a continuous 
function. A <units> is just a term used for the reference 
developer and user. 

ARES supports statements of the form 
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of the 



<statement A> <operator> <statement B>, 

where <operator> is either OR/ AND thereby permitting fuzzy logic 
The confidence of the statement is resolved according to fuzzy 
algebra. 

ARES also support statements of the form 
<noun> is_not <modifier>. 

The confidence value is the fuzzy logic's complement of the 
statement "<noun> is <modifier>. " 

It is the nature of an expert's decision process to decide 
between alternative or conflicting options. If the options 
represent points values along a roughly continuous scale, the 
conflict can be resolved by taking the centroid of the points as 
a compromise value (use confidence values as weights) . But, in 
ARES, the subassemblies only accept a few discrete non-scaler 

values: ON, OFF, STANDBY, ROLLOVER, PLAYDEAD, ARES resolve 

conflicts by choosing the most confident of the conflicting 
values. 

IMPLEMENTATION : 

In APPENDIX 4D is the listing for SAVANT. 3. SAVANT. 3 is a 
program written for NASA's CLIPS expert system development shell 
The first section of SAVANT. 3 's code is the particular knowledge 
base used to verify the system. All the following sections of 
SAVANT. 3 comprise the inference engine that would be transparent 
to the user. 

Expanded knowledge bases are presented in APPENDIXES 4A-4C. 

These knowledge bases were developed by three teams for the C0 2 
Removal Assembly, C0 2 Reduction Assembly, and Oxygen Generation 
Assembly. The assembly knowledge bases together contain rules 
used to identify current operation conditions and diagnose 
problems in the AR. They are based on expertise of how the AR 
sub-assemblies should be controlled and how ARES should respond 
to expected conditions. 

The following sections discussed the rules of each assembly 
knowledge base in depth. 
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4.2.1 C0 2 REMOVAL ASSEMBLY 


The C0 2 Removal Assembly is responsible for removing excessive 
carbon dioxide from the cabin air. This removal is accomplished 
by a four-bed molecular sieve. The sieve contains two beds for 
removing moisture from the incoming air and two beds for removing 
the carbon dioxide. In addition to the four beds, the system 
contains a blower for moving the air, a precooler, a C0 2 pump and 
accumulator, and a variety of flow control valves. A system 
diagram of the assembly is shown in Fig. 4.2-1. 



Figure 4.2-1 C0 2 Removal Assembly 


The purpose of this section is to indicate the work that has been 
done on designing an expert system to monitor and control the 
assembly. The rules were written to be compatible with the CLIPS 
expert system package and are located in Appendix 4B. 

The main function of the designed expert system was to monitor 
the components of the C0 2 Assembly. Different sensor values are 
read and compared to certain parameters and the expert system 
makes a decision based on those sensor values. The focus of the 
rules was to detect component failures. Currently, the rules do 
not include any valves or the four beds. 

Rule # 1 

The expert system monitors the temperature sensor at the output 
of the precooler. If the value does not fall within accepted 
parameters, the expert system makes a decision if failure has 


93 




occurred or not. If a failure decision is reached, the expert 
system instructs the Removal Assembly to shut down. 

Rule # 2 

The expert system monitors the two pressure sensors located on 
either side of the blower and checks for the differential 
pressure. If the pressure is unacceptable, the Removal Assembly 
is instructed to shut down. 

Rule # 3 

The expert system monitors the power sensor on the C0 2 pump. If 
the value is unacceptable, then the Removal Assembly is 
instructed to shut down. 

Rule # 4 

The expert system monitors the pressure sensor located on the C0 2 
pump. If the value is not within the accepted parameters, the 
Removal Assembly is instructed to shut down. 

Rule # 5 

The expert system monitors the gas flow at the C0 2 accumulator 
through two gas flow sensors located on either side of the 
accumulator. If a failure decision is reached, both the Removal 
and the C0 2 Reduction subassemblies are instructed to shut down. 
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4.2.2 C0 2 REDUCTION ASSEMBLY 


The C0 2 Reduction Assembly is comprised of a Sabatier Methanian 
Reactor (SMR) which converts C0 2 and H 2 to methane gas (CH 4 ) and 
water vapor (H 2 0) , a Condenser/ Separator (C/S) to remove product 
water, and a Carbon Formation Reactor (CFR) to break down methane 
into carbon and hydrogen. A block diagram of the C0 2 Reduction 
Assembly is shown in Figure 4.2-2. 



Figure 4.2-2 co 2 Reduction Assembly 


The main focus of this section will be on the steps that were 
taken to construct an expert system for the SMR. The rules that 
were written in CLIPS will be explained in the order that the 
reactants flow through the SMR. See Appendix 4C for the rules as 
found in CLIPS. 

The main purpose of the SMR is to catalytically react C0 2 with H 2 
to form methane gas (CH 4 ) and water vapor (H 2 0) . The SMR receives 
C0 2 from the C0 2 Removal Assembly and the H 2 is recycled from the 
carbon formation reactor and is also received from the C0 2 
Removal and Oxygen Generation Subassemblies. The chemical 
reaction that takes place in the SMR is written as 

C0 2 + 4H 2 = CH 4 + 2H 2 0 + HEAT. 

A 4.0 molar ratio, which is equivalent to a molecular weight of 
10.4 g/gmol of H 2 to C0 2 , must be maintained. To do this a 
molecular weight sensor was already included in the prototype 
SMR, so this was not included in the expert system sensors. The 
actual volume of C0 2 was needed to accurately control the 
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molecular weight sensor which was not readily available and 
difficult to approximate. 

Rule # 9 

The chemical reaction to convert C0 2 and H 2 to methane and water 
is exothermic. Therefore the heaters used in the SMR are 
primarily for reactor start up. Two redundant heaters are used 
to heat up the catalyst bed to initiate the reaction [Forsythe, 
1984]. A sensor to indicate a non-functional heater or heaters 
was included in the expert system. When the reaction is below 
400°K which would normally be during reactor start-up, there is 
no conversion taking place. This would indicate to the system 
that a heater or heaters has broken down and the system is turned 
off so maintenance procedures can take place. 

Rules # 7 and 8 

The chemical reaction is also self-limiting via the reverse 
endothermic reaction, to temperatures below 873°K [8:]. This 
causes a loss in the normally expected 98% conversion of carbon 
dioxide. Thus the efficiency of the reactor decreases. To 
monitor this, two rules were added to the expert system. The 
efficiency level is low below the temperature of 850°K. When 
this occurs the expert system instructs the SMR to shut off since 
low efficiency could be an indication of another problem. A 
normal efficiency level was also included in the rules when the 
temperature of the reaction was between 850°K and 873°K. The SMR 
would continue to run as normal if this rule was executed. 

Rules if 1, 2, 3, and 4 

The next step in the SMR is to cool the chemical reaction 
temperature down to the point that water can condense after it 
has left the reactor. To do this a cooling jacket consisting of 
a number of stainless steel serrated fins covered with a shell of 
stainless steel is wrapped around the reactor bed and cooling air 
is forced through a mid and rear cooling jacket inlets 
[Forsythe, 1984] . During the simulations of zero gravity on the 
SMR it was found that a coolant flow of approximately 0.424 cfm 
would be required for each jacket to compensate for the loss of 
external convection [1:]. Since the only temperature modeled by 
the C0 2 Reduction Assembly modelers was the reaction temperature, 
it was difficult to assume a temperature for the two cooling 
jacket zones. The coolant flow did not seem to be an adjustable 
variable. It was either on or off. Using these assumptions the 
rules were based on having a coolant flow of 0.424 cfm. If 
either of the jacket flows was below 0.2 cfm the heat would stay 
in the reaction and the exit temperature would be too high for 
the condenser/ separator . This indicated a problem with one or 
both of the coolant flows. Therefore the system would be shut 
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down for maintenance. There is also a normal flow mode between 
0.4 and 0.7 cfm when the SMR continues running. 

Rules # 5 and 6 

Although the condenser/ separator has not been modeled, an attempt 
to include it in the expert system was made. Since the condenser 
requires the water to be in a liquefied state it was assumed that 
the exit temperature of the SMR would have to be approximately 
373°K. Therefore, if the exit temperature was above this value 
then the condenser temperature would be too high and the SMR was 
turned off. A high exit temperature could be a result of the 
coolant flows not working. If the exit temperature was below 
373°K then the SMR was considered running normally. 

The Carbon Formation Reactor was not included in the C0 2 
Reduction Assembly expert system. There was limited information 
on the CFR and there had also been no attempts to model it. 


97 



4.2.3 OXYGEN GENERATION ASSEMBLY 


The Oxygen Generation Assembly, OGA, is comprised of a 
electrolysis module and water tanks to convert H 2 0 to H 2 and 0 2 by 
the use of a catalysis, KOH. A block diagram of the OGA is shown 
in figure 4.2-3. 


Electrolysis Module 

Permeable 



To 02 
Supply 


To N2 
Supply 


Figure 4.2-3 Oxygen Generation Assembly 


The main focus of this section will be on the steps that were 
taken to construct an expert system for the OGA. The rules that 
were written in CLIPS will be explained in the order that the 
reactants flow through the OGA. See Appendix 4 A for the rules as 
found in CLIPS. 

The main purpose of the OGA is to produce 0 2 by the electrolysis 
module. The OGA receives water from the C0 2 Reduction Assembly 
and hygiene water. For the purpose of writing an expert system 
to control the Atmosphere Revitalization Assembly, It was assumed 
that all water supply would come from the C0 2 Reduction Assembly. 
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Rules #1/ 2, and 3 


The water in the storage tanks must be moved to the electrolysis 
module by the use of nitrogen pressure on the bellow of the water 
tanks. This pressure is the H 2 0 feed line pressure. From Boeing 
document 683-10011, the H 2 O feed line pressure needed was given 
as 25 ±5 psia. Rule 2 defines the normal pressure to the 
electrolysis module between 20 and 30 psia. Rules 1 and 3 are 
for high and low feed line pressure. If these rules are 
activated, the OGA will be instructed to shutdown for 
maintenance. 

Rules # 4 and 5 

To measure the amount of water in the H 2 0 supply tanks, an 
emitter detector infrared sensor was used. If the bellow was 
obstructing the infrared sensor a binary output of one was the 
output of the sensor. If the water supply in the tank was low, 
the bellow would not obstruct the infrared sensor and the sensor 
outputs a binary output of 0. In this manner the supply of water 
for the electrolysis module was measured. An output of one 
indicates water supply is sufficient. Where an output of zero 
indicates the water supply is insufficient. 

Rules #6, 7, and 8 

After the H 2 0 passes through the electrolysis unit its components 
are H 2 and o 2 . The differential pressure between the H 2 and 0 2 
output lines is an indication of the electrolysis efficiency. 
Boeing document 683-10011 specifies the H2~0 2 differential 
pressure to be 2.83-3.35 psia for normal operation. If the 
differential pressure is lower or higher than this specified 
value it is an indication that one of the cell electrodes has 
failed and the OGA is instructed to shutdown. 

Rules #9, 10, and 11 

The final constraint placed on the OGA is the temperature of 0 2 
returned to the cabin. Although the temperature of 0 2 gas is not 
an output from the models designed by our group, it has been 
included in the expert system. The acceptable output temperature 
listed in Boeing document 683-10011 is 55°-105.8° F. If the 
temperature of 0 2 gas to the cabin is outside this range the OGA 
is instructed to shutdown. 

Rules # 12 and 13 

The electrolysis unit uses a series of voltage cells to control 
and initiate the breakdown of H 2 0. If all cell voltages are high 
meaning greater than 1.7 volts, as listed in NASA document 2- 
H8RG-BMS-194/90, this is an indication of that the water supply 
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to the electrolysis module has been lost. Therefore, if all cell 
voltages are high the OGA is instructed to shutdown. 

An important section of ARES is the Explanation Facility. This 
feature is not necessary for operation of ARES, but very 
important. ARES is only as accurate as the experts of each 
assembly. In other words, it is no smarter than its programmers 
and therefore is suspect to incorrect conclusions. The 
Explanation Facility allow the user to see the steps taken by 
ARES in making a conclusion for each assembly. It does this by 
retracing the paths used to chose responses. Then ARES generates 
a dialogue explaining the responses, allowing the user to refute 
the decision. 
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4.3 VERIFICATION AND TESTING OF EXPERT SYSTEM 


ARES is designed to help monitor and control the Environmental 
Control and Life Support System for the Space Station Freedom 
(SSF) . Currently, the system only has control over limited 
portions of the Air Revitalization System. After rules that 
govern system's behavior were written and placed into the 
knowledge base, the expert system was tested to discover any 
programming, knowledge, or conceptual errors that may have been 
present . 

The final stage of the expert system's development is 
verification - applying a series of inputs and insuring the 
decisions made were correct. Though the current system contains 
only four components that may be varied, verification is a 
difficult task. To verify all possible condition combinations 
would require around 300 individual tests. The complete expert 
system to govern the entire ECLSS will be nearly impossible to 
verify completely. For this report, 12 different sets of 
conditions were tested and checked. These twelve conditions sets 
fall into nine cases as described below. For each case, the 
conditions are stated first, then a brief description of the 
system response is given. The four conditions which can be 
varied are: 

Cabin_0 2 - Indicates the concentration of 0 2 in the cabin. 

Helps control Oxygen Generation Assembly (OGA) 
and C0 2 Removal Assembly 

Inlet_C0 2 - Indicates concentration of C0 2 present in the 
cabin atmosphere. Helps control C0 2 Removal 
Assembly. 

C0 2 Accumulator - Indicates level of C0 2 storage tank 

located between C0 2 removal and C0 2 
reduction assemblies. Helps control 
these same two assemblies. 

H20_Accumulator - Indicates level of water in storage 

tank located between C0 2 reduction and 
OGA assemblies. Helps control these 
same two assemblies. 


CASE 1: 

Condition: All values in normal operating range 
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Response: 

OGA - commanded to turn on (if currently off) or to continue 
in normal operating mode. 

C0 2 Reduction - commanded to turn on (if currently off) or to 
continue in normal operating mode. 

C0 2 Removal - Commanded to move to an efficiency mode. Since 
no supplies are low and cabin C0 2 is normal, C0 2 
Removal Assembly can reduce power requirements by 
decreasing production slightly. 


CASE 2: 

Condition: Both accumulators are empty, all other values in 

normal operating range 


Response: 

OGA - commanded to turn on (if currently off) or to continue 
in normal operating mode. This system will not, 
however, be able to operate if no water is available. 
This system should perhaps turn off for a while until 
water is available. This is the correct response as 
programmed although not really the correct way to 
handle the given situation. 

C0 2 Reduction - commanded to turn on (if currently off) or to 
continue in normal operating mode. This is to refill 
the empty H20 accumulator. 

C0 2 Removal - Commanded to move to high production mode. 

This will help to replenish the C0 2 accumulator and 
provide resources needed by the C0 2 Reduction assembly 
to refill H20 Accumulator. 

CASE 3: 

Condition: Both Accumulators at a low level, all other values 

in normal operating range 


Response: 

OGA - commanded to turn on (if currently off) or to continue 
in normal operating mode. This system will not, 
however, be able to operate if no water is available. 
This system should perhaps turn off for a while until 
water is available. This is the correct response as 
programmed although not really the correct way to 
handle the given situation. 
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C0 2 Reduction - commanded to turn on (if currently off) or to 
continue in normal operating mode. This is to refill 
the empty H20 accumulator. 

C0 2 Removal - Commanded to move to high production mode. 

This will help to replenish the C0 2 accumulator and 
provide resources needed by the C0 2 Reduction assembly 
to refill H20 Accumulator. 

CASE 4: 

Condition: Both Accumulators at a high level, all other values 

in normal operating range 


Response: 

OGA - commanded to turn on (if currently off) or to continue 
in normal operating mode. 

C0 2 Reduction - commanded to turn on (if currently off) or to 
continue in normal operating mode. 

C0 2 Removal - Commanded to move to an efficiency mode. Since 
no supplies are low and cabin C0 2 is normal, C0 2 
Removal Assembly can reduce power requirements by 
decreasing production slightly. 


CASE 5: 

Condition: Inlet_C0 2 at a high level, all other values in 

normal operating range 


Response: 

OGA - commanded to turn on (if currently off) or to continue 
in normal operating mode. 

C0 2 Reduction - commanded to turn on (if currently off) or to 
continue in normal operating mode. 

C0 2 Removal - Commanded to move to a high removal mode. 

This is an attempt to remove the excess C0 2 from the 
cabin and restore a proper C0 2 - 02 balance for the 
astronauts . 


CASE 6: 

Condition: Inlet C0 2 at a low level, all other values in 

normal operating range 


103 



Response: 


OGA - commanded to turn on (if currently off) or to continue 
in normal operating mode. 

C0 2 Reduction - commanded to turn on (if currently off) or to 
continue in normal operating mode. 

C0 2 Removal - Commanded to turn off. This is done to avoid 
overfilling the C0 2 accumulator before the need to 
remove C0 2 becomes important. This also conserves 
power at times when the system is not needed. 

CASE 7: 

Condition: Cabin_02 level is low, all other values in normal 

operating range 

Response: 

OGA - commanded to turn on (if currently off) or to continue 
in normal operating mode. Since there is no high 
production mode, the expert system can only make sure 
the OGA is producing oxygen as the environment 
requires . 

C0 2 Reduction - commanded to turn on (if currently off) or to 
continue in normal operating mode. 

C0 2 Removal - Commanded to move to an efficiency mode. Since 
no supplies are low and Cabin C0 2 is normal, C0 2 
Removal Assembly can reduce power requirements by 
decreasing production slightly. This will not affect 
the rate of oxygen production which is most important 
in this test case. 

CASE 8: 

Condition: Cabin 02 at a critical level, all other values in 

normal operating range 


Response: 

OGA - commanded to turn on (if currently off) or to continue 
in normal operating mode. Again, we would desire a 
high production mode, but none is available. We can 
only make sure the OGA is on and operating. 

C0 2 Reduction - commanded to turn on (if currently off) or to 
continue in normal operating mode. 
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C0 2 Removal - Commanded to move to a high removal mode. This 
is done to help restore C0 2 - 02 balance in cabin and 
to insure that adequate resources are available for the 
OGA to continue operating normally. 


CASE 9: 

Condition: All values at a high level in the operating 

range 


Response : 

OGA - commanded to turn on (if currently off) or to continue 
in normal operating mode. 

C0 2 Reduction - commanded to turn on (if currently off) or to 
continue in normal operating mode. 

C0 2 Removal - Commanded to move to a high removal mode. This 
is done to compensate for excess C0 2 in the cabin. 

SPECIAL NOTE: For this section, one of the input values, the 

C0 2 accumulator, did not get read into the 
system. The results may have been different if 
this value was available. 

Limitations 

The biggest limitation at present is the number of rules which 
the system can operate on. With only four pieces of data to make 
decisions with, those decisions will be limited in scope and 
accuracy. Additionally, some things have not been taken into 
account. For example, the system will continue to remove C0 2 
from the cabin even if the C0 2 accumulator is full. Next, the 
system produces confidence factors although these factors 
currently have no meaning. Finally, the expert system does not 
always produce the best decision. As noted in cases 2 and 3 
above, the OGA is commanded to turn on even though adequate 
resources for operation may not be available. It is evident that 
further options need to be considered and tested before a final 
decision is made as to the operation of the system. 
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4.4 CONCLUSIONS 


Because of the complex nature of the conditions which may affect 
operation of the ECLSS, standard control schemes do not provide 
convenient or adequate means for system development and 
maintenance. The expert system, on the other hand, can handle 
systems with multiple inputs and outputs and non-linear behavior. 
Expert systems can be tolerant of inexact or incomplete 
information permitting rapid prototyping. With an expert system 
in charge ECLSS will not require constant supervision (a premium 
commodity on a space flight) . For these reasons, an expert 
system solution makes a better choice for ECLSS control. 

The expert system constructed this semester used NASA's CLIPS 
expert system development shell. The CLIPS shell program does 
not readily support real-time modeling or color graphics user 
interfacing. Future effort should include a search for 
development tools that can provide this support. 

The rules defined in the previous sections obviously do not 
represent a set which will completely control the Air 
Revitalization assembly. Rather, these rules are a beginning set 
which define the overall or general operation of the assembly's 
main components. Many other, more detailed rules will be needed 
to complete the knowledge bases for the expert system and to 
provide an adequate control system for the ECLSS. Further 
information about the operation and interaction of the various 
sub-assemblies will be required before such rules can be defined. 

The Savant. 3 decision mechanism is the most complete portion of 
the current expert system design. The inference engine, which 
interprets knowledge base data, and the response selector, which 
chooses an action based on the number of conclusions and their 
associated confidence factors, will not be affected by the 
addition of new rules. In fact, these sections could be used 
with minor alteration in an expert system for control of the 
entire ECLSS. 

System testing validated the programming of the SAVANT. 3 
inference engine and its general application of deductive 
reasoning. This testing also validated the technique of 
representing the knowledge base. While the condensed version of 
the control scheme generally worked, the expert system faithfully 
demonstrates certain fallacies in the scheme. But this real 
objective of the ECLSS simulation — to help verify the 
conceptual control schemes. 
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APPENDIX 2 A 

C0 2 REMOVAL ASSEMBLY PROGRAM 


* * 


★ 

★ 

* 

★ 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

★ 


SOURCE FILE: REMOVE 2 • C AT&T 3B1 Version 


* 

it 


DESCRIPTION: 


USES: 


This program models the C02 removal system. * 
It simulates one adsorption and desorbtion * 
cycle. The desorbing desiccant bed has not * 
been included in this model. This program * 
can be used as a design tool to test the * 
C02 removal system equations. * 

* 

main() * 

simple_plot ( ) (S. Dyer) * 

★ 


AUTHORS: Robert Young/ Terry Hon, Joel Sommer 

LAST REVISION: 06 March 1991 


* 

* 

it 

it 


* LIMITATIONS : 

* 

* 

it 

it 

it 

* 

* 

★ 

* 

* 


The desorbing desiccant bed has not been * 
included in this model because the CLIPS * 
Expert System is not concerned with the * 
rehumidification of cabin air. Only the * 
levels of C02 into the Removal System, the * 
C02 released to the cabin, the C02 level of * 
the accumulator and the C02 demand of the * 
Reduction System are currently necessary for * 
CLIPS to make an inference about mode * 
transitions. * 

★ 


*********************************** **************************** ^ 


#include <stdio.h> 
/include <math.h> 
/include "p_plot.h" 
/include "simple_plot . h” 
/define NUM POINTS 100 


main ( ) 

{ 

int dummy, pen; 
long int i ; 
char rite; 
double Time [100], 
xTime, 

Mg [ 100 ] , 
xMg, 

Md[ 100 ] , 
xMd, 

Tb[100] , 
xTb, 

Ua[ 100] , 
xUa, 

Ma[100] , 
xMa, 

Pp[100], 

xPp, 

Tp(100] , 


/** time array **/ 

/** mass of gas in void space **/ 

/** mass of gas in sorbent material **/ 

/** temperature of bed and exiting air **/ 
/** internal energy of accumulator **/ 

/** mass of C02 in accumulator **/ 

/** output pressure of pump **/ 

/** output temperature of pump **/ 
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xTp, 

Ta[ 100 ] , 
xTa, 

Conc_w[ 100 ] , 

xConc_w, 

Ccs_w, 

M_f low, 

Ti, 

T_des [ 100 ] , 
xT_des, 
xT_fan, 
xT_hx, 

Cci_w, 

Cco_w, 

Md_w[ 100 ] , 
xMd_w, 

Ccs, 

Cci, 

Cco, 

Cone [100] , 
xConc , 

T_sorb[100] , 
xT sorb ; 


/** temperature of accumulator **/ 

/** relative humidity of outgoing air **/ 

/** cone, of H20 on sorbent (kgH20/kgair ) **/ 
/** mass flow rate of air (kg/sec) **/ 

/** incoming air temperature (K) **/ 

/** temp, of air leaving desiccant bed (K) **/ 

/** temperature of air after fan (K) **/ 

/** air temp, after heat exchanger (K) **/ 

/** cone, of H20 entering bed (kgH20/kgair ) **/ 
/** cone, of H20 leaving bed (kgH20/kgair ) **/ 
/** mass of H20 on sorbent (kg) **/ 

/** cone. C02 in sorbent (kgC02/kg air) **/ 

/** cone. C02 entering bed (kg C02/kg air) **/ 
/** cone. C02 leaving bed (kg C02/ kg air) **/ 
/** cone. C02 leaving bed (ppm) **/ 

/** temp, of air leaving sorbent bed (K) **/ 


double Pg, /** C02 pressure in void space **/ 

Pb, /** pump exit pressure **/ 

m_o, /** mass flow rate out of bed dMg/dt ** / 

m_out, /** assembly exiting flow rate (bang-bang) **/ 

w, /** angular velocity of pump (rad/sec) **/ 

Q_loss, /** accumulator heat loss (bang-bang) **/ 

Tref, /** reference temp, of desorbent bed **/ 

k_w, /** ratio of Cco w/Ccs_w **/ 

k_s; /** ratio of Cco7ccs **/ 

double Power=1000 . 0, 

Mb=20.0, /** mass of sorbent material (kg) **/ 

Cpa=1003.0, /** heat capacity of air (J/kgK) **/ 

Cvs=844.0, /** heat capacity of sorbent (J/kgK) **/ 

Sc_w=2700.0E3, /** heat of sorption of H20 

into sorbent (J/kgK) **/ 

Cvb=657.0, 

Cp = 846.0, 

R=0 . 1889 , 

To=0.0, /** reference temperature for internal energy **/ 

Vg=1.0, 

Va=2 . 0 , 
kl=0.7, 
k2= . 0006 , 
k3=.0023, 

k4= 2.0, /** accumulator heat transfer coef. (Watts/K) **/ 

k = 1.29, 

Sc=951E3; /** C02 heat of sorption into sorbent (J/kgK) **/ 

double dt, 
dTb, 
dT_des , 
dT_sorb , 
dMd_w, 
dMg, 
dMd, 
dMa, 
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dUa; 


/** Inputting of parameters by user **/ 

printf { "\nlnput mass flow rate of air, M_flow (kg/s) : "); 
scanf (”%lf",£M_flow) ; 

printf ( "\nlnput H20 concentration entering assembly,"); 
printf (" Cci_w (kgH20/kg air) : "); 
scanf ( "%lf " , &Cci_w) ; 

printf ( "\nlnput concentration of C02 entering assembly,"); 
printf (" Cci (kgC02/kg air): "); 
scanf ( "%lf " , &Cci) ; 

printf ( "\nlnput temperature of incoming air, Ti (K) : "); 
scanf ( "%lf " , &Ti) ; 

printf ("\nlnput angular velocity of pump (radians/second) : "); 
scanf ( "%lf " , &w) ; 

printf ("\nlnput desired bed reference temperature, Tref (K) : "); 
scanf ( ”%lf " , &Tref ) ; 

Time(0] = 0.0; /** absolute start time in minutes **/ 

xTime = 0.0; 

/** initial conditions for desiccant bed **/ 

Md_w [ 0 ] = xMd_w =0; 

T_des[0] = xT_des - Ti; 

Ccs_w = Md_w [ 0 ] / Mb; 
k_w = 0.1; 

Cco_w = k_w * Ccs_w; 

Conc_w(0] = xConc_w = Cco_w * 10000.0; 

/** initial conditions for sorbent bed **/ 

Md(0] = xMd = 0.0, 

T_sorb [ 0 ] = xT_sorb = Ti; 

CCS = Md(0] / Mb; 
k_s =0.1; 

Cco = k_s * Ccs; 

Conc[0] = xConc = Cco * 1.0E6; 
dt = 1.0; 

printf ( "\n\nWorking desiccant and sorbent bed forward cycle. .. \n\n” ) ; 

^************************************************* 
it * 

* Desiccant bed and Sorbent bed calculations * 

* * 

************************************************* y 

i = 1; 

while (i<100) { 

for ( dummy= 1 ; dummy<=3 3 ; dummy++ ) < 

/** desiccant bed equations **/ 

dT_des = dt * ( (M_flow * Cpa * (Ti - xT_des)) + 

(M_flow * Sc_w * (Cci_w - Cco_w) ) ) / (Mb * Cvs); 
dMd_w = M_flow * (Cci_w - Cco_w) * dt; 
xT_des = xT_des + dT_des; 
xMd_w = xMd_w + dMd_w; 

Ccs_w = xMd_w / Mb; 

Cco_w = k_w * Ccs_w; 
xConc_w = Cco_w * 10000.0; 

/** fan and heat exchanger "equations" **/ 
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/** this is a kludge **/ 


xT_fan = xT_des + 6.81; 
xT_hx = xT_des - 50.0; 

/** sorbent bed calculations **/ 

dT sorb = ( (M_flow * Cpa * (xT_des - xT_sorb)) + 

(M_f low * Sc * (Cci - Cco) ) ) * dt / (Mb * Cvs); 
dMd = M_flow * (Cci - Cco) * dt; 
xT_sorb = xT_sorb + dT_sorb; 
xMd = xMd + dMd; 

Ccs = xMd / Mb; 

Cco = k_s * Ccs; 
xConc = Cco * 1.0E6; 

xTime = xTime + dt/60.0; 

> 

Time[i] = xTime; 

T_des(i] = xT_des; 

Md_w [ i ] = xMd_w ; 

Conc_w(i] = xConc_w; 

T_sorb(ij = xT_sorb; 

Md[i] = xMd; 

Cone [ i ] = xConc ; 
i = i + 1; 


^***************************************************** 

* * 

* Plotting routine for Desiccant and Sorbent beds * 

* Dr 

***************************************************** J 

pr intf ( " \n\nDone ..."); 

print f ( "\n\nPlease choose which device you wish to use"); 
printf(" for plotting purposes. \n" ) ; 
print f ("\nDISPLAY = 1"); 
print f ( "\nPLOTTER =2"); 
printf ( "\n\n" ) ; 
scanf ("%d",&pen) ; 
if (pen == 1) 
rite=DISPLAY; 
else rite=PLOTTER; 

/** Simple plot calls for desiccant bed data **/ 

simple_plot (rite, i,Time,T_des, "Time" , "minutes" , "Temperature" , 
"degrees Kelvin” , 

"Temperature of Desiccant Bed (T_des[i] ) ",1, CURVE) ; 
scanf ("%d",&dummy) ; 

simple_plot (rite, i, Time, Md_w, "Time", "minutes" , "Mass H20" , "kg" , 

"Mass H20 in Sorbent Bed Material (Md_w[ i] )", 1, CURVE) 

scanf ( "%d" , &dummy) ; 

simple_plot (rite, i,Time,Conc_w, "Time" , "minutes", 

"Relative Humidity" , "%" , 

"Relative Humidity of Air Leaving Desiccant Bed”, 

1, CURVE) ; 

scanf ("%d",&dummy) ; 
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/** Simple plot calls for sorbent bed data **/ 

simple_plot (rite, i,Time,T_sorb, "Time" , "minutes" , "Temperature" , 

"degrees Kelvin" , "Temp, of Air Leaving Sorbent Bed", 
1 , CURVE ) ; 

scanf ( "%d" , &dummy) ; 

simple_plot (rite, i,Time,Md, "Time" , "Minutes" , "Mass C02" , "kg" , 

"Mass C02 in Sorbent Bed Material (Md[i] ) " ,1, CURVE ) ; 

scanf ( "%d" , fidummy ) ; 

simple_plot (rite, i, Time, Cone, "Time" , "minutes" , "Concentration" , 
"ppm", "Cone. C02 Leaving Sorbent Bed" , 1, CURVE) ; 


scanf ( " %d" , fidummy ) ; 


system ( "clear" ) ; /** clear the screen **/ 

printf ( " \n\nWorking on desorbing cycle. . . " ) ; 

/** initial conditions for desorbing cycle **/ 


Time ( 0 ] 
Md[0] = 
Tb[0] = 
Mg[0] = 
xMg 

Ua[0] = 
Ma[0] = 
Ta[0] = 
Pp[0] = 
Tp[0] = 


= xTime = 62; /** reset array for this cycle **/ 
xMd = Md[99]; /** reset array for this cycle **/ 
xTb = T_sorb(99]; /** reset array for this cycle **/ 
kl * Vg * Md[ 0 ] * <Tb[0] - Tref) / (Mb * R * Tb[0]); 
: Mg ( 0 ] ; 


xUa = 407.34; 
xMa = 0.002; 
xTa = 300.0; 
xPp = 0.0567; 
xTp = 300.0; 


/********************************************** 
★ * 

* Calculations for Desorbing cycle * 

* * 


i = l; 

while ( i<100) { 

for ( dummy =1 ; dummy<=28; dummy++ ) { 


if (xTb > 330.0) /** bang-bang bed temp, controller **/ 

Power=0 . 0 ; 

else Power = 1000.0; 


Pg = xMg * R * xTb / Vg; 

Pb = kl * xMd /Mb * (xTb - Tref); 

xPp = xMa * R * xTa / Va; 

xTp = xTb*pow( (xPp/Pg) , (1.0-1. 0/k) ) ; 

xTa = xUa / (xMa * Cp) + To; 

m_o = w*k3*Pg/R/xTb; 

if (xMa < Ma[0]) /** bang-bang flow out control **/ 

m_out = 0.0; 

else 

m_out = 0.75 * m_o; 

if (xTa > 600.0) /** bang-bang accumulator cooler **/ 

Q_loss=k4* (xTa - 300); 
else Q loss=0.0; 
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dMd = (Pg - Pb) * k2 * dt; 

dMg = ((Pb - Pg) * k2 - (w * Pg * k3)/R/xTb) * dt; 
dTb = (dMd * Sc / dt + Power ) /Mb/Cvb * dt; 
dMa = dt * (m_o - m_out); 
dUa = dt * (m_o* (Tp[ i-l]-To) *Cp-m_out* 
(Ta(i-l]-To)*Cp - Q_loss) ; 

xMg = xMg + dMg; 
xMd = xMd + dMd; 
xTb = xTb + dTb; 
xMa = xMa + dMa; 
xUa = xUa + dUa; 
xTime = dt/60.0 + xTime; 

} 


Time [ i ] = xTime ; 
Mg[ i] = xMg; 

Md( i] = xMd; 

Tb[ i] = xTb; 

Ma ( i ] = xMa ; 
Ua[ij = xUa; 

Ta[ i] = xTa; 

Tp[ i] = xTp; 

Pp( i] = xPp; 


i = i + 1; 

> 


printf ( "\nPress number and return to begin plotting..."); 
scanf ( ”%d” , Sdurnmy ) ; 

/** SIMPLE PLOT CALLS FOR DESORBENT CYCLE DATA **/ 

simple_plot ( rite , i , Time , Md , "Time", "minutes”, 

"Mass C02", "kg", "Mass C02 in Sorbent Bed Material (Md[i])", 
,1, CURVE); 

scanf ( "%d" , fidummy ) ; 

simple_plot (rite, i, Time, Ma, "Time", "minutes", 

"Mass C02", "kg", "Mass C02 in Accumulator (Ma[i]) n 
,1, CURVE) ; 

scanf ( " %d" , &dummy ) ; 

simple_plot (rite, i,Time,Ua, "Time" , "minutes" , 

"Internal Energy", "J", 

"Internal Energy in Accumulator (Ua[i])" 

, 1 , CURVE ) ; 

scanf (”%d",&dummy) ; 

simple_plot (rite, i,Time,Tp, "Time", "minutes" , "Temperature" , 

"K", "Temperature of Pump Outlet (Tp[i])", 1, CURVE); 

scanf ( "%d" ,&dummy) ; 

simple_plot (rite, i, Time, Ta, "Time" , "minutes" , "Temperature" , 

"K" , "Temperature of Accumulator (Ta[i] ) ",1, CURVE ) ; 
scanf ( "%d" , fidummy) ; 

simple_plot (rite, i,Time,Pp, "Time" , "minutes” , "Pressure" , 

"kPa" , "Outlet Pressure of Pump (Pp( i] )", 1, CURVE) ; 
scanf ( "%d" , dummy) ; 
printf ( "\nDone sending plots..."); 
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/*************************************************************** 


* PROGRAM: CLPMODEL . H VERSION: 1.0 * 

★ * 

* DESCRIPTION: Global variables for clpmodei.c * 

* Hr 

* USES: none * 

* * 

* PROGRAMMER: Robert W. Young * 

* * 

* LAST REVISION: 03 March 1991 * 

* * 




#define NORMAL 0 /** function main should return a value to DOS **/ 

/** GLOBAL VARIABLES **/ 
double time, tstop; 
int pair; 
int dummy; 
long int i; 


double xTime, /** time in cycle **/ 

xMg, /** mass of gas in void space **/ 

xMd, /** mass of gas in sorbent material **/ 

xTb, /** temperature of bed and gas desorbed **/ 

xUa, /** internal energy of accumulator **/ 

xMa, /** mass of C02 in accumulator **/ 

xPp, /** output pressure of pump **/ 

xTp, /** output temperature of pump **/ 

xTa, /** temperature of accumulator **/ 

xT_fan, /** temperature of air after fan (K) **/ 

xT_hx, /** air temp, after heat exchanger (K) **/ 

xConc_w, /** relative humidity of outgoing air **/ 

Ccs_w, /** cone, of H20 on sorbent (kgH20/kgair) **/ 

M flow, /** mass flow rate of air (kg/sec) **/ 

Ti, /** incoming air temperature (K) **/ 

xT_des, /** temp, of air leaving desiccant bed (K) **/ 

Cci_w, /** cone, of H20 entering bed (kgH20/kgair) **/ 

Cco_w, /** cone, of H20 leaving bed (kgH20/kgair ) **/ 

xMd_w, /** mass of H20 on sorbent (kg) **/ 

Ccs, /** cone. C02 in sorbent (kgC02/kg air) **/ 

Cci, /** cone. C02 entering bed (kg C02/kg air) **/ 

Cco, /** cone. C02 leaving bed (kg C02/ kg air) **/ 

xConc, /** cone. C02 leaving bed (ppm) **/ 

xT_sorb; /** temp, or air leaving sorbent bed (K) **/ 

double Pg, /** C02 pressure in void space **/ 

Pb, /** pump exit pressure **/ 

m_o, /** mass flow rate out of bed dMg/dt **/ 

m_out, /** assembly exiting flow rate (bang-bang) **/ 

w, /** angular velocity of pump (rad/sec) **/ 

Q_loss, /** accumulator heat loss (bang-bang) **/ 

Tref, /** reference temp, of desorbent bed **/ 

k_w, /** ratio of Cco w/Ccs_w **/ 

k_s; /** ratio of Cco7ccs **/ 

double Power, 

Mb=20.0, /** mass of sorbent material (kg) **/ 

Cpa=1003.0, /** heat capacity of air (J/kgK) **/ 

Cvs=844.0, /** heat capacity of sorbent (J/kgK) **/ 

Sc_w=2700.0E3, /** heat of sorption of H20 

into sorbent (J/kgK) **/ 
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Cvb=657 .0, 
Cp = 846.0, 
R=0 . 1889 , 
To=0.0, 
Vg=1.0, 
Va=2.0, 
kl=0.7, 
k2=.0006, 
k3=.0023, 
k4= 2.0, 
k = 1.29, 
Sc=951E3 ; 


/** reference temperature for internal energy **/ 

/** accumulator heat transfer coef. (Watts/K) **/ 
/** C02 heat of sorption into sorbent (J/kgK) **/ 


double dt, 
dTb, 
dT_des , 
dT_sorb , 
dMd_w, 
dMg, 
dMd, 
dMa, 
dUa; 

FILE *clips, 

♦model ; 
char *mode ; 


/** file pointer for storage of sensor data **/ 
/** file pointer for storage of model data **/ 
/** string to hold CLIPS mode instructions **/ 


^*************************************************************** 


PROGRAM: 


CLPMODEL.C 


DESCRIPTION: Simulate C02 removal sub-assembly using math 

modeling equations derived in Fall 1990 report 
of "Automation of Closed Environments in Space 
for Human Comfort and Safety". 

Program will provide five minutes of sensor 
data to an ASCII file to be read by CLIPS and 
acted upon. CLIPS will generate an ASCII com- 
mand file which will direct operation of the 
model . 


USES: 


* 

* 

* 
it 
it 
it 
* 

★ 

★ 

* 

* 

* 

* 

* 
it 
* 

* PROGRAMMER: Robert W. Young 

* 

* LAST REVISION: 04 March 1991 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 


VERSION: 1.1 


stdio.h 
math . h 
"clpmodel.h” 


LIMITATIONS: 


The desorbing desiccant bed has not been 
included in this model because the CLIPS 
Expert System is not concerned with the 
rehumidification of cabin air. Only the 
levels of C02 into the Removal System, the 
C02 released to the cabin, the C02 level of 
the accumulator and the C02 demand of the 
Reduction System are currently necessary for 
CLIPS to make an inference about mode 
transitions. 


* 

★ 

* 

* 

★ 

* 

it 

it 

it 

* 

* 

It 

it 

it 

it 

it 

it 

* 

* 

* 

* 

★ 

★ 

★ 

* 

* 

★ 

* 

★ 

* 


***************************************************************1 
/include <stdio.h> 

/include <math.h> 

/include "model. h" 
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* PURPOSE: 

* INPUTS: 

* OUTPUTS: 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

★ 


* 

* 

* 

* 

* 

it 

it 

it 

it 

it 

it 

it 

it 

* 

* 

★ 

it 

it 


^*************************************************************** 

* FUNCTION: read_model * 

Write data to C02REMOV.DAT for re-starting model * 
none. 

time — current time in half cycle 
xMg — mass of gas in void space 
xMd — mass of gas in sorbent material 
xTb — temperature of bed and gas desorbed 
xUa — internal energy in C02 accumulator 
xMa — mass of C02 in accumulator 
xT_sorb — temperature of sorbent bed 
xT_des — temperature of desiccant bed 
xMd_w — mass of H20 in desiccant bed 
xTa — temperature of accumulator 
xTp — temperature of pump outlet 
xPp — pressure of pump outlet 
w — pump speed 

power — power applied to sorbent bed heaters 
pair — which pair of beds is operating 
Cci — incoming concentration of C02 
Cci_w — incoming concentration of H20 

***************************************************************/ 
void read_model ( ) 

{ 

model = fopen( "co2remov.dat" , "r" ) ; 

fscanf (model, "%d %f %f %f %f %f", pair, time, xMg,xMd,xMd_w, xMa) ; 
fscanf (model, "%f %f %f %f %f %f",xTb,xT_sorb,xT_des,xUa,w, Power) 
fscanf (model, "%f %f %f %f %f M ,xTa,xTp,xPp,Cci,Cci_w) ; 
f close (model) ; 

> 

y*************************************** ****************** -kit it it it it 

* FUNCTION: read_commands * 

Read CLIPS C02REM0V.CMD file for model operation * 
instructions. CLIPS commands are currently * 

ignored. * 

none . * 

mode — current mode of operation as determined * 

by CLIPS * 

*^******#***************************************Hf**************^ 

void read_commands ( ) 

{ 

clips = fopen( "co2remov. cmd" , "r" ) ; 
fscanf (clips, "%s",mode) ; 
f close ( clips) ; 

> 

^*************************************************************** 

* FUNCTION: decode_commands * 

PURPOSE: Decode the command issued by CLIPS into model * 

parameters. Since the CLIPS commands are ignored * 
no decoding is done. * 

INPUTS: mode — current mode of operation as determined * 

by CLIPS * 

OUTPUTS: tstop — end of half cycle time * 

power — power applied to sorbent bed heaters * 

m_out_scale — fraction of mass flow into C02 * 

accumulator released to C02 * 

reduction assembly * 

w — C02 pump speed * 

blow — blower speed * 

*************************************************************** j 

void decode commands ( ) 

{ 


PURPOSE : 


* 

* 

★ 

* INPUTS: 

* OUTPUTS: 

* 
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> 


tstop = 55.0; 


/** model runs on 55 minute half cycles **/ 


y***************** ******************************************** 

* FUNCTION: adsorb_l * 

* PURPOSE: Simulate the operation of Desiccant bed #1, the * 

* blower, heat exchanger and sorbent bed #2. * 

* This provides for later simulation of component * 

* failures. * 

*************************************************************** i 

void adsorb_l() 

{ 

Ccs_w = xMd_w /Mb; 
k_w = 0.1; 

Cco_w = k_w * Ccs_w; 
xConc_w = Cco_w * 10000.0; 

Ccs = xMd / Mb; 
k_s = 0.1; 

Cco = k_s * Ccs; 
xConc = Cco * 1.0E6; 
dt = 0.1; 

for (dummy=l;dumray<=5*60/dt;dummy ++) { 

/** desiccant bed equations **/ 
dT_des = dt * ( (M_f low * Cpa * (Ti - xT_des)) + 

(M_flow * Sc_w * (Cci_w - Cco_w))) / (Mb * Cvs); 
dMd_w = M_f low * (Cci_w - Cco_w) * dt; 
xT_des = xT_des + dT_des; 
xMd_w = xMd_w + dMd_w; 

Ccs_w = xMd_w /Mb; 

Cco_w = k_w * Ccs_w; 
xConc_w = Cco_w * 10000.0; 

/** fan and heat exchanger "equations" **/ 
xT_fan = xT_des + 6.81; 

xT_hx = xT_des - 50.0; /** this is a kludge **/ 

/** sorbent bed calculations **/ 
dT_sorb = ( (M_flow * Cpa * (xT_hx - xT_sorb) ) + 

(M_f low * Sc * (Cci - Cco))) * dt / (Mb * Cvs); 
dMd = M_f low * (Cci - Cco) * dt; 
xT_sorb = xT sorb + dT_sorb; 
xMd = xMd + ctMd; 

Ccs = xMd / Mb; 

Cco = k_s * Ccs; 

xConc = Cco * 1.0E6; 

xTime = xTime + dt/60.0; 


> 


} 


/*************************************************************** 

* FUNCTION: desorb_l * 

* PURPOSE: Simulate the operation of sorbent bed #4, C02 * 

* pump and C02 accumulator. * 

* This provides for later simulation of component * 

* failures. * 

****************************************** ********************* ^ 

void desorb_l ( ) 

{ 

for (dummy=l;dummy<=5*60/dt;dummy++) { 

if (xTb > 330.0) /** bang-bang bed temp, controller **/ 

Power =0. 0; 
else 

Power = 1000.0; 

Pg = xMg * R * xTb / Vg; 

Pb = kl * xMd /Mb * (xTb - Tref); 
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xPp = xMa * R * xTa / Va; 

xTp = xTb*pow( (xPp/Pg) , (1.0-1. 0/k) ) ; 

xTa = xUa / (xMa * Cp) + To; 

m_o = w*k3*Pg/R/xTb; 

if (xMa < 0.002) /** bang-bang flow out control **/ 

m_out = 0.0; 
else 

m_out = 0.75 * m_o; 

if (xTa > 600.0) /** bang-bang accumulator cooler **/ 

Q_loss=k4* (xTa - 300); 
else 

Q_loss=0.0; 

dMd = (Pg - Pb) * k2 * dt; 

dMg = ((Pb - Pg) * k2 - (w * Pg * k3)/R/xTb) * dt; 
dTb = (dMd * Sc / dt + Power) /Mb/Cvb * dt; 
dMa = dt * (m_o - m_out) ; 

dUa = dt * (m_o* (xTp-To) *Cp-m_out* (xTa-To) *Cp-Q_loss) ; 

xMg = xMg + dMg; 

xMd = xMd + dMd; 

xTb = xTb + dTb; 

xMa = xMa + dMa; 

xUa = xUa + dUa; 

xTime = dt/60.0 + xTime; 

> 

> 


* FUNCTION: adsorb_2 * 

* PURPOSE: Simulate the operation of Desiccant bed #3, the * 

* blower, heat exchanger and sorbent bed #4. * 

* This provides for later simulation of component * 

* failures. * 


***************************************************************/ 
void adsorb_2 ( ) 

{ 

Ccs_w = xMd_w /Mb; 
k_w = 0.1; 

Cco_w = k_w * Ccs_w; 
xConc_w = Cco_w * 10000.0; 

Ccs = xMd / Mb; 
k_s = 0.1; 

Cco = k_s * Ccs; 
xConc = Cco * 1.0E6; 
dt = 0.1; 

for (dummy=l;dummy<=5*60/dt;dummy++) { 

/** desiccant bed equations **/ 
dT_des = dt * ( (M_flow * Cpa * (Ti - xT_des)) + 

(M_flow * Sc_w * (Cci_w - Cco_w) ) ) / (Mb * Cvs); 
dMd_w = M_flow * (Cci_w - Cco_w) * dt; 
xT_des = xT_des + dT_des; 
xMd_w = xMd_w + dMd_w; 

Ccs_w = xMd_w / Mb; 

Cco_w = k_w * Ccs_w; 
xConc_w = Cco_w * 10000.0; 

/** fan and heat exchanger "equations" **/ 
xT_fan = xT_des + 6.81; 

xT_hx = xT_des - 50.0; /** this is a kludge **/ 

/** sorbent bed calculations **/ 
dT_sorb = ( (M_flow * Cpa * (xT_hx - xT_sorb)) + 

(M_flow * Sc * (Cci - Cco))) * dt / (Mb * Cvs); 
dMd = M_flow * (Cci - Cco) * dt; 
xT sorb = xT sorb + dT sorb; 
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xMd = xMd +• dMd; 

Ccs = xMd / Mb; 

Cco = k_s * Ccs; 
xConc = Cco * 1.0E6; 
xTime = xTime + dt/60.0; 

> 

> 


* FUNCTION: desorb_2 * 

* PURPOSE: Simulate the operation of sorbent bed #2, C02 * 

* pump and C02 accumulator. * 

* This provides for later simulation of component * 

* failures. * 


*************************************************************** j 

void desorb_2 ( ) 

{ 

for ( dummy=l ; dummy<= 5 * 60/dt ; dummy++ ) < 

if (xTb > 330.0) /** bang-bang bed temp, controller **/ 

Power=0 . 0 ; 

else Power = 1000.0; 

Pg = xMg * R * xTb / Vg; 

Pb = kl * xMd /Mb * (xTb - Tref); 

xPp = xMa * R * xTa / Va; 

xTp = xTb*pow( (xPp/Pg) ,{ 1.0-1. 0/k) ) ; 

xTa = xUa / (xMa * Cp) + To; 

m_o = w*k3*Pg/R/xTb; 

if (xMa < 0.002) /** bang-bang flow out control **/ 

m_out = 0.0; 
else 

m_out = 0.75 * m_o ; 

if (xTa > 600.0) /** bang-bang accumulator cooler **/ 

Q_loss=k4* (xTa - 300); 
else 


Q_loss=0.0; 

dMd = (Pg - Pb) * k2 * dt; 

dMg = ((Pb - Pg) * k2 - (w * Pg * k3)/R/xTb) * dt; 
dTb = (dMd * Sc / dt + Power ) /Mb/Cvb * dt; 
dMa = dt * ( m_o - m_out ) ; 

dUa = dt * (m_o* (xTp-To) *Cp-m_out* (xTa-To) *Cp-Q_loss) ; 

xMg = xMg + dMg; 

xMd = xMd + dMd; 

xTb = xTb + dTb; 

xMa = xMa + dMa; 

xUa = xUa + dUa; 

xTime = dt/60.0 + xTime; 


> 

^*************************************************************** 


* FUNCTION: write_clips 

* PURPOSE: Write~data to C02REM0V . SEN for CLIPS to analyze. 

* INPUTS: time — current time in half cycle 

* Cci — concentration of C02 at assembly inlet 

* CCo — concentration of C02 at assembly outlet 

* xMa — mass of C02 in accumulator 

* m_out — mass C02 released to reduction assembly 

* OUTPUTS: none. 


* 

* 

* 

* 

* 

★ 

* 

* 


void write_clips ( ) 

{ 

clips = fopen( "co2remov. sen" , "w" ) ; 

fprintf (clips, " (sensor C02_REMOVAL INLET_C02 %f kg/kg %f)\n" 
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, Cci,time) ; 

fprintf (clips, "( sensor C02_REM0VAL OUTLET_C02 %f kg/kg %f)\n" 

, Cco,time); 

fprintf( clips, "(sensor C02_REM0VAL C02_ACCUMULATOR %f kg %f)\n" 

, xMa, time) ; 

fprintf{ clips, "(sensor C02_REMOVAL C02_DEMAND %f kg %f)\n” 

,m_out, time) ; 
fclose(clips) ; 

> 

^*************************************************************** 

* FUNCTION: write_model * 

* PURPOSE: Write data to C02REM0V.DAT for re-starting model * 

* INPUTS: time — current time in half cycle * 

* xMg — mass of gas in void space * 

* xMd — mass of gas in sorbent material * 

* xTb — temperature of bed and gas desorbed * 

* xUa — internal energy in C02 accumulator * 

* xMa — mass of C02 in accumulator * 

* xT_sorb — temperature of sorbent bed * 

* xT_des — temperature of desiccant bed * 

* xMd_w — mass of H20 in desiccant bed * 

* xTa — temperature of accumulator * 

* xTp — temperature of pump outlet * 

* xPp — pressure of pump outlet * 

* w — pump speed * 

* power — power applied to sorbent bed heaters * 

* pair — which pair of beds is operating * 

* Cci — incoming concentration of CC2 * 

* Cci_w — incoming concentration of C02 * 

* OUTPUTS: none. * 

*************************************************************** j 

void write_model ( ) 

{ 

model = fopen( "co2remov.dat" , "w" ) ; 

fprintf (model, "%d %f %f %f %f %f \n", pair, time, xMg,xMd,xMd_w, xMa) ; 
fprintf (model, "%f %f %f %f %f %f\n", xTb,xT_sorb,xT_des,xUa,w, Power) ; 
fprintf (model, "%f %f %f %f %f ",xTa,xTp,xPp,Cci,Cci_w) ; 
f close (model) ; 

> 

/*************************************************************** 

* The two halves of the assembly could be modeled using only * 

* two functions, but by using four functions, the model can * 

* by used to simulate component failures and valve change * 

* sequencing. * 

***************************************************************/ 
main( ) 

{ 

read_model ( ) ; /** get info to start model **/ 

read_commands ( ) ; /** read CLIPS command file **/ 

print f ( "\nOperating in %s mode .", mode ) ; 

decode_commands ( ) ; /** decode CLIPS command into model 

parameters **/ 


* Entry point for restarting model. * 

* * 

* First check to see if mode change should be made, and * 

* make if necessary. Then using model data, run model * 

* for five minutes, writing sensor data and model data * 

* and quit to batch file. * 

***********************************************************| 
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if (time >= tstop) { 
time = 0.0; 
if (pair == 1) 
pair = 2; 
else 

pair = 1; 


/** Determine which pair of functions 
/** will be operating. If the model 
/** is about to change half cycles, 
/** make the necessary arrangements. 


> 

printf ( "\nPair %d is operating. pair) ; 


it it j 

* * J 
it it J 

•kltj 


> 


if (pair == 1) { 
adsorb__l ( ) ; 
desorb~l ( ) ; 
write_clips ( ) ; 
write_model ( ) ; 

> 


else { 

adsorb_2 ( ) ; 
desorb_2 ( ) ; 
write_clips( ) ; 
write_model ( ) ; 


return NORMAL; 


/** Run the model for five minutes and **/ 
/** output data to the necessary files.**/ 


/** Run the model for five minutes and **/ 
/** output data to the necessary files.**/ 


/** main() should return a value to DOS **/ 
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APPENDIX 2B 

C0 2 REDUCTION ASSEMBLY PROGRAM 


/ 


******************************************************* 


* 

* Program Name 

* Author : 

* Date: 

* Description: 


* 

Sabatier Methanation Reactor * 

Jeff Newell * 

2/15/91 * 

This program will determine and * 
set a reactor volume depending on * 
the initial carbon dioxide flow * 

rate entered by the user. It will * 
then let the user enter different * 
flow rates to see how conversion * 

is affected. It will also give * 

outputs of C02 and H2 not reacted * 
as well as H20 and CH4 produced. * 
******************************************************** ^ 


h 


INCLUDE THE NECESSARY HEADER FILES 


7 


#include <stdio.h> 
#include <math.h> 


/' 


DEFINE CONSTANTS USED 


7 


#define A 1.112393E6 /* catalyst activity constant */ 

/define E 8553 /* catalyst activity constant = Ea/R */ 


/’ 


START MAIN FUNCTION 


7 


main ( ) 

{ 

/*- 

int 


INITIALIZE VARIABLES */ 

doanother=l, /* variable for doing another run */ 
index; /* index for arrays */ 


double 


double 


actconv , 

actf low, 

ch4prod, 

co2notreac, 

co2rate, 

conversion, 

h2notreac, 

h2oprod, 

integral , 

keq, 

partlkeq, 

part2keq, 

partlreac, 

part2reac, 

part3reac, 

sqrtemp, 

value, 

volcentl, 

volcent2 , 

volgal ; 

reac[ 550 ] > 
pco2 [550] , 
pch4[550] , 


/* 

/* 

/* 


/* actual conversion rate i 
/* actual flow rate of C02 
/* CH4 produced */ 

C02 not reacted */ 
incoming flow of C02 */ 
conversion rate */ 

/* H2 not reacted */ 

/* H20 produced */ 

/* mass balance integral */ 

/* equilibrium constant */ 
/* part of keq */ 
part of keq */ 
part of reaction rate 
part of reaction rate 
part of reaction rate 
square of temperature 
used for integration *- 


/ 

*/ 


/* 

/* 

/* 

/* 

/* 

/* 

/* 

/* 

/* 


*/ 

*/ 

*/ 

*/ 

7 


volume of sabvol in cubic cent. */ 
volume of sabatier in cubic cnt */ 
volume of Sabatier in gallons */ 


/* rate of reaction */ 

/* partial pressure of C02 */ 
/* partial pressure of CH4 */ 
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ph2o[ 550 ] , 
ph2 [ 550 ] , 
sum[550] , 
temp [ 550} ; 


/* partial pressure of H20 */ 

/* partial pressure of H2 */ 
/* integral */ 

/* temperature of reaction */ 


/* ENTER INPUT VALUES */ 

/* This section will let the user input a flow rate of C02 * / 

printf ( "Enter initial flow rate of carbon dioxide (in mol/sec): 
scanf (”%lf",&co2rate) ; 

/* This section calculates the volume required to achieve 
maximum conversion in cubic cent, and gallons */ 

volcentl = 1.652E7 * co2rate; 
volgal = volcentl * 2.642E-4; 

/* This section lets the user know the volume necessary 
to achieve maximum conversion for his flow rate */ 

printf ("In order to achieve a maximum conversion rate of 0.515, 
printf ("the volume of the reactor must be at least %5.21f 
\n" , volcentl) ; 

printf ( "cubic centimeters, or in other words, %2.21f gallons. 
\n" , volgal) ; 

printf ("\n\n\n") ; 

/* This section sets the volume of the reactor to the value 
above * / 

printf ("The volume of the reactor is now set at %2.21f gallons. 
\n", volgal) ; 

printf ("\n\n") ; 


while ( doanother— 1 ) 

{ 

/* This section lets the user input his own volume */ 

printf ("You may now enter a test flow rate of carbon \n"); 
printf ( "dioxide to see the exact conversion rate for \n"); 
printf ("the volume set above. \n"); 
printf ( "\n" ) ; 

printf ( "Enter test carbon dioxide flow rate (in mol/sec): 

\n"); 

scanf ( "%lf " , fiactf low) ; 

/* This section determines what the value of the mass balance 
integral must be for the user's volume and flow rate */ 

integral = volcentl / actflow; 

/* This section compares the value of the mass balance 
integral to the value it should be for maximum conversion. If the integral is 
larger than 1.652E7, then the volume of the reactor is large enough to achieve 
maximum conversion for the user's flow rate. If this is the case, this 
section sets the actual conversion rate to 51.5% and determines the temperature 
of reaction. */ 


\n” ) ; 


\n" ) ; 
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if ( integral>=l . 652E7 ) 

{ 

actconv = .515; 

index = 515; 

temp [index] = 1088 * actconv + 297.59; 

> 

/* This section also compares the value of the mass balance integral to the 
maximum value. In this case, if the integral is smaller than 1.652E7, then the 
reactor is not large enough to achieve maximum conversion. If this is the case, 
this section evaluates the mass balance integral to determine the conversion rate 
for the inputted flow rate of C02. */ 

if ( integraKl. 652E7 ) 

{ 


/* Initialize values */ 

conversion = 0.0; 
index =0; 
actconv = 0.0; 
sqrtemp = 0.0; 
value = 0.0; 

for ( conversion* . 001 ; conversions .515; conversion** . 001 ) 

{ 

index = conversion * 1000; 

temp[index] = 1088 * conversion + 297.59; 

/* partl&2keq used to calculate keq */ 

sqrtemp = temp [index]* temp [ index ] ; 

partlkeq = 5600/sqrtemp + 34633/temp[ index] ; 

part2keq = -16.4* (log(temp[ index] ) ) + 0. 005 57 *temp[ index] ; 

keq * exp( (1.0/1. 987 )* (partlkeq + part2keq) + 33.165); 

/* partial pressures of reactants and products */ 
pco2[index] = 0.2*( ( 1-conversion) /( 1-. 4*conversion) ) ; 
pch4[ index] = 0. 2* ( conversion/ ( 1-0. 4*conversion) ) ; 
ph2o[ index] = 0.2*( (2*conversion)/(l-0.4*conversion) ) ; 
ph2[ index] = 0.2* ( (4-4*conversion) / ( l-0.4*conversion) ) ; 

/* parti, 2&3reac used to calculate reaction rate */ 
partlreac = pow(pco2 [ index] ,0.25) *ph2 [ index] ; 
part2reac = pow(pch4 [ index] , 0. 25 ) *pow(ph2o[ index] , 0. 5 ) /keq; 
part3reac = A/temp [ index] ; 

reac[index] = partireac* (partlreac-part2reac ) *exp(-E/temp[ index ] ) ; 

/* This line multiplies the integrand by the time step */ 

value = ( l/reac[ index ])*. 001; 

/* This line adds the values above to determine conversion */ 
sum [index ] * sum [index-1] + value; 

/* determines if integral has been computed */ 
if ( sum [ index ] >=integral ) 

{ 

actconv=conversion; 
conversion* . 515 ; 

> 

> 

> 

/* This section prints output values */ 

printf ( "Actual rate of conversion = %4.41f \n\n" , actconv) ; 
printf ( "Final temperature of reaction = %6.21f K 
\n\n", temp [index] ) ; 

co2notreac = 44. 0*actf low* ( 1-actconv) ; 

printf ( "Amount of C02 not reacted = %4.41f grams/sec. 
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\n" , co2notreac ) ; 


h2notreac = 2.0*4.0*actflow* (1-actconv) ; 

printf ("Amount of H2 not reacted = %4.41f grams/sec. 

\n" , h2notreac ) ; 

h2oprod = 18 . 0*2 . 0*actf low*actconv; 

printf ( "Amount of H20 produced = %4.41f grams/sec. \n" 
,h2oprod) ; 

ch4prod = 16.0*actflow*actconv; 

printf ( "Amount of CH4 produced = %4.41f grams/sec. 

\n" , ch4prod) ; 

printf("If you wish to do another analysis, enter 1, else 0. 
\n" ) ; 

scanf ( "%d" , fidoanother) ; 

> 

> 

/* END OF FILE */ 
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APPENDIX 2C 

OXYGEN GENREATION ASSEMBLY PROGRAM 


****************************************************************** 

* SOURCE FILE: OGA.C 

* 

* DESCRIPTION: This program models the OGA electrolysis unit, 

H20 tanks, and the H2 and 02 output from the 


electrolysis unit. 
main( ) 

Original in fortran by Cark Perkins, David Kley, 
and Frank Lunkwitz. 


* 

* 

* USES: 

* 

* AUTHORS : 

* 

* 

* Rewrite in C by Kevin Forssberg and William Biehl 

* 

* DATE: March 1, 1991 

************************* **************************************** j 


y************* 

* NOTE TO ROBERT YOUNG AND HAL KRAUS: 

* — KEVIN AND BILL COULD NOT FIGURE OUT HOW TO STOP 

* AND RESTART THE PROGRAM EVERY 5 MIN. THEREFORE, 

* THIS PROGRAM DOES NOT HAVE THIS CAPABILITY. 

************ j 

#include <stdio.h> 

#inc lude <mat h . h> 


#def ine 
#def ine 
#def ine 
/define 


FILE 


F 96487 

FULL 35 /* A FULL WATER TANK IS PRESET AT 35 MOLES */ 

EMPTY 3 /* H20 AND A EMPTY TANK IS AT 3 MOLES H20 */ 

X 18 /* X IS THE NUMBER OF CELLS IN THE ELECTROL- 
YSIS UNIT */ 

*outf ile; /* DECLARAITION OF OUTPUT FILE */ 


main( ) 

{ 

int 1=0, /* ELECTROLYSIS UNIT CURRENT */ 

t=0, /* INCREMENTING VARIABLE IN FOR LOOPS */ 

TIME=0, /* TIME INITIAL SET POINT */ 

LEVEL=0 ; / * THIS IS THE DISCRETE INITIAL SETTING OF 
THE WATER TANK WHERE 1 INDICATES FULL 
AND 0 INDICATES EMPTY */ 


double H20UT, /* MOLES OF H2 OUT OF ELECTROLYSIS UNIT*/ 

/* MOLES OF 02 OUT OF ELECTROLYSIS UNIT*/ 

/* INITIAL SETTING OF H2 TANK 
/* INITIAL SETTING OF 02 TANK 
/* AMOUNT OF WATER IN SYSTEM 
/* OPERATING TEMP OF ELECTROLYSIS UNIT 
02_0UTTAKE_TEMP , /* LINEAR CONVERTION FOR TEMP OF 02 
OUT TO THE CABIN */ 


H20UT, 
020UT, 
H2TANK=0 . 0 , 
O2TANK=0 . 0 , 
H20, 

TEMP, 


*/ 

*/ 

*/ 

*/ 


j ************************************************************* j 
/ * PRAMETERS NOT CALCULATED BY MODEL BUT NEEDED FOR EXPERT */ 
/* SYSTEM. — TO VARY EXPERT SYSTEM OUTPUTS CHANGE THESE */ 

/* VARIABLES AT THIS POINT ONLY. */ 
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^ ************************************************************* ^ 
H20 FEEDLINE_PRESSURE=25 .0, /* H20 PRESSURE TO ELECTROL- 
YSIS UNIT */ 

H _0_DEFEREN TI AL_ p RE s SURE=3 . 0 , / * DERERENTIAL PRESSURE OF 

02 AND H2 OUTPUT LINES */ 

CELL_VOLTAGE=1.67; /* INDIVIDUAL CELL VOLTAGE 

IN ELECTROLYSIS UNIT */ 


outf ile=fopen( "OGA.DAT" , "w” ) ; /* OPEN DATAFILE "OGA.DAT” */ 

H20 = FULL; 

for (1=15; I<=28; I++) 

{ 

TIME++; 

H20UT= (double) I/F*X*60.0; 

020UT=H20UT/2 .0; 

H2TANK=H2TANK+H20UT ; 

02 TANK=02 TANK+02 OUT ; 

H20=H20-H20UT-020UT ; 

TEMP=( (-0.0391*1*1) + (4.38*1) + 67); 

02_OUTTAKE_TEMP=TEMP-68 ; 
if (H20 <= EMPTY) 

{ 

H20=FULL; 

LEVEL=0 ; 

> 

else 

LEVEL=1 ; 

fprintf( out file, "(SENSOR OGA H20_FEEDLINE_PRESSURE %4.2f PSIG %2.0d)\n", 
H20_FEEDLINE_PRESSURE , TIME ) ; 

fprintf (outfile,” (SENSOR OGA H20_PHOTOCELL %d discrete %2.0d)\n", 

LEVEL, TIME) ; 

fprintf (out file, "(SENSOR OGA H/0_DEFERENTIAL PRESSURE %4.2f PSIG %2.0d)\n", 
H _°_DEFERENTIAL_PRESSURE,TIME) ; 

fprintf (out file, "(SENSOR OGA 02_OUTTAKE_TEMP %5.3f F %2.0d)\n", 
02_OUTTAKE_TEMP , TIME ) ; 

fprintf (out file, "(SENSOR OGA CELL_VOLTAGE %4.2f volts %2 . Od) \n\n\n" , 
CELL_VOLTAGE , TIME ) ; 

> 

for(t=l; t<200; t++) 

{ 

TIME++; 

H20UT=H20UT ; 

020UT=020UT; 

H2 TANK=H2 TANK+H2 OUT ; 

02 TANK=02 TANK+02 OUT ; 

H20=H20-H20UT-020UT ; 
if (H20 <= EMPTY) 

{ 

H20 = FULL; 

LEVEL = 0; 

> 

else 

LEVEL = 1; 

fprintf (outfile, "(SENSOR OGA H20_FEEDLINE_PRESSURE %4.2f PSIG %2.0d)\n", 
H20_FBEDLINE_PRESSURE , TIME ) ; 

fprintf (outfile, "(SENSOR OGA H20_PHOTOCELL %d discrete %2.0d)\n", 

LEVEL, TIME) ; 

fprintf (outfile, " (SENSOR OGA H/0_DEFERENTIAL PRESSURE %4.2f PSIG %2.0d)\n", 
H O DEFERENTIAL PRESSURE, TIME) ; 

fprintf (outfile, " (SENSOR OGA 02_0UTTAXE_TEMP %5.3f F %2.0d)\n". 
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02_0UTTAX£_TEMP , TIME ) ; 

fprintf (outfile, " (SENSOR OGA CELL_VOLTAGE %4.2f volts %2 .Od) \n\n\n” , 
CELL_VOLTAGE , TIME ) ; 

> 

fclose (outfile) ; /* CLOSE OUTPUT FILE "OGA. DAT" */ 

> 
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APPENDIX 2D 

AIR REVITALIZATION SYSTEM PROGRAM 


Contained within this appendix is information about the compliling of the 
computer modeling groups into one. Information is given for the benefit of 
giving next semester students background on the computer program and why it 
operates the way it does. Also, included in this appendix is logon procedures 
for operating the program on the AT&t lab computers on second floor of 
Durland. 

Logon on Procedures 

To log onto the 3bl's, AT&T lab — second floor Durland: 

login: <student> type student @ login command and return 

use mouse to select Full screen unix 

To run program: insert floppy w/ program 

program name is AR_model.c 

To look at program — Type: mread -t ar_model.c ar_model.c 
To compile program — Type: cc -G -I/usr/ local/ include arjmodel.c -lksu - 
levdi -lm -o run 

To run program Type: run 

3bl's packet, available from Ruth Dyer, gives editor commands, logon 
procedures, etc. 


Assumptions and Conditions of Program 

Assumptions of each individual portion of this modeling program can be 
examined in detail from previous sections. Each model has its’ input and 
output parameters, which are discussed in their individual sections. 

This program takes all three modeling programs and puts them into one working 
program, each dependent on the other. A few things have been changed and 
added to complete this task. 

1. The C02 is not delivered from the Removal to the Reduction until the 
sorbent bed is reversed and C02 is being removed from the bed. 

2. Certain values are inputs to the system and are asked for when running 
the program. 

3. Between plots of output, any key has to be hit to continue to see next 
plot. 


4. A tank is inserted between the C02 Reduction and the OGA systems. This 

tank is initially set and has two inputs and one output: in — C02 Red. and 

Hygiene Water 

out — OGA system. 

5. The OGA system operates a little differently in this program than it 
does seperately. 

a. The system is already at steady state by the time the C02 Red. 
receives any C02 from the C02 Rem. 

b. The time interval is different in this program, it is made to 
better match the time scale that is seen 

from the C02 Removal. 

c. Since the small H20 tanks in the OGA system will not empty in 

the time period in which this program operates, the water is assumed 

to come directly from the tank that we created between the C02 Red. 
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and the OGA. 


Modeling Program 

The following is a listing of the file that was created and edited to fit all 
three working computer models of the AR system. 


/*************************************************************** 
* * 


* 

* 

* 

it 

it 

it 

it 

it 

* 

it 

it 

it 

it 

it 

it 

it 

it 

it 

it 


SOURCE FILE: MODEL4.C * 

* 

DESCRIPTION: This program models the C02 reduction system * 

during its sorbtion and desorption half * 

cycles. C02 conversion and 02 generation * 
by the REDUCTION and OGA subsystems are also * 
included; however, the REDUCTION model is * 
not time dependent. We have arranged it so * 


its outputs appear time dependent by assuming* 


steady-state operation at all times. * 

* 

USES: main() * 

simple_plot( ) (S. Dyer) * 

* 

AUTHORS: Robert Young, Terry Hon, Joel Sommer, * 

Carl Perkins * 

it 

LAST REVISION: 29 April 1991 * 
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#include <stdio.h> 
#include <math.h> 
#include "p plot.h" 
#include "sImple_plot.h" 
/define NUM POINTS 100 


void input(); 
void co2removal( ) ; 
void co2 reduct ion{ ) ; 
void oga( ) ; 
void h2otank ( ) ; 


int dummy; 
long int i,q; 
double Time [100], 
xTime, 

Mg[ 100] , 
xMg, 

Md[ 100 ] , 
xMd, 

Tb( 100] , 
xTb, 

Ua[ 100] , 
xUa, 

Ma[ 100] , 
xMa, 

Pp[100], 

xPp, 

Tp[100] , 
xTp, 


/** time array **/ 

/** mass of gas in void space **/ 

/** mass of gas in sorbent material **/ 

/** temperature of bed and gas desorbed **/ 
/** internal energy of accumulator **/ 

/** mass of C02 in accumulator **/ 

/** output pressure of pump **/ 

/** output temeprature of pump **/ 
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Ta(100] , 
xTa, 

! ** 

Conc_w[ 100] , 
xConc_w, 

J -k-k 

Ccs_w, 

/kk 

M flow. 

j kk 

Tl, 

jkk 

T_des ( 100 ] , 
xT des. 

f kk 

Ccl_w, 

Jkk 

Cco w. 

Jkk 

Md_w(100] , 
xMd_w , 

Jkk 

Ccs, 

Jkk 

Cci, 

Jkk 

Cco, 

Jkk 

Cone [ 100 ] , 
xConc , 

Jkk 

T sorb (100], 
xT sorb. 

jkk 

Fa ( 100] , 

jkk 

F_water[100] ; 

jkk 

pg. 

Jkk 

Pb, 

Jkk 

m_o, 
xm out. 

jkk 

m out ( 100 ] , 

Jkk 

w. 

Jkk 

Q_loss, 

Jkk 

Tref , 

jkk 

k w. 

jkk 

k_s. 

Jkk 

to_red. 

Jkk 

conv. 

Jkk 

F, 

Jkk 

X, 

Jkk 

H20UT [ 100 ] , 

jkk 

020UT[ 100 ] , 

Jkk 

a. 

jkk 

minimum. 

jkk 

hyg H20, 

Jkk 

oga_H20[ 100 ] , 

jkk 

elapse [100] , 

Jkk 

level (100] ; 

jkk 


temperature of accumulator **/ 

relative humidity of outgoing air **/ 

cone, of H20 on sorbent (kgH20/kgair) **/ 
mass flow rate of air (kg/sec) **/ 
incoming air temperature (K) **/ 
temp, of air leaving desicant bed (K) **/ 

cone, of H20 entering bed (kgH20/kgair ) **/ 
cone, of H20 leaving bed (kgH20/kgair ) **/ 
mass of H20 on sorbent (kg) **/ 

cone. C02 in sorbent (kgC02/kg air) **/ 
cone. C02 entering bed (kg C02/kg air) **/ 
cone. C02 leaving bed (kg C02/ kg air) **/ 
cone. C02 leaving bed (ppm) **/ 

temp, or air leaveing sorbent bed (K) **/ 

C02 flow exiting the REDUCTION ass. (kmoles/s) **/ 
H20 flow exiting REDUCTION ass. (kmoles/s)**/ 


C02 pressure in void space **/ 
pump exit pressure **/ 

mass flow rate out of bed dMg/dt **/ 

assembly exiting flow rate (bang-bang) **/ 
angular velocity of pump (rad/sec) **/ 
accumulator heat loss (bang-bang) **/ 
reference temp, of desorbent bed **/ 
ratio of Cco w/Ccs_w **/ 
ratio of Ccojccs **/ 

C02 flow rate to REDUCTION assembly (kg/s) **/ 
conversion factor for REDUCTION **/ 

Faraday's constant for OGA **/ 
number of electrolysis cells in OGA **/ 

OGA hydrogen production rate **/ 

OGA oxygen production rate **/ 

OGA current in amps **/ 
minimum level in h2otank **/ 

H20 flow rate into h2o_tank from HYGIENE **/ 
H20 flow rate into h2o_tank from OGA **/ 
time array for OGA and h2otank models **/ 
level in h2otank (kmoles)**/ 


double Power= 1000.0, 
Mb=20.0, 
Cpa=1003.0, 
Cvs=844. 0, 

Sc w=2700.0E3, 


/** power input to desorbing C02 bed **/ 
/** mass of sorbent material (kg) **/ 

/** heat capacity of air (J/kgK) **/ 

/** heat capacity of sorbent (J/kgK) **/ 
/** heat of sorption of H20 

into sorbent (J/kgK) **/ 


Cvb=657 .0, 

Cp = 846.0, 

R=0. 1889, 

To=0.0, /** reference temperature for internal energy **/ 

Vg=l .0, 

Va=2 .0, 
kl=0. 7 , 
k2=.0006, 
k3=.0023, 

k4= 2.0, /** accumulator heat transfer coef. (Watts/K) 

k = 1.29, 

Sc=951E3; /** C02 heat of sorption into sorbent (J/kgK) 


** j 
it ic j 
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double dt, 
dTb, 
dT_des , 
dT_sorb , 
dMd_w, 
dHg, 
dMd, 
dMa, 
dUa; 

main ( ) 

{ 

input(); /** subroutine which asks for user input **/ 

co2removal( ) ; /** subroutine modelling the REMOVAL assembly **/ 

co2 reduct ion( ) ; /** subroutine modelling the REDUCTION assembly **/ 

oga(); /** subroutine modelling the OGA assembly **/ 

h2otank(); /** subroutine modelling the H20 tank serving OGA **/ 

return ( 0 ) ; 

> 

/** Inputting of parameters by user **/ 

void input () 

{ 

system ( "clear" ) ; 

print f ( "\nlnput mass flow rate of air, M_flow (kg/s) : "); 
scanf ( "%lf " , &M_f low) ; 

printf ( "\nlnput H20 concentration entering assembly,"); 
printf ( " Cci_w (kgH20/kg air) : "); 
scanf ( "%lf " ,&Cci_w) ; 

printf ( "\nlnput concentration of C02 entering assembly,"); 
printf (" Cci (kgC02/kg air): "); 
scanf ( "%lf " , &Cci) ; 

printf ( "\nlnput temperature of incoming air, Ti (K) : "); 
scanf ("%lf",&Ti) ; 

printf ("\nlnput angular velocity of pump (radians/second) : "); 
scanf ( "%lf " , &w) ; 

printf ("\nlnput desired bed reference temperature, Tref (K) : "); 
scanf ( "%lf " , ScTref ) ; 

printf ( "\nlnput desired conversion of C02 (0.52 max) : "); 
scanf ( "%lf " , Sconv) ; 

} 

void co2removal() 

{ 

Time[0] = 0.0; /** absolute start time in minutes **/ 

xTime = 0.0; 

/** initial conditions for desiccant bed **/ 

Md_w [ 0 ] = xMd_w =0; 

T_des[0] = xT_des = Ti; 

Ccs_w = Md_w(0] / Mb; 
k_w = 0.1; 

Cco_w = k_w * Ccs_w; 

Conc_w[0] = xConc_w = Cco w * 10000.0; 

/** initial conditions for sorbent bed **/ 

Md ( 0 ) = xMd = 0.0, 

T_sorb[0] = xT_sorb = Ti; 

Ccs = Md[0 ] / Mb; 
k_s =0.1; 

Cco = k_s * Ccs; 

Conc[0] = xConc = Cco * 1.0E6; 

dt = 1.0; 

printf ( "\n\nWorking desiccant and sorbent bed forward cycle. .. \n\n" ) ; 

^************************************************* 

* * 

* Desiccant bed and Sorbent bed calculations * 

★ * 

*************************************************1 

i = l; 
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while (i<100) { 

for ( dummy =1 ; dummy<=33 ; dummy++ ) { 

/** desiccant bed equations **/ 
dT_des = dt * ((M_flow * Cpa * (Ti - xT_des)) + 

(M_flow * Sc_w * (Cci_w - Cco_w) ) ) / (Mb * Cvs); 
dMd_w = M_flow * (Cci_w - Cco_w) * dt; 
xT_des = xT_des + dT_des; 
xMd_w = xMd_w + dMd_w ; 

Ccs_w = xMd_w /Mb; 

Cco_w = k_w * Ccs_w; 
xConc_w = Cco_w * 10000.0; 

/** sorbent bed calculations **/ 
dT_sorb = ( (M_f low * Cpa * (xT_des - xT_sorb)) + 

(M_f low * Sc * (Cci - Cco) ) ) * dt / (Mb * Cvs); 
dMd = M_f low * (Cci - Cco) * dt; 
xT_sorb = xT_sorb + dT_sorb; 
xMd = xMd + dMd; 

Ccs = xMd /Mb; 

Cco = k_s * Ccs; 
xConc = Cco * 1.0E6; 
xTime = xTime + dt/60.0; 

} 

Time[i] = xTime; 

T_des ( i ) = xT_des ; 

Md_w [ i ] = xMd_w ; 

Conc_w [ i ] = xConc_w ; 

T_sorb ( i ] = xT_sorb ; 

Md[i] = xMd; 

Cone ( i ] = xConc ; 

i - i + 1; 


/***************************************************** 

* * 

* Plotting routine for Desiccant and Sorbent beds * 

* * 

*******★************★★**★*★***★★**★**★★** ************ ^ 
print f ( ” \n\nDone ..."); 

printf ( "\n\nPress a number and <Return> to plot..."); 
scanf ("%d",&dummy) ; 

/** Simple plot calls for desiccant bed data **/ 
simple_plot (DISPLAY, i, Time, T_des, "Time", "minutes" , "Temperature" , 
"degrees Kelvin", 

"Temperature of Desicant Bed (T_des[ i] )", 1, CURVE) ; 
scanf ("%d",&dummy) ; 

simple_plot (DISPLAY, i, Time, Md_w, "Time" , "minutes" , "Mass H20” , "kg" , 
"Mass H20 in Sorbent Bed Material (Md_w[ i] ) " , 1, CURVE) ; 
scanf ( "%d" , &dummy ) ; 

simple_plot (DISPLAY, i,Time,Conc_w, "Time" , "minutes" , 

"Relative Humidity" , 

"Relative Humidity of Air Leaving Desicant Bed", 

1 , CURVE ) ; 

scanf (”%d",&dummy) ; 

/** Simple plot calls for sorbent bed data **/ 
simple_plot (DISPLAY, i,Time,T_sorb, "Time" , "minutes" , "Temperature" , 
"degrees Kelvin" , "Temp, of Air Leaving Sorbent Bed", 

1, CURVE) ; 

scanf ("%d",&dummy) ; 

simple_plot (DISPLAY, i,Time,Md, "Time" , "Minutes" , "Mass C02 " , "kg” , 
"Mass C02 in Sorbent Bed Material (Md[i])",l,CURVE); 
scanf ("%d",&duromy) ; 

simple_plot (DISPLAY, i, Time, Cone, "Time”, "minutes", "Concentration", 
"ppm", "Cone. C02 Leaving Sorbent Bed", 1, CURVE ) ; 
scanf ( "%d" ,&dummy) ; 

system ( "clear" ) ; /** clear the screen **/ 
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print f ( "\n\nWorking on desorbing cycle. . . " ) ; 

/** initial conditions for desorbing cycle **/ 
Time[0] = xTime = 62; /** reset array for this cycle **/ 

Md[0] = xMd = Md[99}; /** reset array for this cycle **/ 

Tb[0] = xTb = T_sorb(99]; /** reset array for this cycle **/ 
Mg[0] = kl * Vg * Md[0] * (Tb[0) - Tref) / (Mb * R * Tb[OJ); 

xMg = Mg ( 0 ] ; 

Ua(03 = xUa = 407.34; 

Ma(0) = xMa = 0.002; 

Ta[0] = xTa = 300.0; 

Pp[0] = xPp = 0.0567; 

Tp[0] = xTp = 300.0; 

/it********************************************* 
* * 

* Calculations for Desorbing cycle * 

* * 

********************************************** y 

i = l; 

while (i<100) { 

for ( dummy =1 ; dummy<=28 ; dummy++ ) < 

if (xTb > 330.0) /** bang-bang bed temp, controller **/ 

Power=0 . 0 ; 

else Power = 1000.0; 

Pg = xMg * R * xTb / Vg; 

Pb = kl * xMd /Mb * (xTb - Tref); 

xPp = xMa * R * xTa / Va; 

xTp = xTb*pow( (xPp/Pg) , (1.0-1. 0/k) ) ; 

xTa = xUa / (xMa * Cp) + To; 

m_o = w*k3*Pg/R/xTb; 

if (xMa < Ma [ 0 ] ) /** bang-bang flow out control **/ 

xm_out = 0.0; 
else 

xm_out =0.85 * m_o; /** C02 flow rate to REDUCTION **/ 

if (xTa > 860.0) /** bang-bang accumulator cooler **/ 

Q_loss=k4* (xTa - 300); 
else Q_1 os8=0.0; 

dMd = (Pg - Pb) * k2 * dt; 

dMg = ((Pb - Pg) * k2 - (w * Pg * k3)/R/xTb) * dt; 

dTb = (dMd * Sc / dt + Power ) /Mb/Cvb * dt; 

dMa = dt * (m_o - xm_out); 

dUa = dt * (m_o* (Tp[ i-l]-To) *Cp-xm_out* 

(Ta( i-l)-To) *Cp - Q_loss); 

xMg = xMg + dMg; 
xMd = xMd + dMd; 
xTb = xTb + dTb; 
xMa = xMa + dMa; 
xUa = xUa + dUa; 
xTime = dt/60.0 + xTime; 

> 

Time[ ij = xTime; 

Mg[i) = xMg; 

Md[ i] = xMd; 

Tb[ i] = xTb; 

Ma[i] = xMa; 

Ua[ i) = xUa; 

Ta[i] = xTa; 

Tp[ i] = xTp; 

Pp[ i] = xPp; 
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m_out [ i ] = xm_out ; 
i = i + 1; 


printf ( "\nPress* number and return to begin plotting..."); 

scanf ( "%d" , Sdummy ) ; 

/** SIMPLE PLOT CALLS FOR DESORBENT CYCLE DATA **/ 

simple_plot( DISPLAY, i, Time, Ma, "Time", "minutes", 

"Mass C02" , "kg" , "Mass C02 in Accumulator (Ma[i))" 

, 1 , CURVE ) ; 

scanf (”%d”,&dummy) ; 

simple_plot (DISPLAY, i, Time, Ua, "Time" , "minutes" , 

"Internal Energy", "J", 

"Internal Energy in Accumulator (Ua[iJ)" 

, 1, CURVE) ; 

scanf ( "%d" , fidummy) ; 

simple_plot (DISPLAY, i, Time, Tp, "Time" , "minutes" , "Temperature" , 
"K" , "Temperature of Pump Outlet (Tp[i])"» 1, CURVE); 

scanf ("%d",£dummy) ; 

simple_plot (DISPLAY, i,Time,Ta, "Time" , "minutes" , "Temperature", 
"K" , "Temperature of Accumulator (Ta[i] )", 1, CURVE) ; 

scanf ( " %d" , & dummy ) ; 

simple_plot (DISPLAY, i, Time, Pp, "Time", "minutes" , "Pressure" , 

"kPa" , "Outlet Pressure of Pump (Pp[ i) )", 1, CURVE) ; 

scanf ( "%d" , Sdummy ) ; 

printf ( " \nDone sending plots . . " ) ; 

system ( "clear" ) ; 


void co2 reduct ion ( ) 

{ 

for ( i=l ; i<100;i++) 

{ 

Fa(i] = m_out(i] /44.0 * (1-conv); 
F_water[i] =2/44.0 * m_out[i] * conv; 

} 


simple_plot (DISPLAY, i, Time, Fa, "Time" , "minutes" , "C02 Flow Rate", 

"kmoles/s" , "C02 Flow Rate Exiting REDUCTION Assembly (Fa(i])", 
1, CURVE) ; 

scanf ("%d",&dummy) ; 

simple_plot (DISPLAY, i, Time, F_water, "Time" , "minutes" , "H20 Flow Rate" , 
"kmoles/s" , 

"H20 Flow Rate Exiting REDUCTION Assembly (F_water ( i] ) " , 

1, CURVE); 

scanf ( "%d" , fidummy) ; 

printf ("\nDone with REDUCTION. .. \n\n" ) ; 
system( "clear" ) ; 


void oga( ) 

{ 

F=96487 . 0 ; 

X=18.0; 

a=15.0; /** initial current in electrolysis module **/ 

elapse[0]=0.0; /** time array **/ 

H20UT ( 0 ] =0 . 0 ; 

020UT [ 0 ] =0 . 0 ; 
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oga_H20 [ 0 ] =0 . 0 ; 

for (q=l;q<=13;q++) /** this loop models transient operation **/ 

{ /** when the current is increasing w/time **/ 

H20UT(q]= a*X/F; 

020UT[qj= H2OUT[q]/2.0; 
oga_H20[q]= H20UT[q); 
a=a+l; 

> 

while (q<=46) /** this loop models steady-state operation **/ 

{ /** when current is constant ** / 

H20UT[q] =a*X/F; 

020UT [ q ] =H20UT [ q ] /2 . 0 ; 
oga_H20 ( q ] =H20UT [ q ] ; 
elapse[q]=q; 

q=q+i ; 

> 

simple_plot (DISPLAY, q, elapse, 020UT, "Time" , "minutes" , "flow rate" , 
"kmoles/s" , "Oxygen Production by OGA" , 1, CURVE ) ; 
scanf ("%d",&dummy) ; 
system ("clear"); 
printf ( "\nDone with OGA. . .\n\n" ) ; 


void h2otank( ) 

{ /** models the water tank that the OGA draws from **/ 

q=l; 

hyg_H2O=0 . 0 ; 
level[0]=65.0; 
while (q<=45) 

minimum=50 . 0 ; 

if ( level (q]>minimum) /** checks hyg_H20 supply **/ 
hyg_H2O=0.0; 

else hyg_H2O=0 . 1 ; /** turns hyg_H20 supply on **/ 

level [ q] =level ( q-1 ] - ( 60 . 0*oga_H2O [ q ) ) + 

( 60 . 0*F_water ( i ] +hyg_H20 ) ; 

q=q+l; 

} 

simple_plot (DISPLAY, q, elapse, level, "Time" , "minutes” , "tank level" , 
"kmoles" , "h2otank level", 1, CURVE ) ; 
scanf ( "%d" , Sdummy) ; 
system ( "clear" ) ; 

printf ( "\n\nDone with AR-Simulation. . . \n\n" ) ; 
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APPENDIX 3A 

C0 2 REMOVAL ASSEMBLY KNOWLEDGE BASE 


C02 REMOVAL PRECOOLER OUTPUT AIR TEMPERATURE SENSOR (RULE #1) 
PARAMETER MODIFIER DEFINITIONS 

(level: if COOLER TEMP is X then TEMP is HI 100 212 500 900 F) 

(level: if COOLER~TEMP is X then TEMP is NORMAL 0 32 100 212 F) 

(level: if COOLER JTEMP is X then TEMP is LOU -400 -100 0 32 F) 

DEFINE CONDITIONS 

(condition: if COOLER TEMP is HI then COOLER is FAIL 1) 

(condition: if COOLER~TEMP is NORMAL then COOLER is FAIL 1) 

(condition: if COOLER~TEMP is LOW then COOLER is FAIL 1) 

DEFINE RESPONSES 

(response: if COOLER is FAIL then C02_REMOVAL_COMMAND is OFF 0.70) 


C02 REMOVAL BLOWER PRESSURE DIFFERENTIAL SENSORS (RULE #2) 

PARAMETER MODIFIER DEFINITIONS 

(level: if BLOWER is X then BLOW PRES is HI 18 35 100 100 psia) 

(level: if BLOWER is X then BLOWj>RES is NORMAL 8 12 18 35 psia) 

(level: if BLOWER is X then BLOU_PRES is LOU -100 -100 8 12 psia) 

DEFINE CONDITIONS 

(condition: if BLOU_PRES is HI then BLOWER is FAILED 0.2) 

(condition: if BLOU_PRES is NORMAL then BLOWER is FAILED 0.0) 

(condition: if BLOU_PRES is LOW then BLOWER is FAILED 1.0) 

DEFINE RESPONSES 

(response: if BLOWER is FAILED then C02_REMOVAL_COMMAND is OFF 1) 


C02 PUMP POWER SENSOR (RULE #3) 

PARAMETER MODIFIER DEFINITIONS 

(level: if C02_PUMP is X then PUMP POWER is HI 400 600 800 1000 watt) 

(level: if C02 PUMP is X then PUMP~POWER is NORMAL 20 100 300 500 Watt) 

(level: if C02_PUMP is X then PUMP_POUER is LOW 0 5 20 100 watt) 

DEFINE CONDITIONS 

(condition: if PUMP_POUER is HI then C02 PUMP is FAILED 1.0) 

(condition: if PUMP_POUER is NORMAL then C02~PUMP is FAILED 0.0) 

(condition: if PUMP_POUER is LOW then C02~PUMP is FAILED 0.4) 

DEFINE RESPONSES 

(response: if C02_PUHP is FAILED then C02_REMOVAL_COMMAND is OFF 1) 


C02 PUMP PRESSURE SENSOR (RULE #4) 

PARAMETER MODIFIER DEFINITIONS 

(level: if C02_PUMP PRES is X then PUMP PRES is HI 50 100 500 900 psia) 

(level: if C02_PUMP~PRES is X then PUMP~PRES is NORMAL 2 10 30 75 psia) 

(level: if C02_PUMP_PRES is X then PUMP_PRES is LOU 0025 psia) 

DEFINE CONDITIONS 

(condition: if PUMP PRES is HI then C02 PUMP is FAILED 0.8) 

(condition: if PUMP - PRES is NORMAL then C02 - PUMP is FAILED 0.0) 

(condition: if PUMP_PRES is LOW then C02_PUMP is FAILED 0.1) 
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define responses 


(response: if C02_PUMP is FAILED then C02_ 

REMOVAL. 

.COMMAND is 

OFF 

0.9) 


C02 ACCUMULATOR AIR FLOW SENSOR (RULE 

#5) 






PARAMETER MOOIFIER DEFINITIONS 







(level: if C02 ACCUM is X then ACCUM FLOW 

is HI 

50 

100 

500 

900 

Ib/hr) 

(level: if C02~ACCUM is X then ACCUM _ FLOW 

is NORMAL 1 

5 

20 

50 

Ib/hr) 

(level: if C02 ACCUM is X then ACCUM FLOW 

is Low 

-900 

-100 

0 

1 

Ib/hr) 


DEFINE CONDITIONS 


(condition: if ACCUM FLOW is HI then C02 ACCUM is FAILED 0.2) 

(condition: if ACCUM'FLOW is NORMAL then C02~ACCUM is FAILED 0.0) 

(condition: if ACCUM_FLOW is LOW then C02_ACCUM is FAILED 0.8) 

DEFINE RESPONSES 

(response: if C02_ACCUM is FAILED then C02_REMOVAL COMMAND is OFF 0.9) 
(response: if C02_ACCUM is FAILED then C02_REDUCTION_COMMAND is OFF 1.0) 
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APPENDIX 3B 

C0 2 REDUCTION ASSEMBLY KNOWLEDGE BASE 


/** This set of statements defines the sensors, their names, and the units they read in **/ 


(def facts sensor_list 

"(define sensor: Tsubassembly Tsensor ?units)“ 

(def ine - sensor C02 REDUCTION JACKET FLOW 1 CFM) 

(def ine - sensor C02 - REDUCTI0N JACKET - FL0W - 2 CFM) 

(def ine - sensor C02 - REDUCTI0N SYSTEM~EXIT~TEHP K) 

(def ine - sensor C02 - REDUCTI0N REACTION TEMP K) 

) 

/** This section defines the ranges the sensors will be reading if certain problems are occurring. These ar the "rules" 
referred to in the report portion of this paper. **/ 


(def facts parameter_level_def ini tions 


II 

sensor 




parameter 

modifier 

minimum 

lower 

upper 

maximum 

units" 

(level: if 

JACKET FLOW 1 

is 

X 

then 

COOLANT FLOW 1 is 

LO 

0.2 

0.2 

0.4 

0.4 

CFM) 

(Rule 1) 

(level: if 

JACKET - FLOW — 1 

is 

X 

then 

COOLANT - FLOW - 1 is 

NORMAL 

0.4 

0.4 

0.7 

0.7 

CFM) 

(Rule 2) 

(level: if 

JACKET - FLOW - 2 

is 

X 

then 

COOLANT - FLOW - 2 is 

LO 

0.2 

0.2 

0.4 

0.4 

CFM) 

(Rule 3) 

(level: if 

JACKET - FLOW — 2 

is 

X 

then 

COOLANT - FLOW - 2 is 

NORMAL 

0.4 

0.4 

0.7 

0.7 

CFM) 

(Rule 4) 

(level: if 

SYSTEM EXIT TEMP 

is X then CONDENSER TEMP 

is HIGH 

360 

373 

1000 

1000 

K) 

(Rule 5) 

(level: if 

SYSTEM EXIT TEMP 

is X then CONDENSER TEMP 

is NORMAL 

300 

310 

350 

360 

K) 

(Rule 6) 

(level: if 

REACTION TEMP 

is 

X 

then 

EFFICIENCY LEVEL 

is LO 

873 

900 

2000 

2000 

K) 

(Rule 7) 

(level: if 

REACTION - TEMP 

is 

X 

then 

EFFICIENCY - LEVEL 

is NORMAL 

400 

450 

850 

873 

K) 

(Rule 8) 

(level: if 

REACTION - TEMP 

is 

X 

then 

EFFICIENCY - LEVEL 

is NONFUNCTIONAL 

300 

310 

380 400 

K) 

(Rule 9) 


) 


/** This section defines the probable trouble occurring for the sensor indications listed above. **/ 


(def facts condi ti on_def ini ti ons 

"condition: ~ strength" 

(condition: if COOLANT_FLOU_1 is LO then EXIT TEMP is HIGH 1) 

(condition: if COOLANT_FLOU 1 is NORMAL then EXIT~TEMP is OK 1) 

(condition: if C00LANT_FL0W - 2 is LO then EXIT - TEMP is HIGH 1) 

(condition: if COOLANT FLOW 2 is NORMAL then EXIT TEMP is OK 1) 

(condition: if CONDENSER_TEMP is HIGH then CONDENSER is NOT WORKING 1) 

(condition: if CONDENSER TEMP is NORMAL then CONDENSER is WORKING 1) 

(condition: if E F F I C I ENCY_LEVEL is LO then REACTOR is TOO HOT 1) 

(condition: if EFFICIENCY LEVEL is NORMAL then REACTOR is AT RIGHT_TEMP 1) 

(condition: if EFFICIENCY - LEVEL is NOT FUNCTIONAL then REACTOR is TOO_COLD 1) 

) 


/** This section defines the expert system's suggested response to problems occurring in the assembly. **/ 


(deffacts response_def ine 

"?sub_assembly ?mode Tcondition Tstrength" 

(response: if EXIT_TEMP is HIGH then C02 REDUCTION COMMAND is OFF 1) 

(response: if EXIT TEMP is OK then CO2 - REDUCTION - C0MMAND is NORMAL 1) 

(response: if CONDENSER is NOT WORKING then C02_REDUCT I ON COMMAND is OFF 1) 

(response: if CONDENSER is WORKING then C02_REDUCTI0N COMMAND is NORMAL 1) 

(response: if EFFICIENCY LEVEL is LO then C02 REDUCTION COMMAND is OFF 1) 

(response: if EFFICIENCYJ.EVEL is NORMAL then - C02_REDUCTI0N COMMAND is NORMAL 1) 
(response: if EFFICIENCY LEVEL is NONFUNCTIONAL then C02 REDUCTION COMMAND is OFF 1) 
) 


/** This section sends a set of command signals to the system based on the decisions made above. **/ 


(deffacts implementation_commands 

"?sub_assembly ?mode Tcommand" 

(implement: if C02_REDUCTI0N COMMAND is OFF then COMMAND C02_REDUCTI0N to OFF) 
(implement: if C02 REDUCT I ON'COMMAND is NORMAL then COMMAND C02 REDUCTION to NORMAL) 
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APPENDIX 3C 

OXYGEN GENERATION ASSEMBLY KNOWLEDGE BASE 


/** This set of statements defines the sensors, their names, and the units they read in **/ 
(def facts sensor list 


(sensor: 

Tsubassenfaly Tsenor 

?units)" 

(sensor: 

OGA 

H20 FEEDLINE PRESSURE 

PS I 

) 

(sensor: 

OGA 

H20 - PHOTOCELL 

logical) 

(sensor: 

OGA 

H/0 DEFERENTIAL PRESSURE 

PS I 

) 

(sensor: 

OGA 

02 OUTTAKE TEMP - 

F 

) 

(sensor: 

OGA 

CELL VOLTAGE 

volts ) 


) 


/** This section defines the ranges the sensors will be reading if certain problems are occuring. These 
are the "rules" referred to in the report portion of this paper. **/ 


(def facts parameter_level_def i ni t i ons 

" sensor parameter modifier minimum lower upper maximum units" 


(level: 

if 

H20 FEEDLINE PRESSURE is X then MOOULE H20 PRESSURE 

is 

HI 


30 

31 

40 

40 PSI) 

(rule 1) 

(level: 

if 

H20 FEEDLINE PRESSURE is X then MODULE~H20~PRESSURE 

is 

NOMINAL 

19 

20 

30 

31 PSI) 

(rule 2) 

(level: 

if 

H20 FEEDLINE PRESSURE is X then MOOULE~H20~PRESSUR£ 

is 

LOU 


0 

0 

19 

20 PSI ) 

(rule 3) 

(level: 

if 

H20 _ PHOTOCELL is X then H20 SUPPLY LEVEL is HI 




1 

1 

1 

1 BIN) 

(rule 4) 

(level: 

if 

H20~PHOTOCELL is X then H20 - SUPPLY - LEVEL is LOW 




0 

0 

0 

0 BIN) 

(rule 5) 

(level: 

if 

H/0 DEFERENTIAL PRES is X then CROSS CATALYIST PRESSURE 

is 

HI 

3.1 

3.35 

4 

4 PSI) 

(rule 6) 

(level: 

if 

H/0 DEFERENTIAL PRES is X then CROSS~CATALYIST PRESSURE 

is 

NOMINAL 2.83 

2.9 

3.1 

3.35PSI) 

(rule 7) 

(level: 

if 

H/0 DEFERENTIAL~PRES is X then CROSS~CATALYIST~PRESSURE 

is 

LOW 

0 

0 

2.83 

2.9PSI) 

(rule 8) 

(level: 

if 

02 OUTTAKE TEMP is X then OUTTAKE 02 is HOT 




105.8 

110 

125 

125 F) 

(rule 9) 

(level: 

if 

02 OUTTAKE TEMP is X then 0UTTAKE~02 is NOMINAL 




50 

55.4 

105.8 

110 F) 

(rule 10) 

(level: 

if 

02 OUTTAKE~TEMP is X then 0UTTAKE~02 is COLD 




32 

32 

50 

50 F) 

(rule 11) 

(level: 

if 

CELL VOLTAGE is X then CELL VOLTAGE is HI 




1.68 

1.7 

2 

2 V) 

(rule 12) 

(level: 

if 

CELL~VOLTAGE is X then CELL VOLTAGE is NOMINAL 




1.6 

1.64 

1.68 

1.7V) 

(rule 13) 


) 


/** This section defines the probable trouble occuring for the sensor indications listed above. **/ 


(def facts condi t i on_def i ne: 

"condition: 

(condition: if MOOULE H20_PRESSURE is HI then 

(condition: if MODULE~H20 PRESSURE is NOMINAL then 

(condition: if MODULE Jt20~PRESSURE is LOW then 

(condition: if H20 SUPPLY~LEVEL is HI then 

(condition: if H20 SUPPLY LEVEL is LOW then 

(condition: if CROSS CATALYIST PRESSURE is HI then 

(condition: if CROSS~CATAL Y I ST “PRESSURE is NOMINAL 

(condition: if CROSS~CATALYIST PRESSURE is LOW then 

(condition: if OUTTAKE 02 is HOT then 

(condition: if 0UTTAKE~02 is NOMINAL then 

(condition: if 0UTTAKE_02 is COLD then 

(condition: if CELL VOLTAGE is HI then 

(condition: if CELL - VOLTAGE is NOMINAL then 

) 


strength 


H20 DELIVERY PRESSURE is HI 1) 
H20~DELIVERY~PRESSURE is NOMINAL 1) 
H20 DELIVERY PRESSURE is LOW 1) 
WATER TANK is OK T) 
WATER~TANK is EMPTY 1) 
CELL ELECTRODE is BAD 1) 
then"CELL ELECTROOE is GOOD 1) 
CELL ELECTRODE is BAD 1) 
THROTTLE is TOO TIGHT 1) 
THROTTLE is FINE 1) 
THROTTLE is TOO OPEN 1) 
INTAKE H20 is LOST 1) 
INTAKE~H20 IS FINE 1) 


■■ 


/** This section defines the expert systems suggested response to problems occuring in the assembley **/ 


(def facts response_def i ne 

"?sub_assembly ?mode 


(response: 

(response: 

(response: 


f H20 DELIVERY PRESSURE 
f H20 DELIVERY PRESSURE 
f H20 _ 0ELIVERY"PRESSURE 


Tcondition Tstrength" 
is HI then OGA COMMAND is SHUTDOWN 
is NOMINAL then 0GA~C0MMAND is NORMAL 
is LOW then OGA~COMMAND is SHUTDOWN 


(response: 

if 

WATER TANK is OK 

then 

OGA 

COMMAND 

is 

NORMAL 

1) 

(response: 

if 

WATER - TANK is EMPTY 

then 

OGA' 

'COMMAND 

is 

SWITCH 

1) 

(response: 

if 

CELL ELECTROOE is BAD 

then 

oga' 

'command 

is 

SHUTDOWN 

1) 

(response: 

if 

CELL - ELECTRODE is GOOD 

then 

oga' 

’command 

is 

NORMAL 

1) 

(response: 

if 

THROTTLE is TOO TIGHT 

then 

oga command 

is 

SHUTDOWN 

1) 

(response: 

if 

THROTTLE is FINE 

then 

oga 

command 

is 

NORMAL 

1) 

(response: 

if 

THROTTLE is TOO OPEN 

then 

oga 

'command 

is 

SHUTDOWN 

1) 

(response: 

if 

INTAKE H20 is LOST 

then 

oga' 

'command 

is 

SHUTDOWN 

1) 

(response: 

if 

INTAKE~H20 IS FINE 

then 

oga' 

'command 

is 

NORMAL 

1) 


1 ) 

1 ) 

1) 
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/** This section sends a set of command signals to the system based on the decisions made above. **/ 


(def facts implementation_commands 

“ “ response?sub_assembly ?mode 


(implement: 

(implement: 

(implement: 

(implement: 

(implement: 

(implement: 


if OGA COMMAND is NORMAL then 
if OGA'COMMAND is SHUTDOWN then 
if OGA COMMAND is UNPOWERED then 
if OGA~COMMAND is STANDBY then 
if OGA~COMMAND is CURRENT then 
if OGA COMMAND IS SWITCH then 


COMMAND: 

COMMAND: 

COMMAND: 

COMMAND: 

COMMAND: 

COMMAND: 


) 


Tconroand" 

OGA to NORMAL > 
OGA to SHUTDOWN ) 
OGA to UNPOWERED > 
OGA to STANDBY ) 
OGA to CURRENT ) 
OGA to SWITCH ) 
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APPENDIX 3D 

SAVANT.3 EXPERT SYSTEM CODE 


; C02_REM0VAL - KNOWLEDGE BASE 


(deffacts sensor_list 

; "(sensor: ?sub-assembly ?sensor ?units )“ 

(sensor: C02 REMOVAL INLET_C02 ppm ) 

(sensor: C02~REMOVAL C02 ACCUMULATOR g ) 

(sensor: OGA H20_ACCUMULAT0R g ) 

(sensor: OGA CABIN 02 ppm ) 

) 


(deffacts parameterjnodifier_def initions " define triangle membership 


(level: if ’sensor is ?X then Tparameter is ?modifier ?minimum 
(level: if INLET_C02 is X then C02 LEVEL is HI 900 

(level: if INLET_C02 is X then C02 - LEVEL is NORMAL 275 

(level: if INLET_C02 is X then C02 - LEVEL is LOW 0 


functions" 

’lower ?upper ’maximun ?units) 
1100 1000000 1000000 ppm) 

325 900 1100 ppm) 

0 275 325 ppm) 


(level: 

if 

C02 ACCUMULATOR 

is 

X 

then 

C02 

SUPPLY 

is 

HI 

700 

900 

1000 

1000 

9 

) 

(level: 

if 

C02~~ ACCUMULATOR 

is 

X 

then 

C02 

SUPPLY 

is 

OK 

100 

300 

700 

900 

9 

) 

(level: 

if 

C02 ACCUMULATOR 

is 

X 

then 

C02" 

SUPPLY 

is 

LOW 

50 

75 

100 

300 

9 

) 

(level: 

if 

C02~ACCUMULATOR 

is 

X 

then 

C02' 

SUPPLY 

is 

CRITICAL 0 

0 

50 

75 

9 

) 

(level: 

if 

C02~ACCUMULATOR 

is 

X 

then 

C02" 

SUPPLY 

is 

EMPTY 

0 

0 

5 

10 

9 

) 

(level: 

if 

H20 ACCUMULATOR 

is 

X 

then 

H20 

SUPPLY 

is 

HI 

1400 

1800 

2000 

2000 

9 

) 

(level: 

if 

H20~ ACCUMULATOR 

is 

X 

then 

H20 

SUPPLY 

is 

OK 

200 

600 

1400 

1800 

9 

) 

(level: 

if 

H20 _ ACCUMULATOR 

is 

X 

then 

H20" 

SUPPLY 

is 

LOW 

100 

300 

200 

600 

9 

) 

( level : 

if 

H20 ACCUMULATOR 

is 

X 

then 

H20 

’supply 

is 

CRITICAL 0 

0 

100 

300 

9 

) 

(level: 

if 

H20 - ACCUMULAT0R 

is 

X 

then 

H20_SUPPLY 

is 

EMPTY 

0 

0 

5 

10 

9 

) 


(level: 

if 

CABIN 02 

is 

X 

then 02 LEVEL 

is 

HI 

230000 

250000 

1000000 

1000000 

ppm 

) 

(level: 

if 

CABIN 02 

is 

X 

then 02~LEVEL 

is 

OK 

15000 

18000 

23000 

25000 

ppm 

) 

(level: 

if 

CABIN 02 

is 

X 

then 02 LEVEL 

is 

LOW 

12000 

15000 

15000 

18000 

ppm 

) 

(level: 

if 

CABIN 02 

is 

X 

then 02 — LEVEL 

is 

CRITICAL 0 

0 

12000 

15000 

ppm 

) 


(deffacts define_conditions 


(condition: if ’parameter is Tmodifier 
(condition: if C02 LEVEL is HI 
(condition: if 02 LEVEL is CRITICAL 
(condition: if C02_SUPPLY is LOW 

(condition: if C02_LEVEL is NORMAL 

(condition: if C02 LEVEL is LOW 

(condition: if C02~LEVEL is LOW 


then ’condition 
then C02_REMOVAL_DEMAND 
then C02 REMOVAL DEMAND 
or H20 SUPPLY is CRITICAL then C02J(EM0VAL _ DEMAND 

then C02 REMOVAL DEMAND 


and C02 SUPPLY is LOW 


then C02_REM0VAL DEMAND 
then C02 REMOVAL~DEMAND 


is Tmodifier ?strength) 
is HI 1 ) 
is HI 1 ) 
is HI 1 ) 
is NORMAL .8 ) 
is LOW .7 ) 
is NORMAL .9 ) 


(condition: if 02_LEVEL is HI 
(condition: if 02 LEVEL is not HI 


then 02 GENERATION DEMAND is LOW 1 ) 

then 02 GENERATI0N_DEMAND is HI 1 ) 


(condition: if H20_SUPPLY is HI 
(condition: if H20 SUPPLY is not HI 


then C02 REDUCT 1 0N_DEMAND is LOW 1 ) 

then C02~REDUCTI0N_DEMAND is HI 1 ) 


(condition: 

(condition: 

(condition: 

(condition: 


if H20_SUPPLY is EMPTY 
if H20 SUPPLY is CRITICAL 
if C02~ SUPPLY is EMPTY 
if C02~SUPPLY is CRITICAL 


) 


then 0GA_SUB- ASSEMBLY is OUT OF WATER 1 ) 

then OGA SUB-ASSEMBLY is OUT~OF~WATER .2 ) 

then C02 - REDUCTI0N SUB- ASSEMBLY - i s OUT_OF_C02 1 ) 

then C02~REDUCTI0N~SUB- ASSEMBLY is OUT 0F_C02 .2 ) 


(deffacts def i ne^responses 

; (response: Tf ?condition is Tmodifier then ?sub-assembly is Tcommanded Tstrength) 

(response: if C02 REMOVAL DEMAND is HI then C02 REMOVAL is GO TO HIGH REMOVAL_MOOE 1 ) 

(response: if C02 - REM0VAL DEMAND is NORMAL then C02 - REM0VAL is GO - TO - POUER EFFICIENCY_MODE 1 ) 

(response: if C02_REMOVAL j)EMAND is LOW then C02~REMOVAL is TURN_OFF - 1 ) 

(response: if 02 GENERATION DEMAND is LOW then OGA is TURN OFF 1 ) 

(response: if 02~GEN E RAT I ON - D EMAND is HI then OGA is TURN~0N 1 ) 

(response: if C02_REDUCT I ON - DEMAND is LOW then C02 REDUCTION is TURN~OFF 1 ) 

(response: if CO2_REDUCTION“0EMAND is HI then C02 - REDUCTI0N is TURNJ3N 1 ) 

(response: if OGA SUB-ASSEMBLY is OUT 0F_WATER then OGA is TURN_OFF 1 ) 

(response: if C02~REDUCTION_SUB-ASSEMBLY is 0UT_0F_C02 then C02_REDUCTION is TURN_0FF 1 ) 
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) 

(def facts command-implementation 


(implement: 

if ?sub-assembly 


is ’commanded 

then WRITE: ?sub-assembly to Tstring 

) 

(implement: 

if 

C02 REMOVAL 

is 

GO TO HIGH REMOVAL MODE * 

then WRITE: C02 REMOVAL 

to 

HIGH POWER 

) 

(implement: 

if 

C02 — REMOVAL 

is 

GO TO POWER EFFICIENCY MODE 

then WRITE: C02 - REM0VAL 

to 

ECONOMY 

) 

(implement: 

if 

C02~ REMOVAL 

is 

TURN OFF 

then WRITE: C02 - REM0VAL 

to 

OFF 

) 

(implement: 

if 

OGA 

is 

TURN - OFF 

then WRITE: OGA 

to 

OFF 

) 

(implement: 

if 

OGA 

is 

TURN - 0N 

then WRITE: OGA 

to 

ON 

) 

(implement: 

if 

C02 REDUCTION 

is 

turn“off 

then WRITE: C02 REDUCTION 

to 

OFF 

) 

(implement: 

if 

C02~ REDUCTION 

is 

turn“on 

then WRITE: C02~REDUCTI0N 

to 

ON 

) 


) 


SAVANT. 3 


; "INFERENCE" ENGINE (shell above CLIPS) 

(defrule evaluate_parameters "(direct reading of parameters)" 

(sensor: ?sub-assembly1 Tsensor ?units ) 

(sub-assembly: ?sub-assembly2 ?file_id ) 

(time ?time ?time ?iterations) 

(sensor ?sub-assembly1 Tsensor ?value ?units ?time_stamp) 

(test (= ?time ?time_stamp)) 

=> 

(assert (?sensor is ?value ?units =(gensym) )) 

) 

(defrule determine_parameter_modif iers "evaluate membership function to convert value to weighted level" 
(level: if ?sensor is ?X then ?parameter is ?modifier ?minimum ?lower ?upper ?maximum ?units) 
(Tsensor is ?value ?units ?tag ) 

; (test (and (>= ?value ?minimum) (< ?value ?maximum))) 

=> 

(bind Tweight 0) 

(if (and (> ?value ?minimum) (< ?value ?lower )) 

then (bind ?weight (/ (- ?value Tminimun ) (- ?lower ?minimum )))) 

(if (and (>= ?value ?lower ) (<= ?value ?upper)) 
then (bind ?weight 1)) 

(if (and (> ?value ?upper ) (< ?value Tmaximum)) 

then (bind ?weight (/ (- ?value Tmaximum ) (- ?upper ’maximun )))) 

(bind ’confidence ’weight) 

(assert (?parameter is Tmodifier ?confidence =(gensym) )) 

) 

(defrule determine_conditions 

(condition: if $?parameter_statement then $?asserted_condition_statement Tstrength) 
($?parameter_statement ?conf idence ?tag ) 

=> 

(bind Tconfidence (* ?confidence Tstrength)) 

(assert ($?asserted_condition statement Tconfidence =(gensym) )) 

) 

(defrule determine_response 

(response: if $?condition_statement then $?asserted_response_statement Tstrength) 
($?condition_statement ?conf idence ?tag ) 

=> 

(bind Tconfidence (* Vconfidence ’strength)) 

(assert ($?asserted_response_statement Tconfidence =(gensym) )) 

) 

.-mutual renforcement? 

; SAVANT. 3 


; RESOLVE LOGICAL OPERATIONS 

(defrule logical_and "VW" 

(declare (salience -5)) 

?rm <- ($?left S?left statement and $?right_statement $?right) 
($?left_statement ?Tconfidence ?ltag ) 

(test (nunberp ? I confidence )) 

($?right_statement ?rconf idence Trtag ) 

(test (numberp ?rconf idence )) 

=> 

; v 

(bind ?conf i dence (* ?lconf idence ?rconf idence )) 

(retract ?rm) 

(assert (S’left $’left_statement AND $’right_statement STright)) 
(assert ($?left_statement~AND $?right_statement - ?conf idence =(gensym) )) 
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) 


(defrule logical_or "W" 

(declare (salience -5)) 

?rm <- ($?left $?left_statement or $?right_statement $?right> 
($?left_statement ?Tconfidence ?ltag ) 

(test (nunberp ?lconfidence )) 

($?right_statement ?rconfidence ?rtag ) 

(test (nunberp ?rconfidence )) 


(bind ?conf idence (min (+ ?lconfidence ?rconfidence ) 1)) 

(retract ?rm) 

(assert (S’left $?left_statement OR $?right_statement $?right)) 
(assert ($?left statement OR $?right statement - ?conf idence =(gensym) )) 
) 

(defrule logical_not 

(declare (salience -5)) 

?rm <- ($?left ?noun is_not ?modifier $?right> 

(?noun is ?modifier ’confidence ?tag ) 

(test (nunberp ’confidence )) 


(bind ?not_conf idence (-1 ?conf idence )) 

(retract ?rm) 

(assert ($?left ?noun 1S_N0T ?modifier $?right)) 

(assert (?noun IS NOT ?modifier ?not_conf idence =(gensym) )) 

) 

(defrule mutual_renforcement "'OR' confidence values of idential facts 
[those with identical nouns and modifiers]. 

Form a single fact from the two." 

?rm1 <- (?noun1 is ?modifier1 ?confidence1 ?tag1 ) 

?rm2 <- (?noun2 is ?modifier2 ?confidence2 ?tag2 ) 

(test (and (eq ?noun1 ?noun2) 

(eq ’modifierl ?modifier2 ) 

(neq ’tagl ?tag2 ))) 

=> 

(bind ’confidence (min (+ ?confidence1 ’conf idence2) 1)) 

(retract ?rm1 ?rm2) 

(assert (?noun1 is ?modifier1 ?confidence =(gensym) )) 

) 

; SAVANT. 3 

; RESPONSE SELECTION AND IMPLEMENTATION 

(defrule select_response1 

(declare (salience -10)) 

(response: if $?condition_statement then ?sub- assent) l y is ?command ?strength) 
’rml <- (?sub- assembly is ?command1 Tconfidencel ?tag1 ) 

?rm2 <- (?sub- assembly is ?command2 ?confidence2 ?tag2 ) 

(test (neq ’commandl ?command2 )) 

(test (neq ’tagl ?tag2)) 

=> 

(if (< Tconfidencel ?conf idence2) 
then (retract ’rml) 
else (retract ?rm2)) 

) 

.•threshold? 

(defrule implement_response 
(declare (salience -100)) 

(implement: if ?sub-assembly1 is ?mode then WRITE: ?sub-assembly to ’string ) 
(?sub-assembly1 is ?mode ?confidence ’tag ) 

(sub-assembly: ?sub-assembly2 ?file_id) 

(time ’time ?time ?iterations) 

=> 

(bind ?file id. ext (str_cat ?file_id ".CMD")) 

(bind $?comman3_fact (mv-append command ’sub-assemblyl to ?string ?time)) 
(printout ?file~id.ext "(" $?comnand_fact ")" crlf ) 

(printout t “ "(" $?command~fact ")" crlf ) 

) 

; SAVANT. 3 


EXPLANATION FACILITY 
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(defrule echo_parameters 

(sensor: ’sub-assemblyl ?sensor ?units ) 

(?sensor is ?value ’units ?tag1 ) 

(level: if ’sensor is ?X then ’parameter is ?modifier $’ ?units) 

(?parameter ?link ?modifier ?confidence ?tag2 ) 

(test (> ?confidence 0 )) 

=> 

(printout t ?sensor " is 11 ?value 11 " ?units 11 => " ?parameter " " ?link " 11 ?modifier " 11 ?confidence crlf) 

) 

(defrule evaluate_parameters "(direct reading of parameters)" 

(sensor: ?sub-assembly1 ?sensor ?units ) 

(sub-assembly: ?sub-assembly2 ?file_id ) 

(time ?time ’time ?iterations) 

(sensor ?sub-assembly1 ?sensor ?value ?units ?time_stamp) 

(test (= ’time ?time_stamp)) 

=> 

(assert (?sensor is ?value ?units =(gensym) )) 

) 

(defrule determi ne_parameter_modi f iers "evaluate membership function to convert value to weighted level" 

(level: if ?sensor is ?X then ?parameter is ?modifier ’minimun ?lower ?upper ’maximum ?units) 

(?sensor is ?value ?units ?tag ) 

; (test (and (>= ’value ?minimun) (< ?value ?maximum))) 

=> 

(bind ?weight 0) 

(if (and (>= ?value ?minimum) (< ?value ?lower )) 

then (bind ’weight (/ (- ?lower ?value) (- ?lower ?minimun)))) 

(if (and (>= ?value ?lower ) (< ?value ?upper)) 
then (bind ?weight 1)) 

(if (and (>= ?value ?upper ) (< ?value ?maximum)) 

then (bind ?weight (/ (- ’maximum ?value) (- ?maximum ’upper)))} 

(bind ?confidence ?weight) 

(assert (?parameter is ’modifier ?confidence =(gensym) )) 

) 


(def ru l e start_expl anat i on_of _response 
(declare (salience -1000)) 

?rm <- (?response is ?mode ?confidence ?tag ) 

(response: if $?condition_statement then ?response is ’mode ’strength) 

=> 

(retract ?rm) 

(assert (?response is ?mode ?confidence )) 

(assert (explain_response: ?response ?mode)) 

(printout t crlf crlf crlf crlf crlf crlf crlf crlf crlf crlf crlf "EXPLANATION:" crlf) 
(printout t "SAVANT commands the " Tresponse " sub-assembly to " ?mode "." crlf crlf) 

(printout t "This is because SAVANT has determined that the following conditions are " crlf) 
(printout t "likely to exist onboard, and SAVANT has been told that the proper response " crlf) 
(printout t "to these conditions is to command the " Tresponse 11 sub-assembly " crlf) 

(printout t "to " ?mode "." crlf) 

) 

(defrule explain_conditions 
(declare (salience -10)) 

(explain_response: ? response ?mode) 

(response: if Tcondition is Tmodifier then Tresponse is ?mode Tstrength) 

( Tcondition is Tmodifier ?confidence ?tag ) 

(test (> ?confidence 0)) 

=> 

(assert (explain_condition: Tcondition ?modifier )) 

(printout t crlf "CONDITION: " ’condition " is likely to be " Tmodifier crlf) 

(printout t " This condition is likely because: " crlf) 

) 

(defrule list_evidence 

(explain_response: ?response ?mode) 

(response: if (condition is Tcmodifier then ?response is ?mode ?rstrength) 

(explain_condition: ?condition ?cmodifier ) 

(condition: if $?parameter_statement then Tcondition is Tcmodifier ?cstrength) 
($?parameter_statement ?pconf idence ?tag ) 

(test (> ?pconf idence 0)) 

=> 

(printout t " > " $?parameter statement " ( " ?pconfidence " ) " crlf) 

) 

(defrule comp l ete_expl anat i on_of _cond i t i on 
(declare (salience -10))~ 

(explain_response: ? response ’mode) 
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(? response is ?mode ’confidence ) 


(printout t crlf "The cumulative confidence of this conclusion is " ?confidence crlf) 
(printout t "Continue (Y/N)? ") 

(bind ?continue (read)) 

(if (or (eq ?continue n) (eq ’continue no) (eq ?continue M) (eq ’continue NO) ) 
then (halt)) 

) 

SAVANT. 3 


; FILE I/O OVERHEAD 

( def ru l e open_command_f i l e 

(declare (salience 100)) 

(sub-assembly: ?sub- assent l y_name ?file_id) 

=> 

(bind ?file_id.ext (str_cat ?file_id ".CMD")) 
(open ’file idText ?file_id.ext “w") 

) 

(def rule close_command_f i le 
(declare (salience -100)) 

(sub-assembly: ?sub-assembly_name ?file_id) 

=> 

(bind ?file_id.ext (str_cat ?file_id ".CMD")) 
(close ’file id. ext) 

) 

; SAVANT. 3 


(deffacts sub-assembly_l/0_f i le_def initions 

; (sub-assembly: ?sub-assembly ’filename ) 

(sub-assembly: C02_REM0VAL "C02REM0V" ) 

) 

(deffacts conf igure_program 

" sub-assembly file id " 

TO) 

; (load ".CMD") 

(load ".SEN") 

(extract time 0) 

) 

(def rule load_facts_f i le 

(declare (salience 100)) 

(sub-assembly: ?sub-assembly_name ?file_id) 
(load ?ext) 

=> 

(bind ?file_id.ext (str_cat ?file_id ?ext)) 
(load-facts ?fTle id. ext) 

) 

; SAVANT. 3 


; TIME STAMP EXTRACTION AND INSERTION 

(def rule extract_time "what time is it?" 

(declare (salience 100)) 

?rm <- (extract_time ’time) 

(sensor C02_REMOVAL CLOCK ?iterations min ?time_stamp ) 
(test (< ?time ?time_stamp )) 

=> 

(if (!= ?time_stamp ?time) 

then (bind ?time ?time_stamp ) 

(retract ?rm) 

(assert (extract time ?time)) 

)) 

(defrule time_extracted "time is:" 

(dec l are~( salience -100)) 

?rm <- (extract_time ?time) 

(sensor C02_REMOVAL CLOCK ?iterations min ?time) 

=> 

(retract ?rm) 

(assert (time ?time ?time ?iterations)) 

) 


149 





SAVANT. 3 


; EXIT 

(defrule exit 

(declare (salience -1000)) 

(time ’time ?time ?iterations) 

=> 

(save- facts "facts") 

(if (>= ?time 60 ) 
then (halt) 

; block out following line with a semicolon to keep system from looping 

; else (batch "model.bat") 

)) 
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